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In  this  work,  we  investigate  issues  on  the  application  of  a physical  layer  technique, 
sequence  optimization,  in  code  division  multiple-access  (CDMA)  communication  systems. 
The  objective  of  sequence  optimization  is  to  construct  a set  of  optimal  sequences  so  that 
users  in  the  system  are  supported  with  target  signal-to-interference  ratios  (SIRs)  using  a 
minimum  amount  of  total  transmission  power.  We  propose  two  joint  transmitter-receiver 
adaptation  schemes  that  make  use  of  the  optimal  sequences  for  CDMA  systems  over  the  ad- 
ditive white  Gaussian  noise  (AWGN)  channel.  In  both  schemes,  training  symbols  are  sent 
in  blocks  so  that  spreading  sequences  and  minimum  mean  squared  error  (MMSE)  receiver 
weights  are  adapted  simultaneously.  In  the  first  scheme,  spreading  sequences  are  updated 
at  the  receivers  in  a distributed  manner,  and  are  fed  back  to  the  base-station  to  modulate 
the  next  block  of  data.  In  the  second  scheme,  spreading  sequences  are  updated  at  the  base- 
station  in  a centralized  manner  using  a minimum  amount  of  feedback  information  from  the 
receivers.  Our  second  objective  concerns  sequence  optimization  in  CDMA  systems  over 
fading  channels.  We  propose  two  sequence  allocation  schemes,  a minimum-ETSC  scheme 
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for  MMSE  signal  reception  and  a minimum-TSC  scheme  for  matched-filter  signal  recep- 
tion, to  use  in  a CDMA  forward  link  over  slow  flat  fading.  The  MMSE  receiver  weights 
of  the  minimum-ETSC  sequences  turn  out  to  be  the  same  as  those  of  the  matched-filter  re- 
ceiver. There  is  strong  indication  that  the  minimum-ETSC  scheme,  which  has  better  power 
efficiency  than  the  minimum-TSC  scheme,  may  be  the  optimal  scheme  when  the  matched 
filter  is  used  for  signal  reception.  In  addition,  Lagrangian  multiplier  searching  results  indi- 
cate that  the  minimum-ETSC  scheme  is  effective  and  consumes  only  slightly  more  power 
than  the  minimum  power  searched  under  most  channel  conditions  when  the  MMSE  re- 
ceiver is  used.  Adaptive  structure  and  algorithms  in  combination  with  channel  estimation 
are  also  presented  to  facilitate  the  implementation  of  the  minimum-ETSC  sequences  in 
practical  systems.  We  also  investigate  sequence  optimization  in  CDMA  systems  with  mul- 
tipath. We  solve  the  optimization  problem  for  the  case  of  point-to-point  transmission,  in 
which  all  the  receivers  are  co-located  and  undergo  the  same  multipath  channel.  Finally, 
we  propose  a type-I  hybrid  ARQ  protocol  to  make  use  of  the  optimal  sequences  to  allow 
simultaneous  transmission  of  multiple  packets  over  a CDMA  link.  The  transmitted  packets 
are  classified  into  two  classes,  namely  newly  arrived  packets  and  retransmission  packets. 
A resource  allocation  algorithm  is  responsible  for  determining  the  numbers  and  the  target 
SIRs  for  each  class  of  packets  to  be  transmitted  in  each  time  slot  so  that  the  link  throughput 
is  maximized  under  a delay  constraint.  Numerical  results  obtained  from  both  the  Markov 
analysis  and  Monte  Carlo  simulations  are  used  to  evaluate  the  performance  of  the  proposed 
protocol. 
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CHAPTER  1 
INTRODUCTION 


Ever  since  analog  cellular  communication  systems  started  to  reach  their  capacity 
limitation  in  1980s,  development  of  digital  cellular  communication  systems  has  been  car- 
ried out  to  increase  system  capacity.  Frequency  division  multiple  access  (FDMA),  time 
division  multiple  access  (TDMA),  and  code  division  multiple  access  (CDMA)  are  three 
basic  multiple  access  schemes  employed  to  share  the  available  bandwidth  in  cellular  com- 
munication systems.  In  theory,  these  three  multiple  access  schemes  provide  the  same  ca- 
pacity, no  matter  whether  the  spectrum  is  divided  into  frequency  channels,  time  slots,  or 
codes.  However,  in  digital  cellular  communication  systems,  CDMA  exhibits  great  advan- 
tages over  FDMA  and  TDMA.  In  FDMA  or  TDMA,  each  frequency  channel  or  each  time 
slot  is  assigned  to  one  call,  and  no  other  call  can  share  the  same  channel  or  slot  during 
the  call  period.  In  CDMA,  additional  code  sequences  can  be  squeezed  into  the  same  radio 
channel,  which  is  a property  unique  to  CDMA.  It  is  recognized  that  CDMA  capacity  is  in- 
terference limited,  whereas  FDMA  and  TDMA  capacities  are  primarily  bandwidth  limited 
[1-3].  Hence,  any  reduction  of  interference  translates  directly  into  an  increase  of  capacity 
in  CDMA  systems.  It  is  known  that  voice  signals  are  intermittent  with  a duty  factor  of 
approximately  3/8.  Suppression  of  transmissions  during  the  silent  periods  of  each  speaker 
can  increase  capacity  by  an  amount  inversely  proportional  to  this  factor.  Similarly,  spatial 
isolation  through  the  employment  of  multi-beamed  or  multi-sectored  antennas  or  antenna 
arrays,  which  reduces  interference,  also  provides  a proportional  increase  in  capacity.  In 
CDMA  cellular  systems  all  the  cells  reuse  the  same  (entire)  broadband  spectrum  simulta- 
neously, and  isolation  among  cells  is  provided  by  path  loss,  which  typically  increases  with 
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the  fourth  power  of  the  distance  for  terrestrial  UHF  propagation.  Consequently,  employ- 
ment of  CDMA  increases  the  capacity  by  the  frequency  reuse  factor,  as  opposed  to  FDMA 
and  TDMA,  where  different  frequency  bands  have  to  be  allocated  to  contiguous  cells.  Be- 
cause of  its  attractive  properties  for  the  wireless  medium,  there  has  been  substantial  interest 
in  the  CDMA  technology  in  recent  decades.  For  instance,  in  the  US  digital  cellular  mo- 
bile standard  IS-95  [4],  DS-CDMA  is  selected  as  the  technology  to  support  the  second 
generation  mobile  cellular  services.  Also,  Wideband  CDMA  (W-CDMA)  is  recommended 
to  use  in  the  third  generation  mobile  communication  system,  called  International  Mobile 
Telecommunications-2000  (IMT-2000). 

In  this  work,  we  consider  the  direct-sequence  code  division  multiple  access  (DS- 
CDMA)  technique,  which  provides  multiple  access  by  direct-sequence  spread  spectrum 
signaling  [5],  in  wireless  communication  systems.  In  order  to  increase  the  system  capacity, 
advanced  processing  techniques  need  to  be  developed  in  both  the  physical  layer  and  upper 
layers.  We  present  research  results  in  both  the  physical  layer  and  the  data  link  layer  of  a 
DS-CDMA  system  in  this  work. 

1.1  Suppression  of  MAI  in  DS-CDMA  Systems 

The  important  feature  of  FDMA  and  TDMA  is  that  the  various  users  are  operating 
in  separate,  non-interfering  channels.  However,  in  DS-CDMA  systems,  the  information 
bearing  signal  for  each  user  is  spread  over  a wider  bandwidth  by  means  of  a spreading 
waveform  unique  to  that  user,  and  all  the  users  share  the  same  broad  bandwidth  at  the 
same  time.  Hence,  the  receiver  will  actually  observe  the  sum  of  all  users’  signals  embed- 
ded in  additive  background  noise.  Multiple  access  interference  (MAI),  which  refers  to  the 
interference  between  direct-sequence  users  due  to  many  simultaneous  users  in  the  same 
frequency  band,  is  a factor  that  limits  the  capacity  and  performance  of  DS-CDMA  sys- 
tems. Unique  to  CDMA  systems,  the  MAI  caused  by  any  single  user  is  generally  small. 
However,  MAI  becomes  substantial  and  can  cause  severe  channel  imperfections  as  the 
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power  of  interference  increases.  If  the  spreading  waveforms  for  different  users  as  seen 
at  the  receiver  are  orthogonal,  MAI  could  be  eliminated.  However,  in  practice,  signals 
from  different  users  arrive  at  the  receiver  at  different  delays,  and  it  is  impossible  to  de- 
sign spreading  waveforms  that  maintain  orthogonality  over  all  the  possible  relative  delays. 
Hence,  more  effort  focuses  on  designing  waveforms  that  are  nearly  orthogonal,  that  is,  to 
design  waveforms  that  have  low  cross-correlation.  Conventional  system  design  exploits 
this  property  of  the  waveform  and  ignores  MAI  in  receiver  design.  For  instance,  the  IS-95 
[4]  standard  uses  matched  filter  reception.  The  matched  filter  receiver  is  relatively  easy  to 
implement.  However,  it  is  optimal  only  if  there  is  no  MAI.  When  the  signal  from  a distant 
mobile  unit  is  very  weak,  it  might  be  overwhelmed  by  the  contribution  of  a nearby  strong 
interfering  signal  at  the  matched  filter  output.  The  strong  interference  makes  it  impossible 
for  a reliable  reception  of  the  desired  weak  signal.  This  problem  is  usually  addressed  as 
the  near-far  problem  in  wireless  CDMA  systems.  Two  general  approaches  have  been  stud- 
ied extensively  to  overcome  the  near-far  problem  and  improve  the  performance  of  CDMA 
systems:  transmitter  power  control  and  multi-user  detection. 

Transmitter  adaptation  in  the  form  of  power  control  has  been  applied  in  current 
CDMA  systems  [2,  4],  and  is  considered  indispensable  and  mandatory  for  successful  signal 
transmission/reception  in  the  CDMA  reverse  link  [6-10].  In  the  IS-95  standard  [4],  the 
mobiles  adjust  their  power  through  two  methods.  The  first  method,  known  as  open  loop 
power  control,  is  for  a mobile  to  adjust  its  transmission  power  to  be  inversely  proportional 
to  the  power  level  it  receives  from  the  base-station.  The  second  method,  known  as  closed 
loop  power  control,  is  that  the  base-station  periodically  sends  information  to  the  mobiles  in 
its  coverage  area  instructing  them  to  adjust  their  transmission  power  so  that  signals  from  all 
users  arrive  at  the  base-station  at  approximately  the  same  power  level  and  therefore  no  user 
is  substantially  disadvantaged  relative  to  the  other  users.  In  IS-95  systems,  forward-link 
power  control  is  far  less  powerful  than  reverse-link  power  control.  Although  transmitter 
power  control  in  the  CDMA  forward  link  is  generally  thought  not  as  essential  as  in  the 
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reverse  link,  it  has  been  realized  that  an  inefficient  forward-link  power  control  scheme 
limits  the  forward-link  capacity  significantly  [1 1-13].  Several  forward-link  power  control 
approaches  have  been  suggested  [14-18], 

Multi-user  detection,  which  aims  to  develop  near-far  resistant  receivers,  is  another 
approach  to  overcome  the  near-far  problem  as  well  as  to  eliminate  or  alleviate  the  neces- 
sity of  stringent  power  control.  In  multi-user  detection,  code  and  timing  (and  possibly 
amplitude  and  phase)  information  of  multiple  users  is  jointly  used  to  better  detect  each 
individual  user’s  signal.  Verdu  proposes  the  optimum  multi-user  receiver  [19],  which  is 
in  fact  a maximum  likelihood  sequence  receiver.  The  asymptotic  efficiency  of  such  a re- 
ceiver for  DS-CDMA  systems  over  additive  white  Gaussian  noise  (AWGN)  channels  is 
thoroughly  analyzed  [20].  The  asymptotic  efficiency  is  defined  as  the  loss  relative  to  the 
conventional  matched  filter  receiver  operating  in  the  presence  of  only  AWGN  in  the  limit 
as  the  signal-to-noise  ratio  (SNR)  goes  to  infinity.  The  receiver  proposed  by  Verdu  out- 
performs the  matched  filter  receiver,  at  the  expense  of  a significant  increase  in  complexity. 
Specifically,  the  computational  complexity  of  the  optimum  multi-user  receiver  [19]  grows 
exponentially  with  the  number  of  users.  The  optimum  multi-user  receiver  is  too  complex 
for  practical  DS-CDMA  systems.  Therefore,  over  the  last  decade  or  so,  much  research  has 
focused  on  devising  sub-optimal  multi-user  receivers,  which  are  less  complex  and  more 
feasible  to  implement. 

Simpler  linear  multi-user  receivers,  which  have  computational  complexity  that  is 
linearly  proportional  to  the  number  of  users,  have  been  proposed  by  Lupas  and  Verdu 
[21,  22]  to  reject  MAI.  The  major  working  principle  underlying  these  receivers  is  to  ap- 
ply the  inverse  of  the  correlation  matrix  of  the  conventional  receiver  output  to  project  the 
spreading  sequence  of  the  desired  user  onto  the  space  orthogonal  to  the  spreading  se- 
quences of  interferers.  This  decorrelating  detector  completely  eliminates  the  MAI,  and 
provides  a decent  asymptotic  efficiency  with  a very  simple  structure.  A disadvantage  of 
this  receiver  is  that  it  causes  noise  enhancement  (similar  to  a zero-forcing  equalizer).  The 
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power  associated  with  the  noise  term  at  the  output  of  the  decorrelating  detector  is  always 
greater  than  or  equal  to  the  power  associated  with  the  noise  term  at  the  output  of  the  con- 
ventional receiver.  A more  significant  disadvantage  of  the  decorrelating  detector  is  the 
difficulty  to  perform  the  computations  needed  to  invert  the  correlation  matrix  in  real  time. 
For  synchronous  systems,  the  problem  is  somewhat  simplified.  We  can  decorrelate  one 
bit  at  a time  and  the  size  of  the  correlation  matrix  is  small.  For  asynchronous  systems, 
however,  the  size  of  the  correlation  matrix  depends  on  the  message  length  and  channel 
memory.  The  decorrelating  detector  has  received  much  attention,  and  there  have  been  nu- 
merous sub-optimal  approaches  [23-26]  in  the  literature  to  implement  the  decorrelating 
detector. 

The  multi-user  receivers  for  detection  of  a signal  in  the  presence  of  noise  and  MAI 
discussed  above  require  the  information  of  spreading  codes,  timing,  and  phases  of  other 
users.  This  information  may  be  available  at  the  base-station,  but  is  not  likely  to  be  avail- 
able at  the  mobile  units.  In  addition,  these  receivers  may  not  be  feasible  in  a multipath 
environment  where  it  would  be  difficult  to  obtain  the  knowledge  of  timing  of  the  delayed 
paths  of  all  the  received  signals.  If  the  multipath  is  ignored,  the  near-far  resistance  of  these 
receivers  would  be  compromised,  resulting  in  a reduction  in  system  capacity.  One  of  the 
most  commonly  used  antimultipath  techniques  for  DS-CDMA  systems  is  the  RAKE  re- 
ceiver [5],  Pateros  and  Saulnier  [27]  have  shown  that  an  adaptive-equalizer-based  receiver 
can  obtain  performance  very  close  to  that  of  an  ideal  RAKE  receiver,  without  the  knowl- 
edge of  the  delays  of  the  multipath  components.  The  scheme  proposed  [27]  consists  of 
sampling  the  channel  output  at  the  chip  rate,  and  using  an  N-tap  finite  impulse  response 
(FIR)  filter  to  minimize  the  mean  squared  error  (MSE)  between  the  transmitted  and  de- 
tected symbol,  where  N is  the  processing  gain.  The  scheme  makes  use  of  the  fact  that  the 
minimum  mean  squared  error  (MMSE)  linear  detector  just  described  can  be  implemented 
as  a fractionally  spaced  tapped  delay  line  (TDL).  In  the  absence  of  MAI,  the  N-tap  de- 
tector reduces  to  the  conventional  matched  filter  receiver.  The  relative  simplicity  of  this 
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receiver  makes  it  feasible  not  only  at  the  base-station,  but  also  at  the  mobile  units.  Miller 
[28]  published  a thorough  analysis  of  the  MAI  rejection  capability  and  system  capacity  im- 
provement of  this  receiver.  Similar  receiver  structures  based  on  the  MMSE  criterion  also 
have  been  developed  [29-32],  which  have  shown  that  this  type  of  receiver  structure  can 
provide  significant  performance  improvement  in  the  presence  of  multiple  users.  A blind 
adaptive  linear  receiver,  which  eliminates  the  need  for  training  sequences  as  required  in 
the  algorithms  of  Madhow  [30]  and  Rapajic  [31],  has  been  proposed  by  Honig  [33].  The 
MMSE  receiver  generally  provides  a better  error  probability  performance  than  the  decor- 
relating  detector  because  it  takes  the  background  noise  into  account.  As  the  background 
noise  goes  to  zero,  the  MMSE  receiver  converges  in  performance  to  the  decorrelating  de- 
tector. Performances  of  several  important  linear  receivers  including  the  matched  filter, 
MMSE,  and  decorrelator  receivers  are  analyzed  in  a reverse-link  single-cell  “large”  sys- 
tem with  random  sequences  by  Tse  et  al.  [34,  35].  Not  surprisingly,  the  MMSE  receiver  is 
the  most  efficient  one  in  terms  of  maximizing  the  user  capacity  '.  As  has  been  generally 
realized,  one  appealing  property  of  the  MMSE  receiver  is  its  relative  simplicity  from  the 
point  of  view  of  implementation  compared  with  most  other  multi-user  receivers. 

Another  important  group  of  receivers  can  be  classified  as  subtractive  interference 
cancellation  receivers  [36-38],  The  basic  principle  underlying  these  receivers  is  the  cre- 
ation of  separate  estimates  of  the  MAI  contributed  by  each  user  at  the  receiver  so  that  part 
or  all  of  the  MAI  seen  by  the  desired  user  can  be  subtracted  out  from  the  received  signal. 
These  receivers  are  similar  to  decision  feedback  equalizers  used  to  combat  inter-symbol 
interference  (ISI).  In  feedback  equalization,  decisions  of  previously  detected  symbols  are 
fed  back  in  order  to  cancel  part  of  the  ISI.  Thus,  a number  of  these  types  of  multi-user 
receivers  are  also  referred  to  as  decision-feedback  detectors  [39,  40], 

1 User  capacity  refers  to  the  maximum  number  of  users  that  can  be  supported  with  the  desired  QoS 
requirement. 


7 


1.2  Joint  Transmitter-Receiver  Optimization 

Almost  all  the  MAI  suppression  approaches  described  above,  except  adaptive  power 
control,  focus  on  the  receiver  design.  Multi-user  detection  aims  to  develop  near-far  resis- 
tant receivers.  However,  in  a practical  CDMA  system,  multi-user  receivers  are  envisioned 
to  detect  signals  from  a subset  of  the  users  in  the  system.  For  example,  in  the  reverse  link  of 
a cellular  system,  only  signals  from  the  users  within  the  same  cell  can  be  jointly  detected, 
while  signals  from  the  users  other  than  the  desired  cell  are  treated  as  interference.  Hence, 
it  is  conceivable  that  joint  transmitter  power  control  and  multi-user  detection  offers  a per- 
formance advantage  over  power  control  or  multi-user  detection  alone.  Several  transmitter 
power  control  algorithms  in  combination  with  MMSE  signal  reception  have  been  proposed 
[41^4-3].  While  the  combination  of  adaptive  power  control  and  MMSE  signal  reception  is 
useful  to  guarantee  a certain  SIR  at  the  receiver,  it  does  not  exploit  the  full  potential  of 
the  transmitter  design,  which  can  serve  to  pre-equalize  the  channel  and  avoid  the  MAI.  A 
number  of  results  regarding  joint  transmitter  and  receiver  optimization  are  available  in  the 
literature.  One  of  the  approaches  [44]  is  to  precode  all  the  transmitted  signals  by  a linear 
transformation  before  transmitting  to  minimize  MAI.  The  precoding  transform  and  the  lin- 
ear receivers  are  chosen  jointly  to  minimize  the  sum  MSE  at  the  receivers.  This  approach 
can  only  be  applied  to  the  forward  link  because  of  its  centralized  nature.  An  early  work  of 
Rapajic  [45]  suggests  to  iteratively  replace  the  spreading  sequence  of  a user  by  the  weight 
vector  obtained  at  the  MMSE  receiver.  The  potential  of  MAI  suppression  through  joint 
transmitter-receiver  adaptation  is  discussed  and  it  is  shown  that  the  resulting  MSE  at  the 
receiver  can  be  reduced  provided  that  the  transmissions  of  the  other  users  are  fixed.  It  is 
argued  by  Rapajic  and  Vucetic  [45]  and  Wong  [46]  that  the  SIR  achieved  by  the  MMSE 
receiver  depends  on  the  choice  of  the  spreading  sequences.  If  the  spreading  sequences  are 
chosen  or  adapted  suitably  together  with  adaptive  power  control,  the  MAI  level  as  seen 
by  the  MMSE  receiver  can  be  further  suppressed  and  hence  the  performance  of  the  sys- 
tem is  improved.  In  synchronous  DS-CDMA  systems,  orthogonal  spreading  sequences 
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are  generally  employed  to  avoid  MAI.  The  number  of  orthogonal  sequences  is  limited  by 
the  value  of  the  spreading  gain  and  thus  makes  itself  a bottleneck  limiting  the  number  of 
admissible  users.  As  we  mentioned  above,  a unique  property  of  CDMA  technology  is 
that  additional  sequences  can  be  squeezed  into  the  system.  Making  use  of  this  property, 
capacity  restriction  can  be  avoided  via  the  employment  of  additional  sequences  that  are 
not  orthogonal  to  each  other.  Introduction  of  non-orthogonal  spreading  sequences  causes 
MAI,  and  makes  power  control  and  sequence  allocation  essential  to  guarantee  the  QoS  of 
the  system.  In  recent  years,  much  attention  has  been  paid  to  an  approach  of  spreading 
sequence  optimization  in  conjunction  with  MMSE  signal  reception. 

An  early  work  of  Massy  [47]  states  that  a set  of  equal-energy  Welch-bound-equality 
(WBE)  sequences,  which  meets  the  Welch’s  lower  bound  [48]  on  the  total  squared  correla- 
tion (TSC),  maximizes  the  sum  capacity  in  the  synchronous  CDMA  system  over  AWGN. 
Furthermore,  users  in  the  system  can  be  supported  with  a uniform  target  SIR  using  a mini- 
mum amount  of  total  transmission  power  when  the  equal-energy  WBE  sequences  are  used 
as  spreading  sequences  and  the  MMSE  receivers  are  used  for  signal  reception.  Weights  of 
the  MMSE  receiver  turn  out  to  be  the  same  as  those  of  the  matched-filter  receiver,  which 
is  an  attractive  property  to  simplify  the  receiver  design.  A solution  to  this  problem  in 
the  general  case  of  unequal  sequence  power  constraints  [49]  states  that  a class  of  general- 
ized WBE  (GWBE)  sequences  is  the  set  that  maximizes  the  sum  capacity  of  a multi-user 
synchronous  CDMA  system.  The  optimization  problem  of  power  control  and  spreading 
sequence  allocation  in  conjunction  with  MMSE  signal  reception  so  that  the  SIR  targets 
of  all  users  are  met  with  the  minimum  total  received  power  in  a synchronous  single-cell 
reverse-link  CDMA  system  is  solved  by  Viswanath  et  al.  [50].  Orthogonal  sequences  are 
optimal  when  the  number  of  users  is  less  than  the  spreading  gain.  When  the  number  of 
users  is  larger  than  the  spreading  gain,  the  users  in  the  system  are  classified  as  two  classes, 
called  oversized  users  and  nonoversized  users,  according  to  their  SIR  requirements.  Or- 
thogonal spreading  sequences  are  allocated  to  the  oversized  users,  and  generalized  WBE 
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sequences  are  allocated  to  the  nonoversized  users.  Similar  to  the  WBE  sequences,  MMSE 
receiver  weights  are  the  same  as  those  of  the  matched-filter  receiver  with  this  scheme.  A 
distributed  algorithm  that  minimizes  TSC  iteratively  is  suggested  by  Ulukus  [51]  to  obtain 
the  WBE  sequences  when  a uniform  target  SIR  is  required.  A variant  of  the  sequence 
optimization  problem  in  a multicarrier  setting  has  been  solved  independently  by  Lok  and 
Wong  [52],  Third  generation  (3G)  wireless  communications  will  be  dominated  by  multi- 
media  applications.  High  spectral  efficiency,  high  power  efficiency,  and  flexible  data  rates 
are  the  main  focuses  for  3G  wireless  networks.  A simple  way  to  support  variable  trans- 
mission rates  in  CDMA  systems  is  based  on  the  multicode  approach,  in  which  multiple 
spreading  sequences  for  multiple  transmission  streams  are  assigned  to  a user  when  its  rate 
requirement  exceeds  the  basic  level.  The  system  proposed  by  Wong  [46]  explores  the 
potential  of  sequence  optimization  in  a multicode  CDMA  system,  and  solves  the  optimiza- 
tion problem  of  minimizing  the  transmission  power  for  a subset  of  users  in  the  system, 
which  turns  out  to  be  the  general  solution  of  the  sequence  optimization  problem  with  ad- 
ditive colored  noise.  An  iterative  procedure  based  on  the  power  method  is  also  proposed 
[46]  to  construct  the  spreading  sequences  iteratively.  This  algorithm  has  the  advantage 
that  it  is  suitable  to  all  three  application  scenarios  for  a multicode  CDMA  system,  includ- 
ing single-stream  adaptation,  multi-stream  adaptation  and  all-stream  adaptation.  Sequence 
optimization  and  sum  capacity  in  a CDMA  system  over  an  AWGN  channel  have  been 
thoroughly  understood;  more  research  has  been  focusing  on  more  complex  channel  mod- 
els including  asynchronous  and  multipath  fading  channels.  User  capacity  of  a single-cell 
symbol-asynchronous  CDMA  system  with  the  matched  filter  signal  reception  has  been 
addressed  by  Ulukus  [53,  54],  It  is  claimed  that  the  user  capacity  of  a single-cell  symbol- 
asynchronous  CDMA  system  is  the  same  as  that  of  a single-cell  synchronous  system  when 
the  users’  signature  sequences  are  the  class  of  sequences  that  minimizes  the  TSC.  Trans- 
mitter optimization  in  a multicarrier  setting  over  a fading  channel  has  been  attempted  by 
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Lok  and  Wong  [52].  Centralized  and  decentralized  algorithms  based  on  the  method  of  La- 
grange multiplier  [55]  have  been  devised  to  numerically  solve  the  optimization  problem  of 
minimizing  the  total  transmission  power.  Network  capacity  and  admissibility  of  imperfect 
power-controlled  CDMA  systems  with  linear  receivers  in  a fading  environment  has  been 
studied  [56],  assuming  known  received  power  distributions.  Joint  spreading  sequence  and 
receiver  adaptation  in  a point-to-point  DS-CDMA  system  over  multipath  fading  channel 
has  been  examined  by  Rajappan  and  Honig  [57],  A distributed  adaptation  in  which  each 
user  adapts  to  optimize  its  own  performance  and  a centralized  adaptation  in  which  all  re- 
ceivers are  jointly  adapted  to  minimize  the  sum  MSE  over  all  users  are  described.  Adaptive 
algorithms  are  also  introduced  [58]  to  sequentially  update  spreading  sequences  to  search 
“good  sequences”  for  interference  avoidance  in  CDMA  systems  with  multipath.  The  opti- 
mality criteria  employed  include  effective  TSC,  which  is  related  to  the  receiver  output  SIR, 
sum-MMSE,  and  information-theoretic  sum  capacity. 

Adaptive  procedures  need  to  be  developed  to  facilitate  the  application  of  the  pro- 
posed optimal  sequences  in  practical  communication  systems.  The  adaptive  algorithm 
proposed  by  Wong  and  Lok  [46]  has  the  advantage  that  it  is  suitable  to  many  application 
scenarios  for  multicode  CDMA  systems.  The  drawback  of  this  algorithm  is  that  it  re- 
quires a perfect  estimate  of  the  correlation  matrix  of  the  received  signal  and  assumes  that 
the  linear  receiver  always  works  with  the  optimal  weight  coefficients.  These  assumptions 
limit  the  application  of  the  algorithm  in  practical  systems.  One  objective  of  this  work  is  to 
develop  an  adaptive  CDMA  structure  that  eliminates  these  stringent  requirements,  and  is 
more  practical  to  use  in  real  communication  systems.  Sequence  optimization  over  complex 
channel  models  such  as  general  fading  and  multipath  has  been  attempted  using  numerical 
simulation.  However,  theoretical  solutions  are  not  established  in  the  current  literature.  A 
second  objective  of  this  work  is  to  explore  the  sequence  optimization  problem  in  CDMA 
systems  with  more  complex  channel  conditions. 
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1.3  Data  Link  Layer  Protocols 

All  the  works  discussed  above  focus  primarily  on  physical  layer  issues.  Due  to 
the  existence  of  MAI,  packet  error  control  protocols  have  to  be  applied  in  the  data  link 
layer  to  ensure  error-free  transmission  of  packets.  It  is  known  that  random  access  packet 
switching  protocols,  such  as  slotted  ALOHA,  are  more  efficient  than  fixed  access  schemes 
such  as  FDMA  or  TDMA  for  bursty  traffic.  The  key  performance  measure  of  random 
access  networks  is  their  throughput-delay  characteristics.  However,  the  throughput-delay 
performance  for  random  access  protocols  decreases  dramatically  at  high  traffic  loads  when 
collisions  happen.  In  CDMA  systems,  employment  of  spreading  sequences  allows  simul- 
taneous transmission  of  multiple  users.  An  early  work  of  Raychaudhuri  [59]  suggests  the 
possibility  of  using  CDMA  systems  in  a random  access  mode  of  operation.  It  has  been 
shown  that,  at  high  traffic  loads,  it  is  more  efficient  in  terms  of  network  utilization  to  use 
CDMA  in  conjunction  with  random  access  protocols  than  to  use  only  random  access  pro- 
tocols or  fixed  access  schemes  alone  [59].  It  has  been  widely  realized  that  DS-CDMA 
signaling  is  an  attractive  multiple  access  technique  in  packet  radio  networks  due  to  its 
anti-multipath  and  narrowband  jamming  rejection  capabilities  [60], 

Pure  automatic  repeat  request  (ARQ),  type-I  hybrid  ARQ,  and  type-II  hybrid  ARQ 
are  data  link  layer  protocols,  with  increasing  complexities,  which  try  to  provide  error- 
free  transmission  for  packet  radio  systems  over  noisy  communication  channels.  Because 
of  the  incorporation  of  forward  error  correcting  capability  into  the  feedback  error  con- 
trol protocols,  hybrid  ARQ  schemes  usually  give  higher  average  system  throughputs  and 
lower  average  system  delays  than  the  pure  ARQ  scheme.  Hence,  hybrid  ARQ  schemes 
are  often  employed  to  improve  the  throughput-delay  performance  and  to  increase  network 
reliability.  An  early  work  of  Hanratty  and  Stuber  [61]  evaluates  the  throughput-delay  per- 
formance of  a type-I  hybrid  ARQ  protocol  in  a slotted  DS-CDMA  network  operating  in  a 
hostile  jamming  environment.  It  is  shown  that,  for  a given  population  size,  traffic  intensity, 
and  bit  energy/jammer  noise  ratio,  there  is  an  optimal  probability  of  retransmission,  code 
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rate,  and  processing  gain  that  optimizes  network  performance  in  the  presence  of  worst-case 
pulse  jamming.  Since  then  much  research  [62-65]  has  been  done  to  evaluate  the  delay  and 
throughput  performance  of  CDMA  packet  radio  networks  under  various  protocols  and  er- 
ror correction  schemes.  All  these  works  show  an  improved  throughput-delay  performance 
and  reliability  of  the  CDMA  systems.  One  issue  of  great  concern  in  modem  wireless  infor- 
mation systems  is  to  support  multimedia  services  (voice,  data,  video,  and  others).  These 
services  are  usually  partitioned  into  different  traffic  classes  according  to  the  information 
rate  (bandwidth)  and  quality  of  service  (QoS)  requirements.  Given  a fixed  bandwidth,  QoS 
requirements  can  be  satisfied  by  an  appropriate  assignment  of  transmission  power  and  pro- 
cessing gain  to  packets  in  each  class.  In  the  works  of  Manji  [65],  Shen  [66]  and  Kim  [67], 
the  effect  of  power  assignment  is  analyzed,  and  power  assignment/control  algorithms  are 
proposed  to  obtain  the  best  performance  of  the  integrated  transmission  system. 

In  order  to  fully  utilize  the  advantage  of  the  technique  of  sequence  optimization 
developed  in  the  physical  layer,  we  propose  a data  link  protocol  that  can  integrate  the 
usage  of  optimal  sequences  in  a CDMA-based  wireless  link.  The  objective  is  to  make 
the  most  efficient  use  of  the  transmission  power  to  maximize  the  link  throughput  under  a 
transmission  delay  constraint. 

1.4  Chapter  Organization 

In  Chapter  2,  we  describe  a general  synchronous  DS-CDMA  system  model,  as  well 
as  the  structure  and  property  of  the  optimal  linear  MMSE  receiver.  In  Chapter  3,  we  present 
two  joint  transmitter-receiver  adaptation  schemes  to  utilize  the  optimal  sequences  devel- 
oped earlier  [46,  50,  52]  in  CDMA  systems  over  AWGN  channel.  One  scheme  adapts  in 
a distributed  manner,  and  the  other  one  adapts  in  a centralized  manner,  In  both  schemes, 
training  symbols  are  transmitted  in  blocks  so  that  the  transmitter  spreading  sequences  and 
the  MMSE  receiver  weights  are  adapted  simultaneously.  In  Chapter  4,  we  study  the  power 
and  spreading  sequence  allocation  in  a CDMA  forward  link  over  a slow  flat  fading  channel. 
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We  propose  two  sequence  allocation  schemes,  the  minimum-ETSC  scheme  for  MMSE  sig- 
nal reception  and  the  minimum-TSC  scheme  for  matched  filter  signal  reception,  to  support 
users  in  the  system  with  a uniform  target  SIR.  Adaptive  structure  and  algorithms  are  also 
developed  to  facilitate  the  implementation  of  the  minimum-ETSC  sequences  in  practical 
communication  systems.  In  Chapter  5,  we  investigate  sequence  optimization  in  a CDMA 
point-to-point  transmission  over  multipath  fading  channels,  in  which  all  signals  undergo 
the  same  multipath  channel  response.  A set  of  optimal  spreading  sequences  is  constructed 
to  support  users  in  the  system  with  a uniform  target  SIR,  when  a user  capacity  constraint 
is  satisfied.  In  Chapter  6,  we  introduce  a type-I  hybrid  ARQ  protocol  that  can  integrate 
the  usage  of  the  optimally  selected  spreading  sequences  in  the  physical  layer  to  allow  si- 
multaneous transmission  of  multiple  packets  over  a CDMA-based  wireless  link.  Finally, 
we  conclude  this  work  in  Chapter  7 by  exploring  some  possible  further  research  directions 
based  on  the  results  obtained  in  previous  chapters. 


CHAPTER  2 

SYNCHRONOUS  DS-CDMA  SYSTEM  AND  MMSE  SIGNAL  RECEPTION 


We  describe  the  model  of  a synchronous  DS-CDMA  system  in  this  chapter.  We 
consider  a general  CDMA  forward-link  transmission  over  multipath  fading  channel.  Each 
user’s  signal  is  detected  using  optimal  MMSE  receivers  independently  and  simultaneously. 


2.1  Signal  and  System  Model 

We  assume  that  there  are  K users  in  the  system.  Each  user  generates  one  stream  of 
data  symbols  which  is  given  by 

( b ^ b ^ b ^ ) 

for  1 < k < K.  For  binary  communication,  b^  £ {+1,  — 1}-  We  assume  that  the  data 
symbols  are  independent  random  variables  with  zero  mean  and  unit  variance.  The 
periodic  spreading  sequence  of  the  kth  user,  given  by 

/ ( k ) (/c)  (fc)  \ 


is  used  to  spread  the  data  stream  of  the  kt h user.  For  convenience,  we  use  the  notation  a k to 
denote  the  column  vector  [a[,fe\  a[k\  . . . , which  contains  one  period  of  the  spreading 

sequence  of  the  kth  user.  The  fcth  data  stream  is  spread  with  a^,  and  then  modulated  to  the 
carrier  frequency  cec  to  give  the  transmitted  signal 


Sk{t)  = Re 


(2.1) 


where  Tc  is  the  chip  interval,  and  ip(t)  is  the  chip  waveform.  We  assume  that  Tc  = 1,  ip(t) 
satisfies  the  Nyquist  criterion  for  zero  inter-chip  interference,  and  f ^ \ip(t)\2dt  = 1. 

We  now  describe  the  channel  model.  We  consider  a CDMA  forward  link  transmis- 
sion with  multipath,  and  assume  that  the  multipath  response  from  the  base-station  to  the 
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A;th  receiver  is 

Lk 

hk(t)  = hk,i6{t  - Tkj),  (2.2) 

z=i 

where  hkj  and  rkyi  are  the  channel  coefficient  and  the  transmission  delay  associated  with 
the  Ith  path,  in  respect,  and  Lk  is  the  number  of  resolvable  paths  observed  at  the  kth  re- 
ceiver. The  received  signal  at  the  kth  receiver  can  be  expressed  as  the  convolution  of  the 
transmitted  signal  s(t)  and  the  multipath  channel  response  hk(t)\ 

K Lk 

rk{t)  = EE  hk,ism  (t  T~kj)  T Tik  (i) , (2.3) 

m= 1 Z=1 


where  nk(t)  represents  AWGN  with  power  spectral  density  (PSD)  rjk  at  the  receiver.  For 
convenience  of  analysis,  we  assume  that  all  receivers  observe  the  same  number  of  resolv- 
able paths,  that  is,  Lk  = L,  for  k = 1,  2,  • • • , K,  and  the  paths  are  chip  synchronous,  that 
is,  Tk}i  = (l  - 1 )TC,  for  l = 1, 2,  • ■ • , L.  Hence,  the  complex  baseband  representation  of  the 
received  signal  at  the  kth  receiver  is  given  by 

K L oo 

hm  = E E E -(«+!-  ina + <2.4) 

m= 1 Z=1  i— — 00 

Received  signal  is  passed  through  a chip-matched  filter  and  the  filter  output  is  sam- 
pled at  the  chip  rate.  Due  to  multipath  convolution,  there  are  (N  + L — 1)  samples  for 
every  data  symbol,  which  overlap  with  the  samples  of  the  neighboring  symbols  by  (L  - 1) 
samples.  For  instance,  the  (N  + L - 1)  samples  observed  at  the  output  of  the  chip-matched 
filter  in  the  interval  [0,  (N  + L - 1 )TC)  account  for  the  0th  symbol  of  all  the  data  streams. 
For  simplicity,  we  assume  that  L « N and  there  is  no  significant  ISI.  To  detect  the 
0th  symbol  of  the  users,  we  arrange  the  (N  + L — 1)  samples  observed  in  the  interval 
[0,  (N  + L — 1 )TC)  into  an  (N  + L - l)-dimensional  column  vector  rfc,  which  can  be 
expressed  in  matrix  form  as 

K 

rk  = Hfc  ^2  f)am  + nfc, 

m= 1 


(2.5) 
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where  6o"^am  is  the  contribution  from  the  rath  user  for  1 < m < K,  nk  is  an  AWGN  noise 
vector  with  elementwise  variance  r]k,  and 


1 

?r 

0 

...  o 

0 

0 

• 0 

1 

O 

hk,2 

hk,  i 

...  o 

0 

0 

• 0 

0 

hk,L 

hk,L- i 

• • • hk,  i 

0 

0 

• 0 

0 

0 

hk,L 

• • • hk,  2 

hk,  i 

0 

• 0 

0 

• • O 

0 

...  o 

hk,L 

hk,L- i • 

• hk,\ 

0 

• O 

0 

...  o 

0 

hk,L 

■ hk,  2 

hk,  i 

0 

0 

...  o 

0 

0 

• hk,L 

hk,L- i 

1 

o 

0 

...  o 

0 

0 

■■  0 

hk,L  . 

is  the  multipath  channel  response  matrix. 

2.2  MMSE  Signal  Reception 

We  employ  the  adaptive  MMSE  receiver  structure  proposed  by  Pateros  and  Saulnier 
[27],  Miller  [28]  and  Madhow  [30].  All  data  streams  are  detected  simultaneously  and 
independently.  The  chip-matched  filter  output  rfc  is  fed  into  an  FIR  filter,  which  can  be 
modeled  as  an  (N  + L - l)-tap  fractionally  spaced  tapped  delay  line  (TDL),  and  the  filter 
weights  are  chosen  to  minimize  the  MSE  of  the  FIR  filter  [46]: 

w*  = R^Hfeafc,  (2.7) 

where  RTfc  is  the  (N  + L - 1)  x (Ar  + L - 1)  total  correlation  matrix  observed  at  the 
receiver,  which  is  given  by 

RTk  = E[rfcrf] 

K 

= H*  + T]kI 

m= 1 

= H*ArAfH?  + »fcL  (2.8) 

In  above,  AT  = [ax  a2  • • • a.K]  is  an  N x K matrix  formed  by  grouping  all  K spreading 
sequence  vectors.  The  decision  statistic  Zk  — w^r k is  hard-limited  to  obtain  the  estimate 
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of  the  Oth  symbol  of  the  fcth  user.  A simple  analysis  shows  that  the  optimal  MMSE  receiver 
weights  also  maximize  the  receiver  output  SIR  [46],  which  is  given  by 

SIRfe  = af  Hf  Rfc-:Hfcafc,  (2.9) 

where  Rfc  is  the  noise-plus-interference  correlation  matrix  observed  by  the  kih  user,  which 
is  defined  by 

Rfc  = Rrfc  — Hk^-k^k  H k 

= HfcA*A?Hf + TfcI.  (2.10) 

In  above,  A*  = [ax  • • • afc_x  afc+i  • • • aK\  is  an  N x (I<  - 1)  matrix  formed  by  grouping 
all  K spreading  sequence  vectors  except  a^.  For  ease  of  analysis,  we  derive  another  useful 
expression  of  the  SIR  using  the  matrix  inversion  formula  [68]: 

(2.U) 

Eq.  (2.9)  indicates  that  the  SIR  achieved  by  the  MMSE  receiver  depends  on  the 
choice  of  the  spreading  sequences  as  well  as  the  transmission  powers  of  the  users.  When 
the  system  resources  (characterized  in  terms  of  the  bandwidth  and  transmission  power)  are 
limited,  a good  power  and  sequence  allocation  scheme  is  vital  to  satisfy  the  QoS  require- 


ments of  all  users. 


CHAPTER  3 

JOINT  TRANSMITTER-RECEIVER  ADAPTATION  OVER  AWGN  CHANNEL 

Sequence  optimization  is  a physical  layer  technique  that  can  be  characterized  as 
choosing  a set  of  spreading  sequences  so  that  SIR  targets  are  achieved  for  signals  of  the 
selected  users  with  a minimum  amount  of  total  transmission  power.  In  this  chapter,  we 
propose  a joint  transmitter-receiver  adaptive  structure  to  utilize  the  optimal  sequences  for 
a CDMA  system  over  an  AWGN  channel.  Details  of  the  system  model  and  the  structure  of 
the  MMSE  receiver  are  given  in  Section  3.1.  For  convenience  of  description,  a brief  sum- 
mary of  the  transmitter  sequence  optimization  in  CDMA  systems  over  an  AWGN  channel 
based  on  the  results  of  Wong  and  Lok  [46,  52]  and  Viswanath  et  al.  [50]  is  presented  in 
Section  3.2.  The  sequence  optimization  problem  is  approached  from  three  different  view- 
points, namely,  single-user  viewpoint,  multi-user  viewpoint  and  all-user  viewpoint,  which 
represent  three  different  application  scenarios  in  wireless  communication  systems.  In  Sec- 
tion 3.3,  we  propose  a joint  transmitter-receiver  adaptive  structure.  The  receiver  weights 
are  adapted  based  on  the  recursive  least  square  (RLS)  algorithm.  Information  collected  at 
the  receiver  is  used  to  update  the  spreading  sequences  in  a distributed  manner.  The  updated 
spreading  sequences  are  then  fed  back  to  the  transmitter  and  are  used  for  the  transmission 
of  the  next  block  of  data.  A modified  structure  of  this  joint  adaptation,  which  incorporates 
the  channel  estimation,  is  provided  in  Section  3.4.  Transmitter  spreading  sequences  are  up- 
dated at  the  base-station  in  a centralized  manner  with  this  modified  approach.  Conclusions 
are  drawn  at  Section  3.5. 

3.1  System  Model 

We  describe  the  system  model  in  this  section.  We  consider  a synchronous  CDMA 
system  over  an  AWGN  channel  with  assumptions  that  channel  gains  from  the  base  station 
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to  all  the  users  are  the  same,  and  power  spectral  densities  (PSD)  of  the  thermal  noise  at  all 
the  receivers  are  the  same,  such  that,  Hfc  = 1 and  rjk  = r]  for  k = 1, 2,  • • • , K.  With  these 
assumptions,  the  general  signal  model  described  in  Chapter  2 can  be  simplified.  There  are 
N samples  for  every  data  symbol  at  the  output  of  the  chip-matched  filter,  which  can  be 
arranged  into  an  ^-dimensional  column  vector 

r = 5Z^m)am  + n'  (3-1) 

m— 1 

Accordingly,  the  optimal  MMSE  receiver  weights  are  given  by 

wfc  = IVr1ak,  (3.2) 

where  Rr  is  the  total  correlation  matrix  observed  at  the  receivers,  which  is  given  by 

Rj1  = AtAt  + 77I  - (3.3) 


We  note  that  the  total  correlation  matrices  observed  by  all  the  users  are  the  same  according 
to  the  system  model  we  adopted.  With  the  optimal  MMSE  weights,  the  receiver  output 
SIR  is  given  by 

SIRfc  = afR^a*,  (3.4) 

for  A:  = 1,  2,  • • • , K,  where  Rfc  is  the  noise-plus-interference  correlation  matrix,  which  is 
defined  as 

Rfc  = Rr  - afcaf . (3.5) 


An  equivalent  expression  of  the  SIR  using  the  matrix  inversion  formula  can  be  derived  as 
well: 


SIRfc 
1 + SIRfc 


a^R  Tak- 


(3.6) 


3.2  Transmitter  Optimization  and  Adaptation 


Before  the  description  of  the  joint  transmitter-receiver  adaptation,  a brief  summary 
of  the  transmitter  optimization  in  CDMA  systems  over  an  AWGN  channel  based  on  the 
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results  by  Wong  and  Lok  [46,  52]  and  Viswanath  et  al.  [50]  is  necessary.  Transmitter 
optimization  can  be  addressed  as  the  optimization  problem  of  choosing  an  optimal  set 
of  spreading  sequences  with  the  minimum  total  transmission  power  to  satisfy  a prede- 
termined SIR  requirement.  For  simplicity,  we  assume  a uniform  target  SIR  for  all  the 
data  streams.  Transmitter  optimization  can  be  approached  from  three  different  viewpoints, 
namely,  single-user  optimization,  multi-user  optimization  and  all-user  optimization.  In  the 
single-user  optimization  case,  only  one  of  the  users  is  allowed  to  optimize  its  transmission 
sequence  and  the  transmission  sequences  of  all  other  users  are  fixed.  In  the  multi-user  op- 
timization case,  more  than  one  users  are  allowed  to  optimize  their  transmission  sequences. 
In  the  all-user  optimization  case,  all  users  in  the  system  are  allowed  to  optimize  their  trans- 
mission sequences.  These  three  viewpoints  represent  three  different  application  scenarios 
in  wireless  communication  systems. 

3.2.1  Single-user  Optimization 

A possible  practical  scenario  described  by  this  problem  is  that  a new  user  with  high 
priority  is  admitted  into  the  system  and  is  about  to  optimize  its  transmission.  Without  loss 
of  generality,  we  assume  spreading  sequence  optimization  of  the  first  user.  The  objective  of 
transmitter  optimization  is  to  achieve  a target  SIR  for  the  first  user  with  a minimum  amount 
of  transmission  power.  Employing  the  optimal  MMSE  receiver  discussed  in  Section  2.2, 
the  single-user  optimization  problem  can  be  expressed  mathematically  as 

min||a1||2 

ai 

subject  to 

af  R^ai  = 7.  (3.7) 

A closed  form  solution  for  this  optimization  problem  can  be  readily  obtained  by 
the  method  of  Lagrange  multiplier  [55],  The  spreading  sequence  ai  should  be  chosen  as 
the  eigenvector  associated  with  the  smallest  eigenvalue,  A m;n,  of  Ri . With  this  choice,  the 
minimum  transmission  power  is  given  by  ||ai  ||  — Amm7- 
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There  is  a simple  intuitive  inteipretation  for  this  optimal  choice  of  the  spreading 
sequence.  Each  eigenvector  of  Ri  represents  a “channel”  in  the  CDMA  system  and  the 
corresponding  eigenvalue  indicates  the  amount  of  interference  in  that  channel.  A large 
eigenvalue  indicates  a channel  with  a large  amount  of  interference.  To  minimize  the  trans- 
mission power,  a new  spreading  sequence  should,  of  course,  be  allocated  into  the  channel 
that  has  the  least  amount  of  interference.  When  the  number  of  the  existing  transmission 
sequences  is  less  than  the  processing  gain  N , the  optimal  spreading  sequence  is  chosen  as 
orthogonal  to  the  existing  sequences.  The  admission  of  the  new  sequence  does  not  degrade 
the  performance  of  the  original  data  signals  due  to  orthogonality.  When  the  number  of 
the  existing  transmission  sequences  is  equal  to  or  larger  than  the  processing  gain  N,  the 
additional  code  sequence  can  still  be  squeezed  into  the  CDMA  system.  Addition  of  the  ad- 
ditional sequence,  of  course,  degrades  the  performance  of  the  existing  transmissions.  The 
optimal  choice  of  the  spreading  sequence  makes  this  degradation  the  least  severe  [46], 
3.2.2  Multi-user  Optimization 

Multi-user  optimization  describes  the  situation  when  two  or  more  new  users  are 
admitted  into  the  system  and  are  about  to  optimize  their  transmissions.  Without  loss  of 
generality,  we  consider  transmitter  optimization  of  the  first  M users.  The  goal  of  the 
transmitter  optimization  is  to  achieve  a predetermined  uniform  target  SIR  7 for  all  the  M 
new  data  streams  with  a minimum  amount  of  transmission  power.  Mathematically,  the 
optimization  problem  can  be  expressed  in  the  following  form: 


M 


mm 

ai,a2,...,a  M 


m=  1 


subject  to 


= 7 


(3.8) 


for  1 < k < M. 


Compared  with  the  single-user  optimization  problem,  the  multi-user  optimization 
problem  in  (3.8)  is  much  more  complex.  Due  to  its  complexity,  the  matrix  optimization 
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problem  in  (3.8)  is  simplified  into  a scalar  optimization  problem  that  is  much  easier  to 
solve.  The  solution  of  the  scalar  optimization  problem  is  then  used  to  construct  the  optimal 
spreading  sequences. 

Details  of  the  discussion  of  the  optimal  solution  for  the  case  of  M < N can  be 
found  in  Wong  and  Lok  [46].  Here,  we  just  restate  the  optimal  solution  briefly.  Let  Ai..m  = 
[ai  • • • aM]  be  an  N x M matrix  by  grouping  all  the  M new  spreading  sequence  vectors 
together,  and  let  Ri..m  be  the  noise-plus-interference  matrix  observed  by  the  M new  users, 
which  can  be  expressed  as  Ri..m  = Rr  — We  ta^e  t^ie  sPectral  factorization 

of  Ri..m  as  Ri..m  = UAUH  , and  let  the  eigenvalues  An,  for  n = 1, 2,  • • • , N,  be  arranged 
in  a descending  order  in  the  diagonal  matrix  A.  The  original  matrix  optimization  problem 
can  be  simplified  as  a scalar  optimization  problem  as 

M 

min  V \N-M+k/3k 

01,02, -,0m  f-f 

k= 1 


subject  to 


E 


ft 

1 + (3k 


0 < A < /?2  < 


Me(  7) 
< 


(3.9) 


where  e (7)  = 

Let  fa,  for  k = 1, . . . , M,  be  the  solution  of  the  scalar  optimization  problem  in 
(3.9).  Construct  aniVxM  matrix  U from  the  last  M columns  (keeping  the  original  order) 
of  U and  construct  the  M x M unitary  matrix  W such  that 


Then, 


diag 


1 

+ Arij 

e(7) 

w 

WH 

= 

0M 

1 ~P~0M  - 

e(7)  . 

yA 


N—M+l 


Pi 
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H 


Ai  ..m  — U 


v/A  nPm 


(3.10) 
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is  a solution  to  the  optimization  problem  in  (3.8)  and  the  resulting  minimum  transmis- 
sion power  is  Y^k=  1 A N-M+kPk • The  intuition  obtained  from  the  single-user  optimization 
case  extends  to  the  multi-user  optimization  problem  above.  The  optimal  sequences  for 
the  multi-user  optimization  problem  lie  in  the  subspace  spanned  by  the  eigenvectors  cor- 
responding to  the  M smallest  eigenvalues  of  Rl.mv  which  corresponds  to  the  M CDMA 
channels  that  have  less  amount  of  interference. 

The  multi-user  optimization  for  the  case  of  M < N given  by  Wong  and  Lok  [46] 
can  be  readily  extended  to  the  case  of  M > N.  All  the  M new  users  can  be  supported  with 
a uniform  target  SIR  7 if  and  only  if  a user  capacity  constraint 

Me(7)<Ar  (3.11) 


is  satisfied.  The  proof  of  this  user  capacity  constraint  can  be  easily  worked  out  following  an 
approach  similar  to  the  one  employed  by  Wong  and  Lok  [52],  Following  the  definition  and 
approach  in  Wong  and  Lok  [46],  we  can  simplify  the  original  matrix  optimization  problem 
as  a scalar  optimization  problem  as 


N 


min 

0\,01y,0N 


a kPh 


k= 1 


subject  to 


E 


ft 

1 + Pk 


0 < ft  < & < 


Me(  7) 
< Pn- 


(3.12) 


Let  pk,  for  k = 1, . . . , N,  be  the  solution  of  the  scalar  optimization  problem  in  (3.12),  and 
construct  an  M x M unitary  matrix  W such  that 


~ 

r /3>  1 

■ 

1+01 

0 N 

w 

WH 

0 

0 _ 

e(7) 


e(7) 


diag 
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Then, 


-A-i  ,,m  — U 


V^iPi 


0 


WH 


(3.13) 


V^nPn 

is  a solution  to  the  optimization  problem  in  (3.8)  and  the  resulting  minimum  transmis- 
sion power  is  %2kLi  ^k(3k-  We  can  observe  from  this  sequence  allocation  that  the  M new 
spreading  sequences  are  squeezed  into  the  entire  space  occupied  by  all  the  N eigenvectors 

of  Ri..m- 

3.2.3  All-user  Optimization 


We  now  consider  the  problem  of  finding  the  optimal  spreading  sequences  from  the 
point  of  view  of  the  whole  system.  All  users  are  allowed  to  optimize  their  transmitters. 
A possible  scenario  described  by  this  problem  is  that  after  a user  leaves  the  system,  all 
the  remaining  users  need  to  adjust  their  spreading  sequences  in  order  to  reduce  the  overall 
transmission  power.  In  this  case,  the  goal  of  optimization  is  to  achieve  the  target  SIRs  for 
all  the  data  signals  in  the  system  with  the  minimum  amount  of  total  transmission  power. 
Mathematically,  the  optimization  problem  can  be  expressed  as 

K 
1=  1 


subject  to 

afR^a*;  = 7 (3.14) 


for  1 < k < K. 

A general  solution  to  this  optimization  problem  is  available  in  Viswanath  et  al. 
[50].  The  solution  in  a multicarrier  setting  is  also  available  in  lok  and  Wong  [52],  For  the 
case  of  K < N,  the  constraint  in  (3.14)  can  always  be  satisfied  and  the  optimal  choice  of 
sequences  is  characterized  by 

A£  Ar  = 7771.  (3.15) 

This  means  that  the  optimal  choice  of  the  spreading  sequences  is  orthogonal  sequences. 
The  minimum  total  transmission  power  needed  is  AT77.  For  the  case  of  K > N,  which 
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is  the  case  that  we  are  more  interested  in,  it  is  argued  that  K users  are  admissible  into  the 
system  with  the  SIR  constraint  7 if  and  only  if  the  SIR  target  7,  the  processing  gain  N,  and 
the  total  number  of  users  K satisfy  the  constraint 


Ke( 7)  < N. 


(3.16) 


Here,  we  follow  the  notation  in  Tse  and  Hanly  [35]  to  denote  e(7)  = ^77-  The  quantity 
e(7)  represents  the  effective  bandwidth  [35]  that  a user  takes  when  achieving  the  MMSE 
receiver  output  SIR  target  7.  There  is  a simple  interpretation  for  this  constraint.  All  K 
users  are  admissible  if  and  only  if  the  sum  of  their  effective  bandwidth  is  less  than  the 
spreading  gain  of  the  system.  If  the  condition  in  (3.16)  is  satisfied,  the  optimal  choices  of 
the  spreading  sequences  are  characterized  by 

rye  (7) 


A^Ah  = 

T N/K  — e('j) 


I. 


(3.17) 


The  minimum  total  transmission  power  needed  is  then 


K 


Nr/e (7) 

: N/K  - e(7) ' 

3.2.4  An  Iterative  Algorithm  for  Transmitter  Adaptation 


= 51  IN 


(3.18) 


k- 1 


From  the  transmitter  optimization  solutions  stated  above,  the  base-station  trans- 
mitter needs  to  estimate  the  channel  parameters  for  the  users  in  the  system,  estimate  the 
correlation  matrix,  and  solve  the  complex  eigen-problem  before  it  can  construct  the  opti- 
mal sequences.  In  a practical  communication  system,  some  users  may  leave  the  system 
and  new  users  may  enter  the  system  at  anytime.  It  is  difficult  for  the  base-station  to  esti- 
mate the  channel  parameters  and  construct  a new  set  of  spreading  sequences  every  time  the 
system  state  changes.  To  solve  the  problem  practically,  an  adaptive  algorithm  based  on  the 
power  method  was  proposed  in  Wong  and  Lok  [46]  to  construct  the  spreading  sequences 
iteratively.  Here,  we  give  a short  illustration  of  this  simple  and  effective  algorithm: 


For  the  A:th  user,  we  update  its  spreading  sequence  vector  in  the  jth  iteration  by 

a k[j  + 1]  = gki^^ijWij],  (3-19) 
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where  the  constant  gain  gk\j]  is  chosen  so  that 

af  [j  + 1]R'71b']afc[i  + 1]  = 7 (1  - af  [j]^[jW[j])  . (3.20) 

Discussion  of  the  convergence  property  of  this  algorithm  can  be  found  in  Wong  and 
Lok  [46].  One  of  the  most  attractive  feature  of  this  algorithm  is  that  it  can  be  employed 
for  sequence  adaptation  in  all  the  three  cases  described  above.  The  drawback  of  this  trans- 
mitter adaptation  algorithm  is  that  it  requires  explicit  knowledge  of  the  correlation  matrix 
Rr[j],  and  assumes  that  the  receiver  always  works  with  the  optimal  weights.  Both  as- 
sumptions are  impractical  in  real  communication  systems.  In  order  to  solve  the  problem 
practically,  we  are  going  to  develop  an  iterative  procedure,  which  adapts  the  transmitter 
spreading  sequence  and  the  receiver  weights  simultaneously  when  the  channel  conditions 
and  interference  parameters  are  unknown. 

3.3  RLS-based  Joint  Transmitter-Receiver  Adaptation 

In  practice,  there  should  be  a central  controller  that  is  responsible  of  admitting  users 
into  the  system  and  releasing  users  from  the  system.  The  central  controller  also  determines 
whether  single-user  optimization,  multi-user  optimization,  or  all-user  optimization  should 
be  performed.  Once  this  decision  is  made  and  the  users  are  informed,  all  a user  needs 
to  do  is  to  update  (if  it  is  allowed  to)  its  transmission  sequence  according  to  an  adaptive 
algorithm  to  achieve  the  target  SIR. 

3.3.1  Joint  Transmitter-Receiver  Adaptive  Structure 

The  general  structure  of  the  adaptive  CDMA  system  is  illustrated  in  Fig.  3.1.  Data 
symbols  of  each  user  are  transmitted  in  blocks  using  the  current  spreading  sequence.  We 
assume  that  the  data  signal  of  each  user  in  the  system  is  demodulated  separately  by  an  adap- 
tive MMSE  receiver.  During  the  initial  training  period,  FIR  filter  weights  are  adapted  every 
symbol  according  to  the  MMSE  criterion.  The  spreading  sequence  is  updated  iteratively 
every  block  based  on  the  information  collected  at  the  receiver  during  the  transmission  of 
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Figure  3.1:  Adaptive  CDMA  system 

the  current  block  of  symbols.  The  updated  spreading  sequence  is  then  fed  back  to  the  trans- 
mitter for  the  transmission  of  the  next  block  of  data  symbols.  The  receiver  keeps  adapting 
its  filter  weights,  updating  the  spreading  sequence  and  feeding  back  the  updated  spreading 
sequences.  This  adaptation  process  continues  until  the  receiver  output  SIRs  converge  to 
the  desired  level. 

3.3.2  RLS-Based  Joint  Transmitter-Receiver  Adaptive  Algorithm 

An  RLS  algorithm  [69]  is  used  to  update  the  MMSE  receiver  weights.  The  esti- 
mate of  the  total  correlation  matrix  is  updated  during  the  iteration  of  the  RLS  algorithm. 
A modified  version  of  the  algorithm  in  (3.20)  makes  use  of  this  estimated  correlation  ma- 
trix as  well  as  the  MSE  information  collected  at  the  output  of  the  receiver  to  update  the 
transmitter  spreading  sequences.  The  iterative  procedure  of  the  £;th  user’s  signal  consists 
of  following  steps. 

Step  1.  The  transmitter  transmits  the  j th  block  of  training  symbols  (contains  L 

symbols)  of  the  A:th  user  with  the  current  spreading  sequence  a*  [j] . (The  initial  spreading 
sequence  afc[0]  can  be  set  arbitrarily.) 

Step  2.  The  receiver  employs  the  RLS  algorithm  to  adapt  the  filter  weights  using 
the  jth  block  of  training  symbols.  Receiver  weights  are  adjusted  according  to  following 
equations, 

£k,j,l+ 1 = h,j,l+x  - wk,j,lzk,j,l+ 1, 
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3,1+ 1 
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RT,j,f+ 1—  ^-i+lZA:J,J+l^T,i,i]’ 


A 

wfc,j,m  = wfc,j,(  + kj,  m6'*,  j,i+i> 

to  detect  the  Ith  symbol,  for  l = 1,  2,  • • • , L,  during  the  reception  of  the  jth  block  of  the 
Ath  user,  for  the  A;  = 1, 2,  • • ■ , K. 

Step  3.  At  the  end  of  reception  of  the  jth  block  of  data,  the  MSE  for  the  jth 
block  of  transmission  of  the  Ath  user  is  estimated  as 

L 

L 


1 L 

Ek,j  — Y ^k,j,l  ~ Wk,j,LZk,j,l 


(3.21) 
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The  spreading  sequence  is  updated  according  to 


a*[j  + 1]  9k,j^k,j,L+ 1) 


(3.22) 


where  the  coefficient  is  chosen  so  that 

a k[j  + l]R^,L+iafc[j  + 1]  = 7 kikd-  (3-23) 

The  new  spreading  sequence  a k[j  + 1]  is  fed  back  to  the  transmitter  for  the  transmission 
of  the  (j  + 1 ) th  block  of  the  training  symbols.  The  adaptation  process  repeats  from  Step  1 
until  the  receiver  output  SIR  and  the  transmission  power  converge. 

There  are  several  points  we  would  like  to  address  here.  First,  it  is  noted  that  this 
joint  adaptation  procedure  does  not  require  perfect  knowledge  or  estimation  of  the  total 
correlation  matrix,  and  assumes  no  prior  knowledge  of  the  spreading  sequences  of  other 
users  or  the  noise  floor.  Second,  the  fast  transient  convergence  behavior  of  the  RLS  algo- 
rithm [69]  allows  an  effective  use  of  limited  training  data.  The  RLS  algorithm  can  converge 
very  fast  even  when  the  eigenvalue  spread  of  the  correlation  matrix  is  large  [69],  Further- 
more, in  order  to  further  improve  the  ability  to  track  the  inverse  correlation  matrix  R7A ; in 
the  adaptive  implementation,  we  use  an  exponentially  weighted  RLS  algorithm  [69]  with  a 
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forgetting  factor  A.  General  choices  of  A range  from  0.8  to  1.  Third,  the  main  computation 
complexity  burden  of  the  RLS  algorithm  lies  in  the  update  of  the  inverse  correlation  ma- 
trix. The  computation  complexity  of  this  joint  transmitter-receiver  adaptation  is  similar  to  a 
general  MMSE  receiver  since  the  inverse  correlation  matrix  obtained  in  the  RLS  algorithm 
is  employed  to  update  the  spreading  sequences.  A remaining  question  is  how  frequently 
the  spreading  sequences  should  be  updated,  that  is,  what  should  be  a good  choice  of  the 
symbol  block  length  L.  A discussion  on  the  choice  of  the  block  length  is  given  through 
computer  simulations.  Finally,  we  would  like  to  mention  that  if  a continuous  update  of 
the  spreading  sequences  is  desired  after  the  training  period,  symbol  decisions  made  by  the 
receiver  can  be  employed  to  replace  the  training  symbols. 

3.3.3  Simulation  Analysis 


We  study  the  performance  of  the  proposed  adaptive  implementation  using  computer 
simulations.  The  performance  of  the  system  is  measured  via  the  receiver  output  SIR  that 
is  defined  as 

OTr>  instanteous  power  of  signal  ^ 0/n 

0 1 il  — ~ j ■ f,  j '/ 


average  power  of  inter f erence  and  noise 

With  this  definition,  the  receiver  output  SIR  of  the  llh  symbol  during  the  j th  block  of 

transmission  of  the  kth  signal  is  given  by 


SIRk,j,i  - 
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(3.25) 


where  R k[j]  = a./ [j ] a./^ [j ] + f/I-  The  total  transmission  power  during  the  jth  block 

of  transmission  is  given  by 

K 

Pli\  = £ IKMf.  (3.26) 


k= 1 

In  the  following  simulation  examples,  we  assume  a spreading  factor  of  N = 31.  Originally 
there  are  K — 37  active  users  in  the  system.  The  PSD  rj  of  AWGN  satisfies  ^~°°  - — = 

5. 


In  the  first  simulation  example,  we  study  the  effect  of  the  block  length  L on  the 
convergence  property  of  the  proposed  sequence  update  algorithm.  A training  sequence  of 
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2000  symbols  is  sent  in  blocks  with  block  length  L = 4,  20,  and  100  symbols,  respectively. 
We  consider  the  case  of  single-user  optimization.  The  designated  SIR  is  7 = 4dB.  Fig.  3.2 
shows  the  receiver  output  SIR  and  the  transmission  power  of  the  adapting  user  during  the 
adaptation  process.  With  small  block  lengths,  the  SIR  reaches  the  designated  SIR  very 
quickly  but  with  large  overshoots.  As  the  block  length  increases,  the  overshoots  become 
less  severe.  We  also  observe  that  the  convergence  rate  in  terms  of  the  number  of  training 
symbols  needed  does  not  show  much  difference  for  different  choices  of  block  length.  This 
observation  is  reasonable.  The  transmitter  spreading  sequence  is  updated  more  frequently 
with  a small  block  length.  On  the  other  hand,  each  update  with  large  block  length  is  more 
effective  since  it  employs  a better  estimate  of  the  correlation  matrix.  The  receiver  output 
SIR  can  reach  the  desired  target  very  quickly  through  more  frequent  updates  with  a smaller 
block  length.  However,  frequent  update  means  a waste  of  bandwidth  in  the  feedback  chan- 
nel. Overshoots  of  transmission  power  also  waste  power.  With  a large  block  length,  we  can 
reduce  the  waste  of  feedback  channel  bandwidth  and  transmission  power,  but  the  output 
SIR  reaches  the  target  more  slowly.  A second  observation  is  that  the  transmission  power 
and  the  receiver  output  SIR  converge  at  similar  rates.  This  is  explicit  since  we  are  adapting 
the  receiver  weights  and  transmitter  spreading  sequence  simultaneously  and  the  adaptation 
of  each  of  them  affects  the  adaptation  of  the  other  directly.  Finally,  we  notice  that  the  re- 
ceiver output  SIR  approaches  the  target  SIR,  but  does  not  reach  the  target  SIR  7 = 4dB. 
This  is  because  the  proposed  adaptive  algorithm  can  reach  the  target  SIR  only  when  re- 
ceiver works  with  the  optimal  weights.  With  RLS  adaptation,  due  to  the  inaccuracy  of 
the  estimate  of  inverse  correlation  matrix,  the  weights  are  not  the  optimal  weights,  which 
cause  a deduction  of  SIR  of  about  ldB. 

The  second  simulation  example  examines  the  convergence  property  of  the  case  of 
multiple-user  optimization.  We  assume  that  there  are  M = 3 new  users  admitted  into  the 
system.  We  set  a uniform  target  SIR  of  7 = 4dB  to  all  of  the  3 new  data  signals.  From 
the  discussion  of  the  effect  of  block  length  in  previous  simulation  example,  we  select  the 
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Figure  3.2:  Effect  of  transmission  block  length  in  single-user  optimization 
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block  length  as  L — 50  symbols.  Fig.  3.3  shows  the  adaptation  process.  We  can  see  that  the 
receiver  output  SIRs  of  all  the  3 users  converge  close  to  the  target  SIR  and  the  transmission 
power  approaches  the  minimum  transmission  power. 

The  third  simulation  example  examines  the  convergence  property  of  the  case  of  all- 
user optimization.  We  set  a uniform  target  SIR  of  7 = 4dB  to  all  K = 37  users.  Fig.  3.4 
shows  an  example  of  the  adaptation  process.  We  can  see  that  the  transmission  power  con- 
verges during  the  adaptation  process  and  it  approaches  the  minimum  transmission  power. 
The  receiver  output  SIRs  of  all  the  37  users  also  converge  close  to  the  target  SIR,  but  having 
oscillations. 

In  the  simulation  example  above,  we  assume  that  the  transmitter  can  make  use  of 
the  new  spreading  sequence  immediately  after  it  is  updated.  This  assumption  is  not  valid 
in  practical  systems  since  there  is  always  a transmission  delay  for  the  updated  spreading 
sequence  to  be  fed  back  to  the  transmitter  due  to  the  propagation  and  protocol  delay  and  the 
processing  time.  To  accommodate  this  feedback  delay,  the  joint  adaptation  process  needs 
to  be  modified  slightly.  We  assume  that  there  is  a feedback  delay  of  5 blocks  of  symbols. 
At  the  very  beginning,  the  transmitter  sends  the  first  block  of  symbols  with  the  spreading 
sequence  a*[l],  and  the  receiver  is  synchronized  to  process  the  signal.  The  receiver  updates 
and  feeds  back  the  spreading  sequence  a*; [2]  to  the  transmitter  when  it  finishes  processing 
the  first  block  of  symbols.  The  transmitter  continues  to  use  a*[l]  to  send  5 + 1 blocks 
of  symbols  until  a* [2]  arrives.  By  then,  the  receiver  has  finished  updating  the  spreading 
sequence  to  afc[5  + 2j.  After  this  initial  transient  period,  the  transmitter  receives  a spreading 
sequence  update  a k[j]  every  block  time  and  uses  a k[j\  to  spread  the  (5  + j) th  block  of 
symbols.  The  receiver  receives  the  (5  + j)th  block  of  signal  and  updates  the  spreading 
sequence  to  a^  [5  + j + 1],  In  this  last  simulation  example,  we  assume  a delay  of  5 = 2,  that 
is,  it  takes  100  symbol  transmission  time  for  the  feedback  spreading  sequence  to  arrive  at 
the  transmitter.  The  simulation  result  is  shown  in  Fig.  3.5.  We  can  see  that  the  existence 


Figure  3.3:  Performance  of  multi-user  optimization 
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Figure  3.4:  Performance  of  all-user  optimization 
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of  transmission  delay  does  not  have  serious  impact  on  the  performance  of  the  proposed 
adaptive  structure. 

3.4  Joint  Transmitter-Receiver  Adaptation  with  Channel  Estimation 

The  approach  in  Section  3.4  updates  spreading  sequences  at  the  receivers  in  a dis- 
tributed manner  and  sends  the  newly  updated  sequences  back  to  the  base-station.  This 
sequence  feedback  approach  causes  a waste  of  bandwidth  in  the  feedback  channel,  espe- 
cially when  the  spreading  gain  is  large.  An  alternative  approach  to  be  proposed  in  this 
section  aims  to  reduce  the  bandwidth  requirement  in  the  feedback  channel  by  updating 
spreading  sequences  at  the  base-station  in  a centralized  manner.  A minimum  amount  of 
feedback  information  is  required  to  facilitate  sequence  update.  This  approach  makes  the 
adaptive  scheme  more  practical  and  feasible  to  use  in  real-life  communication  systems. 
We  consider  the  problem  from  the  point  of  view  of  the  whole  system,  and  assume  that  all 
users  are  allowed  to  adapt  their  sequences.  Each  user  is  assumed  to  have  no  information 
from  other  users.  Again,  training  symbols  are  required  and  receiver  weights  of  each  user 
are  adapted  independently  and  simultaneously  based  on  the  MMSE  criterion.  An  estimate 
of  the  received  signal  power  at  each  receiver  is  fed  back  to  the  base-station  and  is  used  for 
the  estimation  of  power  spectral  density  (PSD)  of  the  AWGN,  which  is  then  employed  to 
update  the  spreading  sequences. 

3.4.1  Adaptive  Algorithm 

In  this  approach,  the  base-station  collects  the  feedback  information  from  the  re- 
ceivers and  adapts  the  spreading  sequences  of  all  users  in  a centralized  manner.  The  itera- 
tive procedure  of  the  klh  user  consists  of  following  steps. 

Step  1.  The  base-station  transmits  the  j th  block  of  training  symbols  (contain  L 
data  bits)  of  the  £:th  user  with  the  current  spreading  sequence  a*  [7],  for  1 < k < K.  (The 
initial  spreading  sequence  AT[0]  = [ai[0],  a2[0],  • • • , aA'[0]]  can  be  assigned  arbitrarily.) 
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Figure  3.5:  Multi-user  optimization  with  transmission  delay 


37 


Step  2.  The  weight  vector  of  the  receiver  of  the  kth  user  is  adapted  using 
the  exponentially  weighted  RLS  algorithm  described  in  the  previous  section,  for  1 < l < 
L during  the  jth  block  of  transmission.  The  received  signal  power  for  the  kth  user  is 
estimated  as 


Li 


£ 

;=i 


\zk,j,l\ 


(3.27) 


where  zk,j,i  is  the  chip  matched  filter  output  vector  at  the  kth  receiver  for  the  Ith  symbol 
during  the  jth  block  of  transmission,  This  estimates  is  sent  back  to  the  base-station  through 
a feedback  channel. 

Step  3.  According  to  our  design  approach,  the  base-station  knows  the  current 
spreading  sequence  a k[j)  for  1 < k < I\,  as  well  as  the  current  total  transmission  power 

<3-28> 

fc=i 


Then  estimates  of  the  PSD  of  the  AWGN  for  each  user  can  be  obtained  as 


f)k  = 


N 


(3.29) 


The  PSDs  of  all  users  are  supposed  to  be  the  same  under  the  AWGN  channel.  Hence,  the 
estimate  of  PSD  of  the  AWGN  can  be  taken  as  an  average  of  % over  all  users: 


fj  _ ^fc=l  Vk 

1 K ' 


(3.30) 


Step  4.  Using  the  estimate  fj,  the  base-station  updates  the  spreading  sequences 
of  all  users  sequentially  starting  from  the  first  user.  For  the  kth  user,  its  spreading  sequence 
is  updated  according  to 

a k[i  + 1]  = 9kR^[j]sik\j  + 1],  (3.3  i) 


where  the  coefficient  gk  is  chosen  so  that 


a S\j  + l]R^[#K(j  + 1]  = 7 (1  - ajE'bWbla/ibl)  ■ (3.32) 

In  above,  RT,  I)]  = Ar,b]AT,  ()]+’)  I.andAT,b]  = [aib+1],  - ■ ’ . at-bj+l],  afc[b],  • • • ,ajtb]]- 
At  the  end  of  the  jth  iteration,  the  updated  spreading  sequences  A T[j  + 1]  = [ai[j  + 


38 


1],  a.2 [j  + 1],  • • • , a*:[j  + 1]]  are  employed  to  spread  the  next  block  of  data.  The  adaptation 
process  repeats  from  Steps  1 to  4 until  the  training  period  finishes. 

3.4.2  Simulation  Analysis 

In  this  section,  we  study  the  performance  of  the  proposed  adaptive  implementation 
through  computer  simulations.  We  employ  the  same  sample  system  as  in  Section  3.4. 
According  to  the  discussion  in  Section  3.4,  a block  length  L = 50  symbols  is  chosen  in 
the  following  simulation  examples.  Starting  from  a set  of  independently  chosen  random 
sequences,  all  users  adapt  their  transmissions.  Fig.  3.6  shows  the  receiver  output  SIRs  and 
the  total  transmission  power  during  the  adaptation  process.  It  is  shown  that  the  receiver 
output  SIRs  of  all  users  converge  close  to  the  SIR  target  and  the  total  transmission  power 
approaches  Pmin  as  calculated  in  (3.18)  after  a number  of  iterations.  We  also  observe  that 
due  to  imperfect  channel  parameter  estimation,  there  are  oscillations  for  both  the  SIR  and 
the  transmission  power  curves. 

Compare  Figs.  3.6  to  3.4,  we  notice  that  the  algorithm  proposed  in  this  section  out- 
performs the  algorithm  proposed  in  Section  3.4,  in  the  sense  that  it  converges  faster  and 
exhibits  less  oscillations.  There  are  several  factors  that  cause  this  performance  improve- 
ment. First,  spreading  sequences  are  updated  in  a centralized  manner  in  this  algorithm. 
However,  in  the  previous  algorithm,  all  the  spreading  sequences  are  updated  parallelly  in 
a distributed  manner.  Second,  both  the  receiver  weights  and  MSE  estimate  of  the  algo- 
rithm in  Section  3.4  are  imperfect  estimates,  which  may  affect  the  convergence  of  algo- 
rithm severely.  In  this  algorithm,  the  base-station  has  the  explicit  knowledge  of  the  current 
spreading  sequences,  and  only  the  estimate  of  the  PSD  is  noisy.  Inaccuracy  in  the  es- 
timation of  AWGN  power  spectral  density  does  not  affect  the  eigenvalue  distribution  of 
correlation  matrix  severely  since  its  contribution  to  the  total  correlation  matrix  is  an  iden- 
tity matrix. 
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Figure  3.6:  Performance  of  joint  adaptation  with  channel  estimation 
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3.5  Summary 

Two  joint  transmitter-receiver  adaptation  structures  in  a DS-CDMA  system  over 
the  AWGN  channel  are  described,  and  their  convergence  properties  are  examined  through 
computer  simulations.  In  both  schemes  proposed,  the  MMSE  receiver  weights  and  the 
spreading  sequences  of  all  users  are  updated  simultaneously  to  meet  a uniform  SIR  target. 
In  the  first  scheme,  the  spreading  sequences  are  updated  at  the  receivers  in  a distributed 
manner,  and  are  sent  back  to  the  base-station  for  spreading  the  next  block  of  transmission. 
An  important  advantage  of  this  algorithm  is  its  simplicity.  Compared  to  an  adaptive  CDMA 
system  that  only  provides  MMSE  receiver  adaptation,  this  approach  is  only  slightly  more 
complex.  In  the  second  scheme,  the  spreading  sequences  are  updated  at  the  base-station 
in  a centralized  manner.  The  advantage  of  this  scheme  over  the  first  one  is  that  only  a 
minimum  amount  of  information  is  sent  back  to  the  base-station  from  receivers  to  facilitate 
the  adaptation.  This  approach  reduces  the  bandwidth  requirement  in  the  feedback  channel, 
and  makes  the  adaptive  structure  more  feasible  to  use  in  a practical  communication  system. 
Computer  simulation  results  indicate  that  the  second  scheme  gives  better  performance  than 
the  first  one  in  the  sense  that  it  converges  faster  and  exhibits  less  oscillations  in  the  presence 


of  channel  estimation  error. 


CHAPTER  4 

SEQUENCE  ALLOCATION  IN  A CDMA  FORWARD  LINK 


Most  of  the  previous  work  regarding  sequence  optimization  concerns  either  the 
AWGN  channel  or  the  reverse  link.  Results  about  joint  power  control  and  spreading  se- 
quence allocation  in  the  CDMA  forward  link  are  lacking.  In  this  chapter,  we  consider 
power  control  and  spreading  sequence  allocation  to  support  users  in  a multi-user  single- 
cell forward-link  CDMA  system.  We  assume  that  the  signals  of  different  users  undergo 
independent  slow  flat  fading.  Details  of  the  signal  and  system  model  is  defined  in  Sec- 
tion 4.1.  Two  power  and  sequence  allocation  schemes,  a minimum-ETSC  scheme  with 
MMSE  signal  reception  and  a minimum-TSC  scheme  with  matched-filter  signal  reception, 
are  presented  in  Section  4.2  and  Section  4.3,  respectively.  With  the  proposed  schemes, 
all  users  in  the  system  can  be  supported  with  a uniform  target  SIR.  In  Section  4.4,  perfor- 
mance of  the  proposed  schemes  is  evaluated  and  compared  in  terms  of  power  efficiency. 
We  seek  the  minimum  total  transmission  power  required  to  support  all  users  numerically 
using  a Lagrangian  multiplier  search  method.  In  Section  4.5,  two  adaptive  algorithms  are 
introduced  to  construct  the  spreading  sequences  iteratively.  Convergence  properties  of  the 
adaptive  algorithms  are  examined  using  computer  simulations.  In  Section  4.6,  an  adaptive 
CDMA  system,  including  joint  channel  parameters  estimation,  spreading  sequence  adap- 
tation, and  adaptive  MMSE  signal  reception,  is  presented  to  facilitate  the  implementation 
of  the  minimum-ETSC  sequences  in  real-life  communication  systems.  Numerical  exam- 
ples are  presented  to  examine  the  performance  of  the  proposed  joint  adaptation  structure. 
Proofs  of  the  propositions  and  convergence  of  the  proposed  algorithms  are  provided  in 
Section  4.7.  Conclusions  are  drawn  in  Section  4.8. 
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4.1  System  Model 

We  consider  a synchronous  CDMA  forward-link  transmission  and  assume  that  the 
signals  of  different  users  undergo  independent  slow  flat  fading,  that  is,  Hfc  = ck  for  k = 

1,2  ,■■■  ,K.  With  these  assumptions,  the  general  signal  model  described  in  Chapter  2 can 
be  simplified.  There  are  N samples  for  every  data  symbol  at  the  output  of  the  chip-matched 
filter,  which  can  be  arranged  into  an  A-dimensional  column  vector  r^: 

K 

Tk  = ck  S bon)3^n  + rifc,  (4.1) 

m— 1 

where  Ck  denotes  the  channel  gain  from  the  base-station  to  the  kth  user’s  receiver,  and  nk 
represents  the  contribution  from  the  AWGN. 

The  vector  rk  is  fed  into  a linear  FIR  filter  with  weights  The  decision  statistic 
Zk  = w^r/c  is  hard-limited  to  obtain  the  estimate  of  the  Oth  symbol  of  the  /cth  user. 

4.2  Spreading  Sequence  Allocation  with  MMSE  Signal  Reception 

In  this  section,  we  propose  a transmission  power  and  spreading  sequence  allocation 
scheme  to  support  all  users  with  a uniform  target  SIR  in  the  forward-link  CDMA  system 
described  in  Section  4.1  when  MMSE  receivers  are  employed  for  signal  reception. 

4.2.1  MMSE  Signal  Reception  and  User  Capacity 

The  optimal  MMSE  reciever  weights  are  chosen  as 

wfc  = — R^a*,  (4.2) 

Ck  k 

where  Rrfc  is  the  effective  total  correlation  matrix  observed  by  the  kth  user  given  by 

R-t*.  = -jEjzfcZj^]  = AtAj>  + fjkI.  (4.3) 

Ck 

In  above,  AT  = [aj  a2  • ■ • a.K]  is  an  N x K matrix  formed  by  grouping  all  spreading 
sequence  vectors  and  fjk  = % is  the  effective  noise  density  at  the  kth  user’s  receiver.  With 
the  optimal  MMSE  weights,  the  receiver  output  SIR  is  given  by 


SIR*  = afR/afc 


(4.4) 
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where  R*  is  the  effective  noise-plus-interference  correlation  matrix  defined  by 


Rfc  — R-rfc  — a*ajf . 


(4.5) 


An  equivalent  expression  of  the  SIR  using  the  matrix  inversion  formula  can  be  derived  as 
well: 


SIR* 


= afR^a*. 


(4.6) 


1 + SIR* 

We  assume  that  all  users  are  to  be  supported  with  a uniform  target  SIR  7.  This 
implies  af  R^a*  > 7,  or  equivalently,  afR^a*  > ^77*  for  1 < k < K.  Due  to  limited 
bandwidth,  there  may  not  be  a feasible  allocation  to  support  users  if  there  are  too  many 
users  in  the  system  or  the  SIR  targets  of  the  users  are  too  high.  Similar  to  the  case  of 
AWGN  channel,  we  define  the  system  user  capacity  as  the  maximum  number  of  users  the 
system  can  support  with  the  target  SIR.  It  is  proved  by  Viswanath  et  al.  [50]  that  a necessary 
and  sufficient  condition  for  the  maximum  number  of  users  K that  can  be  supported  in  the 
forward-link  CDMA  system  with  spreading  gain  N and  SIR  target  7 is 


Ke{ 7)  < N,  (4.7) 

where  e(7)  = ^77  as  defined  in  previous  sections.  There  is  a similar  interpretation  of  this 
user  capacity  constraint  in  terms  of  effective  bandwidth  e{pf)  as  in  the  AWGN  channel, 
that  is,  all  K users  can  be  supported  in  the  system  if  and  only  if  the  sum  of  their  effective 
bandwidth  is  less  than  the  spreading  gain  of  the  system. 

4.2.2  Minimum-ETSC  Allocation  Scheme 

The  transmission  power  of  the  kth  user  is  denoted  as  pk  = ||a*||2,  for  1 < k < K. 
For  the  case  of  K < N,  the  constraint  in  (4.7)  is  always  satisfied  and  we  can  always 
allocate  orthogonal  sequences  to  the  users.  With  this  allocation,  each  user  does  not  suffer 
interference  from  other  users.  The  transmission  power  of  the  kth  user,  for  1 < k < K,  can 
be  set  as 


Pk  = IPk 


(4.8) 
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to  satisfy  the  SIR  target  7.  The  total  transmission  power  is  7 J^k=i  Vk-  Our  main  interest 
lies  in  the  case  when  the  number  of  users  is  larger  than  spreading  gain,  that  is,  K > N. 

A simple  analysis  of  the  structure  of  effective  total  correlation  matrices 

Rr*  = ATA^  + fjk I,  k — 1, 2,  • • • , K 

indicates  that  the  contribution  from  the  spreading  sequences  ATA%  is  shared  by  all  the 
correlation  matrices.  If  we  denote  eigenvectors  of  the  correlation  matrices  as  the  “ channel 
space ” spanned  by  spreading  sequences  in  the  CDMA  system,  eigenvalue  decomposition 
of  the  effective  correlation  matrices 


RTk  = UAU*  + fjkI,  k = 1,2,  - ■ ■ , K 

indicates  that  all  correlation  matrices  have  the  same  eigenvectors  and  thus  share  the  same 
“channel  space”.  In  general,  the  effective  noise  densities  of  users  may  be  widely  spread. 
Following  the  idea  of  user  classification  in  Viswanath  et  al.  [50],  we  classify  all  users  into 
different  classes,  and  employ  different  power  and  sequence  allocation  strategies  to  these 
different  classes.  It  is  intuitively  reasonable  to  allocate  orthogonal  “channels”  to  users 
whose  channel  conditions  are  much  worse  than  those  of  the  other  users,  and  let  the  other 
users  share  the  rest  of  the  “channel  space”,  so  that  the  performance  of  the  whole  system 
would  not  be  ruined  by  those  users  with  poor  channel  conditions.  Based  on  this  obser- 
vation, we  present,  below,  a power  control  and  spreading  sequences  allocation  scheme 
to  support  all  users  with  a uniform  target  SIR  in  the  forward  link  of  the  CDMA  system 
modeled  by  (4.1)  when  the  condition  in  (4.7)  is  satisfied. 


Proposition  4.2.1  Without  loss  of  generality,  assume  a descending  order  of  the  effective 
noise  densities,  that  is,  fji  > > • • • > fjK-  Define  fjo  = 00.  There  exists  a unique 

number  0 < k*  < N(1  + j)  — Ky  such  that 


Vk’+l 


< 


Yji=k*+i  hi 


N(  1 + 7)  — K 7 — k* 


< hk*- 


(4.9) 
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We  denote  the  first  k*  users,  that  is,  Users  1,  2,  • • • , k*  as  overfaded  users,  and  Users  k*  + 
1,  k*  + 2,  • • • , K as  nonoverfaded  users.  When  k*  = 0,  it  means  that  there  is  no  overfaded 
user. 

Now  we  introduce  a power  and  sequence  allocation  scheme  based  on  this  classifi- 
cation of  users.  Overfaded  users  and  nonoverfaded  users  are  applied  different  stratigies. 
Transmission  power  for  the  users  can  be  allocated  as 


With  this  power  allocation,  we  note  that  the  k*  overfaded  users  should  be  allocated  with 
k*  orthogonal  “channels”  to  satisfy  the  SIR  requirement,  which  can  be  done  as  described 
below.  Generate  an  arbitrary  N x N unitary  matrix  U = [ui,  U2,  • • ■ , U/v]  with  column 
vectors  u„  for  n = 1,  2,  • • • , N.  The  spreading  sequence  of  the  Arth  user,  for  1 < k < k* , 
is  assigned  as 


where  pk  is  the  transmission  power  in  (4.10).  The  remaining  K — k * nonoverfaded  users 
share  the  remaining  N - k*  channels,  that  is,  the  N - k*  dimensional  subspace  spanned  by 
Ufc«+i,  ufc*+2,  • • • , u.v  When  K < N,  all  users  are  overfaded  and  are  allocated  orthogonal 
“channels”.  This  result  is  consistent  with  the  general  orthogonal  sequence  allocation  when 
K < N.  Propostion  4.2.2,  whose  proof  is  provided  in  Section  4.7  justifies  the  existence  of 
the  spreading  sequences  of  the  nonoverfaded  users,  and  gives  one  method  of  constructing 
the  sequences  using  the  results  of  theory  of  majorization  [70], 


for  k = 1,  2,  • • • , k*, 

for  k = k*  + 1,  k*  + 2,  ■ • • , K. 


3-k  y/Pk^-kj 


(4.11) 


Proposition 


we  can  construct  an  ( N — k*)  x 


(K  - k*)  real  matrix  A = [afc.+i,  afc.+2,  • • ■ , a K]  that  is  characterized  by 


diag{ArA)x  = [pfc.+1,  • • • ,Pk]T  ■ 


(4.12) 
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We  note  that  A can  be  constructed  recursively  using  the  algorithm  of  Lok  and  Wong  [52]. 
The  spreading  sequence  a*;  of  the  /cth  user,  for  A;  = fc*  + l,  /c*+2,  ■ • • , K,  can  be  constructed 
by  padding  k*  zeros  at  the  beginning  of  the  column  vector  and  then  left  multiplying  by 
the  unitary  matrix  U, 


ak  = U[0  — OalY,  k = k*  + l,k*  + 2,---,K. 

k*  zeros 

With  the  proposed  allocation  scheme,  the  total  transmission  power  is 


(4.13) 


K 


N(  1 + 7)  — A7  — k* 


Pgwbe  — ^2/  m — 7 ^2  Pl  + 


k= 1 


Z=1 


(4.14) 


The  spreading  sequence  matrix  is 


A t — [a-i5  • ' ' ; afc*i  afc*+i)  ■ ' • 5 a/r] 

VpT 


= U 


0 

yJPk * 

0 A 


(4.15) 


and  the  effective  total  correlation  matrix  Rr,  observed  at  the  fcth  user’s  receiver  is 


— AtAj,  + 77*1 


U 


7m 


0 


7m* 


0 


EI'C  - 

i=fc*+l  Vi  I 


+ 77*1. 


(4.16) 


ZV(l+7)  — K'y— k* 

Recall  that  the  SIR  constraint  of  the  hi h user  is  a^ R^1  a^  > e(7).  We  proceed  to 
show  that  all  users  satisfy  the  SIR  constraint  with  this  allocation  scheme: 


(i)  For  k = 1,2, 


afR^afc 


7 Vk 

7 Vk  + Vk 

= e(7). 


(4.17) 


'We  use  the  notation  diag{AT  A)  to  indicate  the  main  diagonal  of  square  matrix  ATA. 
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(ii)  For  k — k*  + 1,  k*  + 2,  • • • , K, 


afR^a*  = a[(AAJ’  + %I)  lhk 

Pk 


= e(7). 


(4.18) 


Proposition  4.2.3  below  indicates  that  the  sequences  constructed  above  minimize 
the  extended  total  squared  correlation  (ETSC)  of  the  unequal  power  spreading  sequences, 
which  is  defined  as 


Proposition  4.2.3  When  K < N,  given  the  transmission  powers  pk,  for  1 < k < K,  as 
defined  in  (4.8),  orthogonal  spreading  sequences  give  the  minimum  ETSC.  When  K > N, 
given  the  transmission  powers  pk,  for  \ < k < K,  as  defined  in  (4.10),  the  spreading 
sequence  set  constructed  in  (4.15)  is  the  set  of  sequences  that  minimizes  the  ETSC. 

Hence,  we  call  this  allocation  scheme  as  minimum-ETSC  scheme. 

4.2.3  Properties  of  the  Minimum-ETSC  Scheme 

There  are  some  observations  pertaining  to  this  well-structured  allocation  scheme. 
When  the  effective  noise  densities  of  users  are  not  widely  spread,  there  will  not  be  any 
overfaded  user,  that  is,  k*  = 0.  In  this  case,  all  K users  share  the  N channels  and  the 
spreading  sequence  matrix  is  characterized  by 


K K 


(4.19) 


k= l 1=1 


(4.20) 


where 


and 
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Due  to  the  characteristic  of  the  matrix  ArAp  = ppT.BE  I,  the  ETSC  of  the  unequal  power 


spreading  sequences  satisfies  the  Welch’s  lower  bound  [48]  with  equality,  that  is, 

K K rs—iK  n_  1 1 2 \ 2 


EE  iia^a'ii2  = 


(£f=il|a*ll2)s 


^=1  (=i 


N 


Following  Viswanath  et  al.  [50],  we  call  this  set  of  spreading  sequences  as  generalized 
WBE  (GWBE)  sequences  with  different  powers,  to  distinguish  from  the  WBE  sequence 
set  (defined  in  Massy  [47]),  in  which  all  sequences  are  of  unit  power.  With  the  GWBE 
sequences,  all  the  N eigenvalues  of  the  matrix  ATAP  are  the  same.  This  means  that  all 
the  “channels”  in  the  CDMA  system  are  evenly  occupied  and  equally  “crowded”.  When 
the  effective  noise  densities  of  all  users  are  equal,  that  is,  fji  — ■ • • = t}k  — fj,  the  channel 
model  reduces  to  an  AWGN  channel.  It  is  easy  to  show  that  the  set  of  spreading  sequences 
constructed  above  reduces  to  the  set  of  optimal  spreading  sequences  for  an  AWGN  channel 
given  in  Chapter  3. 

Using  this  well-structured  power  and  sequence  allocation  scheme,  the  MMSE  re- 
ceiver weight  vectors  can  be  easily  obtained: 


(i)  For  k = 1,2,---,  A:*, 


W,  = 


— u(a  + f}kiyl\jH  ^PkUk 
1 

Cjt(l  + 7 )fjkak' 


(4.21) 


(ii)  For  k — k*  + 1,  k*  + 2,  ■ • • , K, 


W £ — R- 71.  tu 

Cfc  k 

1 


-U(A  + 7pr1UffU[0---0a£; 

Cfc 

1 

f 7Ei=t»+ im_  , - \ 

0fcVjV(l+7)-iC7-fc*  ^ 'Ik) 


(4.22) 
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We  notice  that  the  MMSE  receiver  weights  are  the  same  as  the  matched-filter  re- 
ceivers (up  to  a scaling  factor).  This  property  is  important  in  the  sense  that  it  makes  the 
implementation  of  the  receiver  much  easier,  that  is,  using  the  simple  matched  filter,  the 
system  can  achieve  the  performance  of  the  optimal  MMSE  receiver. 

Finally,  we  notice  that  the  choice  of  the  spreading  sequences  is  not  unique  (although 
the  power  allocation  scheme  is  fixed).  For  example,  any  unitary  rotation  of  the  spreading 
sequence  set  has  the  same  eigenvalue  distribution  and  is  still  a valid  allocation  to  support 
all  users.  If  the  unitary  matrix  U is  chosen  to  be  an  identity  matrix,  the  system  becomes 
a mixed  TDMA-CDMA  system.  Elementary  vectors  are  assigned  as  spreading  sequences 
of  the  overfaded  users.  As  a result  of  this  allocation,  the  signal  of  the  kth  overfaded  users 
is  transmitted  only  during  the  &;th  chip  interval  without  spread.  Nonoverfaded  users  do  not 
transmit  during  the  first  k*  chip  intervals,  and  the  actual  spreading  gain  of  the  nonoverfaded 
users  reduces  to  N — k*. 

4.3  Spreading  Sequence  Allocation  with  Matched  Filter  Reception 

An  important  feature  of  the  class  of  minimum-ETSC  sequences  proposed  in  Sec- 
tion 4.2  is  that  the  MMSE  receiver  weights  are  the  same  as  the  matched-filter  receiver 
weights.  In  this  section,  we  propose  another  power  and  sequence  allocation  scheme,  which 
can  support  users  when  matched-filter  receivers  are  employed. 

4.3.1  Matched  Filter  Signal  Reception  and  User  Capacity 

Signals  of  all  K users  in  the  system  are  detected  simultaneously  and  independently 
using  the  matched-filter  receiver.  The  data  stream  of  the  £;th  user,  for  1 < k < K,  is 
allocated  power  pk  and  is  spread  with  the  unit-energy  signature  sequence  sk.  Here,  we 
normalize  the  power  of  the  spreading  sequences,  and  denote  sk  = as  the  unit-energy 
signature  sequence  of  the  klh  user  to  distinguish  it  from  the  spreading  sequence  All 
users  are  supported  with  a uniform  target  SIR  7.  This  implies  that  the  output  SIR  at  the  kth 
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user’s  MF  receiver,  for  1 < k < K,  satisfies 

orp  ckPk P 'k ^ (4  21) 

'k  Ei^kclpi\sksi\2 + Vk  Ei^kPi\sksi\2  + Vk  ~ 

where  Ck  denotes  the  channel  gain  from  the  base-station  to  the  fcth  user’s  receiver,  77*  is 

the  power  spectral  density  (PSD)  of  the  thermal  noise,  and  fjk  — ^ f is  the  effective  noise 

density  at  the  kth  user’s  receiver. 

The  SIR  constraints  in  (4.23)  can  be  equivalently  expressed  as 


Pk  > e(j) 


' K 

^PlC-kl  + Pk  , 
.1=1 


(4.24) 


where  e(7)  = “ and  Cki  — |Sjfs;|2.  Because  is  of  unit  energy,  we  have  Ckk  = 1 and 
0 < Cki  < 1 for  k ^ l.  Written  in  matrix  form,  we  have 


p > 2e(7)(Cp  + f7),e(7)(Cp  + f7),  (4.25) 


where  p = [px  p2  ■ ■ -Pk]t,  V = [.Vi  V2  • • -fj k]t  are  K -dimensional  column  vectors,  and 
C = (cki)  is  a K x K square  symmetric  nonnegative  matrix.  All  K users  are  supported 
with  the  target  SIR  7 if  (4.25)  can  be  solved  with  equality 


1 

Al) 


I- 


C 


P = V 


(4.26) 


Again,  due  to  limited  system  capacity,  there  is  a user  capacity  constraint  of  the 
system,  which  is  the  maximum  number  of  users  that  can  be  supported  with  the  uniform 
target  SIR  7 for  all  the  possible  allocation  of  sequence  power  and  signature  sequences 
when  matched-filter  receivers  are  employed. 


Proposition  4.3.1  Given  spreading  gain  N,  and  a uniform  SIR  target  7,  the  maximum 
number  of  users  the  system  can  support  satisfy: 


Ke( 7)  < N. 


(4.27) 


2Inequality  in  (4.25)  means  the  elementwise  inequality  between  the  vectors. 
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We  notice  that  the  system  user  capacity  when  employing  the  matched-filter  re- 
ceiver is  the  same  as  that  when  employing  MMSE  receiver.  This  indicates  that  there  are 
valid  schemes  to  allocate  the  signature  sequences  and  transmission  power  for  the  users  to 
achieve  the  target  SIR  using  either  matched-filter  receiver  or  MMSE  receiver  when  the 
user  capacity  constraints  in  (4.27)  is  satisfied.  However,  different  allocation  schemes  em- 
ploy different  set  of  signature  sequences  and  consume  different  transmission  power.  The 
minimum-ETSC  scheme  proposed  in  Section  4.2  is  a valid  allocation  scheme  for  both  types 
of  recivers,  since  MMSE  receiver  weights  are  the  same  as  matched-filter  receiver  weights 
of  each  user  with  the  minimum-ETSC  scheme. 

4.3.2  Minimum-TSC  Allocation  Scheme 


In  synchronous  CDMA  systems,  orthogonal  sequences  are  generally  employed 
when  K < N.  Our  main  interest  lies  in  the  case  of  N < K < A reasonable  de- 
sign procedure  is  to  start  with  the  construction  of  a set  of  unit-energy  signature  sequences, 
and  then  proceed  to  obtain  the  minimum  user  powers  required  to  support  all  users  using 
this  signature  set  of  sequences.  Given  a choice  of  signature  sequences,  the  maximum  tar- 
get SIR  that  can  be  achieved  by  all  users  is  restricted  by  the  largest  eigenvalue  of  C,  as 
indicated  by  Lemma  4.7.5.  In  other  words,  the  choice  of  unit-energy  signature  sequences 
is  essential  to  support  users  even  when  the  user  capacity  constraint  in  (4.27)  is  satisfied. 
With  valid  choice  of  signature  sequences,  the  minimum  power  required  to  support  users  is 
given  by 


P = 


1 

_e(7) 


-i  -1 


I-C 


V- 


(4.28) 


Hence,  different  choices  of  signature  sequences  require  different  transmission  power.  Given 
a target  SIR  7,  and  the  total  number  of  users  K in  the  system,  the  structure  of  the  opti- 
mal set  of  signature  sequences,  which  resuires  a minimum  amount  of  total  transmission 
power  among  all  possible  choice  of  signature  sequences,  is  currently  unknown.  From  the 
proof  of  the  user  capacity,  which  is  provided  in  Section  4.7,  we  know  that  the  set  of  unit- 
energy  WBE  sequences  is  one  valid  choice  of  signature  sequences  when  the  user  capacity 
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constraint  is  satisfied.  The  WBE  sequences  are  characterized  by  STS^  = ^1,  where 
Sr  = [si,  s2,  • • • , s k]  is  an  ./V  x K matrix  formed  by  grouping  all  K unit-energy  signa- 
ture sequence  vectors.  The  total  transmission  power  required  to  support  users  with  the 
target  SIR  7 when  the  set  of  unit-energy  WBE  sequences  is  employed  is  given  in  Proposi- 
tion 4.3.2,  the  proof  of  which  is  provided  in  Section  4.7. 


Proposition  4.3.2  The  minimum  total  transmission  power  required  when  employing  the 
WBE  sequences  as  signature  sequences  is 

Ne(7)  £f=i  fjk 


K 


P 


WBE 


= 


k=l 


N-Ke{i)  ' 


(4.29) 


To  distinguish  from  the  minimum-ETSC  scheme,  we  call  this  allocation  scheme  as 
minimum-TSC  scheme,  since  the  set  of  unit-energy  WBE  sequences  minimizes  the  total 
squared  correlation  (TSC)  that  is  defined  as 

K K 

TSC  = EE  |s"s;|2.  (4.30) 

fc=i  /=i 

WBE  sequences  can  be  constructed  iteratively  using  the  TSC  reduction  algorithm  of  Ulukus 
[51]. 

4.4  Power  Efficiency  of  the  Proposed  Sequence  Allocation  Schemes 

High  power  efficiency  is  one  of  the  main  design  concerns  because  transmission 
power  is  a precious  resource  in  wireless  systems.  Also,  reduction  of  the  transmission 
power  in  one  cell  can  reduce  the  co-channel  interference  (CCI)  to  users  in  other  cells.  In 
this  section,  we  try  to  compare  and  gauge  the  performance  of  the  two  power  and  sequence 
allocation  schemes  described  previously  in  terms  of  the  total  transmission  power  required 
by  each  of  the  schemes. 

4.4.1  Lower  Bound  on  the  Transmission  Power 

In  the  absence  of  MAI,  the  minimum  transmission  power  of  the  kth  user  to  meet 
the  SIR  requirement  is  pk  = jqk.  It  is  straightforward  that  the  user  needs  more  power  in 
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the  existence  of  MAI,  that  is,  pk  > 7 fjk.  Summing  up  the  lower  bounds  of  transmission 
powers  of  all  users,  we  obtain  a lower  bound  of  the  total  transmission  power: 


Pgwbe  > 7Ef=i^  - pB, 
Pwr  r > 7 Vk  = Ps 


(4.31) 


'WBE  ET  7 Z^fc=l  'Ik  — rB- 
Because  of  the  possible  large  spread  of  the  effective  noise  densities  of  all  users  and  the 
“crowdedness”  of  the  system,  this  trivial  lower  bound  may  be  loose. 

4.4.2  Minimum-ETSC  Scheme  vs  Minimum-TSC  Scheme 


A simple  algebra  derivation  shows  that 

>-  + 


NlEhm  ^AA  , (AT -fb  £?=*•+!% 


AT(l+7)  — Ky  ' Ar(l  + 7)  - iCy  — k* 


(4.32) 


with  equality  when  k*  = 0.  The  inequality  in  (4.32)  shows  that  the  minimum-ETSC 
scheme  is  more  power  efficient  than  the  minimum-TSC  scheme. 

4.4.3  Lagrangian  Search  of  Minimum  Power  with  MF  Receiver 

Matched  filter  can  be  employed  for  signal  reception  in  both  the  minimum-ETSC 
and  minimum-TSC  allocation  schemes.  We  seek  the  minimum  total  transmission  power 
required  to  support  users  when  using  matched-filter  receiver  numerically  to  gauge  their 
power  efficiencies.  We  consider  the  Lagrangian  multiplier  method  similar  to  the  one  in  Lok 
and  Wong  [52]  to  minimize  the  transmission  power  and  incorporate  the  SIR  requirements 
as  a penalty  function.  The  Lagrangian  function  L is  formed  as  follows: 


K 


K 


£ = X>fa,  + y>  (t-= , if1  a')‘  . „ ] , (4.33) 

A tt  V E am)2  + i?1afa, 

where  A;,  for  l = 1,  2,  • • • , K are  Lagrange  multipliers.  Given  Given  spreading  gain  N, 
number  of  users  K,  target  SIR  7,  and  effective  noise  densities  77/,  for  1 < l < K,  the 
Lagrangian  function  L is  a function  of  A / and  components  of  a;  for  1 < l < K.  We 
consider  a gradient  search  approach  to  seek  a stationary  point  of  the  Lagrange  function. 

The  derivatives  of  L with  respect  to  A^  and  are  given  by 
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(4.34) 


At  each  iteration  of  the  gradient  search,  a^  is  updated  using  a gradient  descent  algorithm 
and  Afc  is  updated  using  a gradient  ascent  algorithm 


j afc  ^ afc  ^ , A,'  1 , 2 , • • • , K ", 

\ A*;  «—  Afc  + , A;  = 1, 2,  • ■ • , FT, 


(4.35) 


where  ua  is  the  step  size  of  updating  the  spreading  sequences  and  u\  is  the  step  size  of 
updating  the  Lagrange  multipliers.  The  gradient  descent  algorithm  may  converge  to  local 
minima.  Only  the  local  minima  that  meet  the  target  SIR  of  all  users  are  treated  as  valid 
searching  results.  After  convergence,  the  total  transmission  power  can  be  calculated  as 

PLag  = Ef=l  Ikll2- 

We  consider  a sample  system  with  spreading  gain  N = 7 and  SIR  target  7 = 4dB. 
The  user  capacity  of  this  sample  system  is  9 users  according  to  (4.27).  We  conduct  two 
simulation  examples  with  K = 8 and  K — 9 to  evaluate  the  power  efficiencies  of  the  two 
allocation  schemes.  The  inverse  of  effective  noise  densities  J-,  for  1 < k < K,  are  gen- 

7lk 

erated  independent  exponential  distributions  with  E[^]  = lOdB.  For  each  example,  we 
conduct  100  independent  realizations  and  compare  the  best  Lagrangian  multiplier  search- 
ing results  PLag  to  the  total  transmission  power  Pwbe  required  by  the  minimum-TSC  al- 
location scheme  and  the  total  transmission  power  Pgwbe  required  by  the  minimum-ETSC 
allocation  scheme.  For  ease  of  comparison,  we  normalize  Pgwbe  and  Pwbe  with  PLag- 
Simulation  results  are  shown  in  Figs.  4.1  and  4.2.  We  can  observe  that  the  minimum-ETSC 
scheme  is  more  power-efficient  than  the  minimum-TSC  scheme,  which  is  consistent  with 
the  analysis  in  (4.32).  The  power  obtained  using  the  Lagrangian  searching  method  is  very 
close  to  Pgwbe  for  all  the  realizations  conducted. 


'd/d 
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Figure  4. 1 : Power  efficiency  with  MF  receiver:  N=7,  K=8,  7 = 4dB 


Figure  4.2:  Power  efficiency  with  MF  receiver:  N=7,  K=9,  7 = 4dB 
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4.4.4  Lagrangian  Search  of  Minimum  Power  with  MMSE  Receiver 

The  minimum-ETSC  scheme  is  designed  for  MMSE  signal  reception,  although  it 
turns  out  that  the  MMSE  receiver  weights  are  the  same  as  those  of  the  matched-filter  re- 
ceiver. Again,  we  consider  a Lagrangian  multiplier  method  to  seek  the  minimum  total 
transmission  power  required  numerically  to  gauge  the  power  efficiency  of  the  minimum- 
ETSC  scheme  when  using  MMSE  receiver.  The  Lagrangian  function  L is  formed  as  fol- 
lows: 

£ = E IWI2  + E A<h  - a/'E'a,)2,  (436) 

Z=1  l—l 

where  A i,  for  1 < l < K,  are  Lagrange  multipliers.  The  Lagrangian  function  L is  a func- 
tion of  A i and  components  of  a/  for  1 < l < K,  given  the  system  and  channel  parameters 
including  the  spreading  gain  N,  number  of  users  K , target  SIR  7,  and  effective  noise  den- 
sities rji,  for  1 < l < K.  We  consider  a gradient  search  approach  to  seek  a stationary  point 
of  the  Lagrangian  function.  The  derivative  of  L with  respect  to  Xk  is 

dL 


d\ k 


(7  ~ afc  R^  afe)  , k = 1,2,  ■ ■ ■ , K. 


(4.37) 


-1 


The  derivative  of  L with  respect  to  ak  can  be  obtained  as  follows.  First  notice  that  Rfc 
does  not  include  information  of  ak,  while  R;_1  for  l ^ k can  be  expressed  explicitly  as  a 
function  of  ak  via  the  matrix  inversion  formula: 


\-i 


R;  1 = (ATAf  - a,af  + ^1)' 

= (ArA"  - a,af  - afcaf  + fj{l  +afcaf  )_1 
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1 + aj!  Kikak 


(4.38) 


Therefore,  the  derivative  of  L with  respect  to  ak  is  given  by 
dL 


dak 


2afc  -4Afc(7-afRit1a)t)Rfc1a/t  + 4^A/(7-a/HR,  1al) 

l^k 


az//Rifclaz  ft- 1_  !az//RZfclafc|2  n-l_ 
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1 +az  Rlk  az 


(1  + af  Rifc  ak)2 


(4.39) 
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At  each  iteration  of  the  gradient  search,  a.k  and  \k  are  updated  using  the  algorithm  in  (4.35). 
After  convergence,  the  total  transmission  power  can  be  calculated  as  PLag  = J2k= 1 llafcl|2- 

Using  the  same  sample  system  as  in  Section  4.4.3,  we  present  three  numerical  ex- 
amples to  examine  the  power  efficiency  of  the  minimum-ETSC  allocation  scheme.  The 
user  capacity  of  this  sample  system  is  9 users  per  cell  according  to  (4.7).  We  conduct  three 
simulation  example  with  K = 7,  8,  9,  which  represent  three  cases  when  the  number  of 
users  is  smaller  than  or  equal  to  the  spreading  gain,  the  number  of  users  is  larger  than  the 
spreading  gain,  and  the  system  is  maximally  loaded,  respectively.  For  all  three  cases,  we 
conduct  100  independent  realizations  and  compare  the  best  Lagrangian  multiplier  search- 
ing results  to  the  total  transmission  power  Pgwbe  given  in  (4.14).  For  the  case  of  K < N, 
as  shown  in  Fig.  4.3,  Fagrangian  search  always  gives  the  orthogonal  sequences,  which  are 
the  optimal  sequences  for  this  case.  For  the  two  cases  of  K > N , as  shown  in  Figs.  4.4  and 
4.5,  the  Fagrangian  searching  results  indicate  the  existence  of  spreading  sequences  that 
satisfy  the  SIR  requirement  with  less  total  transmission  power  than  the  power  required  by 
the  minimum-ETSC  scheme.  The  trivial  lower  bound  is  also  presented  in  Figs.  4.3,  4.4, 
and  4.5,  from  which  we  can  observe  that  the  lower  bound  becomes  very  loose  when  the 
system  is  heavily  loaded. 

Both  the  theoretical  analysis  and  Fagrange  multiplier  searching  results  indicates 
that  the  minimum-ETSC  scheme  provides  better  power  efficiency  than  the  minimum-TSC 
scheme.  It  appears  that  the  minimum-ETSC  scheme  is  optimal  in  terms  of  requiring  the 
minimum  total  transmission  power  when  the  matched-filter  receiver  is  employed.  How- 
ever, analytical  verification  of  the  optimality  of  this  scheme  remains  open  at  this  moment. 
In  spite  of  the  observation  that  the  minimum-ETSC  scheme  is  not  optimal  in  terms  of  min- 
imizing the  total  transmission  power  when  K > N,  results  from  experiments  under  dif- 
ferent channel  conditions  indicate  that  the  minimum-ETSC  scheme  only  consume  slightly 
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ure  4.3:  Power  efficiency  with  MMSE  receiver:  N=7,  K=7,  7=4  dB 


realizations 


Figure  4.4:  Power  efficiency  with  MMSE  receiver:  N=7,  K=8,  7=4  dB 
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realizations 


Figure  4.5:  Power  efficiency  with  MMSE  receiver:  N=7,  K=9,  7=4  dB 

more  power  than  the  best  Lagrangian  searching  results  except  for  when  the  system  is  max- 
imally loaded.  In  following  sections,  we  only  focus  on  the  application  of  the  minimum- 
ETSC  allocation  scheme  in  CDMA  forward-link.  In  terms  of  application,  an  important 
advantage  of  the  minimum-ETSC  sequence  set  over  the  one  searched  by  Lagrangian  mul- 
tiplier method  after  many  iterations  is  that  the  set  of  minimum-ETSC  sequences  can  be 
easily  constructed  using  algorithms  that  are  simple  and  amenable  to  adaptive  implementa- 
tion. We  are  going  to  present  these  algorithms  in  the  next  section. 

4.5  Iterative  Construction  of  the  Minimium-ETSC  Sequence 

Given  the  effective  noise  densities  of  the  users  in  the  system,  the  matrix  Ay  = 
[ai,  a2,  • ■ • , a k\  can  be  constructed  using  the  procedure  in  Section  4.2.  This  construction 
is  not  feasible  for  practical  implementation,  especially  when  the  channel  conditions  change 
with  time.  From  the  analysis  above,  we  notice  two  important  features  of  the  sequences  con- 
structed in  Section  4.2.  First,  the  sequences  belong  to  a class  of  sequences  that  minimize 
the  extended  total  squared  correlation  (ETSC).  Second,  the  MMSE  receiver  weights  are 
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the  same  as  those  of  the  matched  filter  up  to  a scaling  factor  for  this  class  of  sequences. 
Making  use  of  these  two  observations,  we  introduce  adaptive  algorithms  to  construct  the 
spreading  sequences  iteratively. 

4.5.1  An  Iterative  ETSC  Reduction  Algorithm 

Different  from  the  TSC  reduction  algorithm  of  Ulukus  [51],  which  assumes  equal 
power  of  all  users,  the  algorithm  we  introduce  here  assumes  that  each  user  has  differ- 
ent effective  noise  density  and  is  allocated  different  transmission  power.  Starting  with 
random  spreading  sequences  AT[0]  = [ajO],  a2[0],  • • • , a^-[0]]  with  the  predefined  power 
allocation  in  (4. 10),  the  algorithm  replaces  the  spreading  sequences  by  their  corresponding 
MMSE  receiver  weight  vectors  (up  to  a scaling  factor  to  guarantee  fixed  transmission  pow- 
ers) at  each  iteration.  Within  each  iteration,  spreading  sequences  of  all  users  are  updated 
sequentially  starting  from  the  first  user  as  follows: 


Algorithm  4.5.1  At  the  jth  iteration,  for  the  kth  user,  update  the  effective  total  correlation 
matrix, 


AjjJy]  — [ai[j  -1- 1],  • ■ • , &k-i[j  + 1])  afe[y],  ■ • • , a^L?’]], 

R-iy  [j]  = A Tk  [j]^Tk  \j\  + Vkl- 


Update  the  sequence  by 


a k[j  + 1]  = 


Pk 

i-2r 


r&nljWlj]- 


(4.41) 


An  equivalent  expression  of  Eq.  (4.41),  which  makes  the  analysis  of  convergence 
of  the  algorithm  easier,  is 


a k[j  + 1]  = 


Pk 


,#r 


mrijWij] 


Rfc  1 [;>*[?], 


(4.42) 


where  R k[j]  — KTk [j]  - a.k[j]a.j?[j].  The  proof  of  the  convergence  of  Algorithm  4.5.1  is 
provided  in  Section  4.7.  During  the  adaptation,  the  ETSC  among  the  unequal  power  se- 
quences decreases  monotonically  until  it  reaches  a fixed  point.  Calculation  of  the  transmis- 
sion powers  requires  the  exact  knowledge  of  the  effective  noise  densities  of  users.  Errors 
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in  the  estimates  of  these  channel  parameters  will  cause  inaccuracy  in  power  allocation  and 
hence  performance  degradation. 

4.5.2  Modified  Algorithm  with  Power  Adaptation 

A modified  approach  is  to  adjust  both  the  transmission  powers  and  sequences  of  all 
users  at  each  iteration. 

Algorithm  4.5.2  At  the  jth  iteration,  for  the  kth  user,  update  the  spreading  sequence  by 

a k[j  + 1]  = 9k[j]'Rfl[jW[j],  (4-43) 

where  the  coefficient  gk[j ] is  chosen  so  that 

af  [j  + 1]R^  [jW  [j  + 1]  = e(7) . (4.44) 


It  remains  open  to  analytically  establish  the  convergence  of  Algorithm  4.5.2.  Simulation 
results  show  that  Algorithm  4.5.2  always  converges  to  the  fixed  point  giving  the  minimum 
ETSC  starting  from  random  sequences. 

4.5.3  Convergence  Properties  of  Algorithms  4.5.1  and  4.5.2 

We  study  the  convergence  properties  of  the  proposed  adaptive  algorithms  via  com- 
puter simulations.  We  assume  the  same  sample  system  described  in  Section  4.4.  For  each 
algorithm,  we  conduct  two  simulation  examples  with  number  of  users  in  the  system  K = 7 
and  K = 9.  Simulation  results  using  Algorithm  4.5.1  are  shown  in  Figs.  4.6  and  4.7,  and 
simulation  results  using  Algorithm  4.5.2  are  shown  in  Figs.  4.8  and  4.9.  Simulation  results 
indicate  that  the  receiver  output  SIRs  of  all  users  converge  to  the  target  SIR  using  both 
algorithms  and  the  total  transmission  power  converge  to  the  total  power  Pgwbe  given  in 
(4.14)  using  Algorithm  4.5.2.  Eigenvalue  decomposition  of  the  spreading  sequences  after 
convergence  indicates  that  the  eigenvalues  always  converge  to  the  optimal  eigenvalue  dis- 
tributions given  in  (4.60),  that  is,  both  Algorithms  4.5.1  and  4.5.2  always  converge  to  the 
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Figure  4.6:  Algorithm  1:  N=7,  K=7,  7=4  dB 

fixed  point  giving  the  minimum  ETSC  starting  with  random  spreading  sequences.  We  also 
observe  that  Algorithm  4.5.1  converges  faster  than  Algorithm  4.5.2. 

4.6  Joint  Transmitter-Receiver  Adaptation  with  Channel  Estimation 

Algorithms  4.5.1  and  4.5.2  assume  explicit  knowledge  of  the  effective  noise  den- 
sities to  calculate  the  correlation  matrices,  and  assume  that  the  MMSE  receivers  always 
work  with  the  optimal  weights.  In  a practical  communication  system,  the  channel  con- 
ditions may  change  with  time.  It  is  impractical  to  estimate  the  channel  parameters  and 
construct  a new  set  of  spreading  sequences  every  time  the  channel  conditions  change.  It 
is  also  difficult  for  the  MMSE  receivers  to  adjust  their  filter  weights  when  a completely 
different  set  of  spreading  sequences  is  assigned  to  the  users.  In  order  to  solve  the  prob- 
lem practically,  we  develop  an  adaptive  structure  with  joint  channel  estimation,  transmitter 
power  and  sequence  adaptation,  and  adaptive  MMSE  signal  reception  in  the  forward-link 
CDMA  system. 
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Figure  4.7:  Algorithm  1:  N=7,  K=9,  7=4  dB 


4.6.1  Adaptive  Structure  and  Algorithms 

The  joint  adaptive  structure  illustrated  in  Fig.  3.1  can  be  employed  in  this  situation 
in  well.  Similar  to  the  joint  adaptation  with  channel  estimation  in  Chapter  3,  the  base- 
station  collects  the  feedback  information  from  the  receivers  and  estimates  the  effective 
noise  densities.  The  spreading  sequences  of  all  users  are  updated  at  the  base-station  in  a 
centralized  manner  and  then  employed  to  transmit  the  data  of  the  next  block.  The  iterative 
procedure  consists  of  several  steps  presented  below. 

Step  1.  The  base-station  transmits  the  jth  block  of  training  symbols  (contain  L 
data  bits)  of  the  kth  user  with  the  current  spreading  sequence  a k[j]  for  1 < k < K.  (The 
initial  spreading  sequence  Ar[0]  can  be  assigned  arbitrarily.) 

Step  2.  The  weight  vector  of  the  receiver  of  the  kth  user  is  adapted  using 
an  exponentially  weighted  RLS  algorithm,  for  1 < l < L,  as  was  done  in  Chapter  3.  At 
the  end  of  reception  of  the  jth  block  of  data,  the  MSE  for  the  jth  block  of  transmission  is 


receiver  output  SIR  (in  dB)  total  transmission  power  (normalized  by  P 
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Figure  4.8:  Algorithm  2:  N=7,  K=7,  7=4  dB 


receiver  output  SIR  (in  dB)  total  transmission  power  (normalized  by  P( 
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Figure  4.9:  Algorithm  2:  N=7,  K=9,  7=4  dB 
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estimated  as 

1 „ 

£k,j  = £ - WM,LZA:J,/|25  (4-45) 

and  this  estimates  is  sent  back  to  the  base-station.  Here.  Zk,j,i  is  the  chip-matched  filter 
output  vector  for  the  Ith  symbol,  and  w k,j,L  is  the  Lth  (last)  updated  weight  vector  at  the 
kth  receiver  within  the  j th  block  of  transmission. 

Step  3.  The  base-station  knows  the  current  spreading  sequences  A T[j]  and  ap- 
plies eigenvalue  decomposition  to  AT[j]A^[j]: 


AT[j}A»[j}  = U[j}A[j]XJH\j}. 


(4.46) 


Define  ak[j]  = UH[j]ak[j]  = [a^j,  . . . , a^,_ld]T.  Then  estimates  of  the  effective 

noise  densities  of  the  users  can  be  obtained  by  solving  the  following  set  of  equations, 
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(4.41) 


Step  4.  Using  estimates  i)k,  the  base-station  calculates  the  transmission  powers 
pk,  for  1 < k < K,  of  the  users  according  to  the  minimum-ETSC  scheme  proposed  in 
Section  4.2  and  updates  the  spreading  sequences  of  all  users  sequentially  starting  from  the 
first  user.  The  spreading  sequences  can  be  updated  with  both  Algorithms  4.5.1  and  4.5.2 
proposed  above.  For  the  kih  user,  its  spreading  sequence  is  updated  according  to 


3-k[j  + 1] 


Pk 


a k[j]Rrk  [ iW[j } 


-Rr^afcL)], 


(4.48) 


as  in  Algorithm  4.5.1,  or  according  to 


a k[j  + 1]  = gk&rlifakii], 


(4.49) 


where  the  coefficient  gk  is  chosen  so  that 

af  [j  + l]Rrfcl[j]afc[j  + 1]  = efr),  (4.50) 

as  in  Algorithm  4.5.2.  In  above,  R Tk[j]  = ArJj]AfJj]  + fjk  I,  and  A Tk[j]  = [ai[j  + 
1],  • • • , a k-i[j  + 1],  a.k[j],  • • • , a K\j]\-  At  the  end  of  the  jth  block  adaptation,  the  updated 
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spreading  sequences  A T[j  + 1]  = [ai [j  + 1],  a2[j  + 1],  • • • , a K[j  + 1]]  are  employed  to 
transmit  the  next  block  of  data.  The  adaptation  process  repeats  from  Steps  1 to  4 until  the 
training  period  finishes. 

If  continuous  updating  of  the  spreading  sequences  and  receiver  weights  are  desired 
after  the  training  period,  the  system  can  be  switched  into  a decision  feedback  mode.  Sym- 
bol decisions  made  by  the  receivers  are  employed  to  replace  the  training  symbols. 

4.6.2  Performance  Evaluation 

We  study  the  proposed  adaptive  system  through  computer  simulations.  With  the 
definition  of  the  receiver  output  SIR  in  (3.24),  the  receiver  output  SIR  of  the  Zth  sym- 
bol in  the  jth  block  of  the  kth  user’s  signal  can  be  calculated  by  using  (3.25),  in  which 
R k[j]  — Yhi+k  a; [j]  + The  total  transmission  power  P[j]  during  the  jth  block  of 
transmission  can  be  calculated  using  (3.26). 

A training  sequence  of  5000  symbols  for  each  user  is  sent  in  blocks  during  the  train- 
ing period.  A block  length  of  L — 100  symbols  is  employed  according  to  the  discussion  in 
Section  3.4.  We  study  the  performance  when  there  are  K = 9 users  in  the  system.  Starting 
from  random  sequences,  all  users  adapt  their  transmitters  and  receivers.  Figs.  4.10  and 
4.1 1 shows  the  receiver  output  SIRs  and  the  total  transmission  power  during  the  adaptation 
process.  It  is  shown  that  the  receiver  output  SIRs  of  all  users  converge  close  to  the  SIR 
target  and  the  total  transmission  power  approaches  Pqwbe  as  calculated  in  (4.14).  Due 
to  the  imperfect  estimates  of  channel  parameters,  there  are  oscillations  for  both  the  SIR 
and  the  transmission  power  curves.  Compare  the  convergence  using  Algorithms  4.5.1  and 
4.5.2,  we  observe  less  oscillations  when  using  Algorithm  4.5.2  and  there  is  no  substantial 
difference  on  convergence  speed  in  this  joint  transmitter-receiver  adaptation. 


receiver  output  SIR  (in  dB)  total  transmission  power  (normalized  by  P 


68 


Figure  4. 10:  Performance  of  joint  transmitter-receiver  adaptation:  Algorithm  1 


receiver  output  SIR  (in  dB)  total  transmission  power  (normalized  by  P 
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Figure  4.1 1:  Performance  of  joint  transmitter-receiver  adaptation:  Algorithm  2 
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4.7  Proof  of  Propositions  and  Convergence  of  Algorithm  4.5.1 

4.7.1  Proof  of  Proposition  4.2.1 

We  assume  a decreasing  order  of  the  effective  noise  densities,  that  is,  rji  > f/2  > 
• • • > f/K , and  prove  the  Proposition  4.2.1  for  two  cases  when  K < N and  K > N: 

4.7.1. 1 Case  1:  K < N 

When  K < N,  we  always  have  N{  1 + 7)  — K 7 — (k  — 1)  > K — (A;  — 1)  for 
k — 1,2 ,■■■  ,K.  So  the  inequality 

~ > 2^;=fc  7fc  

/fc  — A'  — (A;  — 1)  IV  (1  + 7)  — K'y  — (k  — 1) 
is  valid  for  all  Ar  = 1,  2,  • • • , K because  the  largest  element  in  a set  is  always  larger  than 
the  average  of  the  set.  This  indicates  that  all  the  users  are  overfaded  when  K < N. 

4.7.1.2  Case  2:  K > N 


Define  fj0  = 00.  No  matter  how  the  effective  noise  densities  of  all  users  are  dis- 
tributed, there  always  exists  a non-empty  set 

T = {k  : 0 < — T < w,  0 < k < N{1  + 7)  - K 7}. 

N(  1 + 7)  — K 7 — k 

Let 

/c*  = maxjfc  : k G J7}.  (4.51) 


If  k*  + 1 < iV(l  + 7)  — A'7,  from  the  definition  of  £:*,  we  know  that 

< 2^t_fc»+2  ^ 

7fc*+i  _ ^ + _ Kl  _ + ’ 


(4.52) 


This  is  equivalent  to  r/*.+i  < • If  k*  + 1 > N(1  + 7)  - A'7,  we  have 

1 < Ei=**+i  ^ because  0 < iV(l  + 7)  - K'y  - k*  < 1.  This  proves 

the  existence  of  k * in  Proposition  4.2.1.  An  equivalent  expression  of  < fifc* 

is  iv(i+7)ilc7-(fc»-i)  < %*•  Wbb  tbe  decreasing  order  of  effective  noise  densities,  we 
always  have  the  inequality 


v ' K ~ 

rll 

iV(l  + 7)  — K'y  — ( k * — 1) 


< -1, 
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which  means  that  k*  — 1 £ T . Inductively,  we  know  that  {1,2,---,  A:*  — 2}  G T . Since 
Vk-+2  < 7fe*+i>  we  always  have  the  inequality 

Y^i=k*+ 2 Vi 


dk*  + 2 < %*+!  < 


This  is  equivalent  to 


%*+2  < 


iV(l  + 7)  — Kj  — (k*  + 1) 
S/=fc*+3  ^ 


7V(1  + 7)  — Kj  — (/c*  + 2)  ’ 
which  means  that  k*  + 2 tfi  T . Inductively,  we  know  that  {k*  4-  3,  k*  + 4,  • ■ • , A'}  ^ T . 

This  proves  the  uniqueness  of  k*  and  the  set  of  overfaded  users  T — {1,  2,  • • • , k*}. 

The  proof  above  validates  a simple  procedure  that  can  be  used  to  generate  the  set 
of  overfaded  users. 


step  1:  Assume  k = 0 at  the  beginning  of  the  procedure. 

y'K’  - 

step  2:  If  rjk+i  > , then  increase  k by  1,  that  is,  k + 1 =>-  k,  and  repeat  the 

step  2. 


step  3:  Else,  set  k*  = k,  and  terminate  the  procedure. 

4.7.2  Proof  of  Proposition  4.2.2 

It  is  well  known  that  the  sum  of  the  diagonal  elements  of  a square  matrix  is  equal  to 
the  sum  of  its  eigenvalues.  When  a matrix  is  symmetric,  the  precise  relationship  between 
the  diagonal  elements  and  the  eigenvalues  can  be  described  using  the  theory  of  majoriza- 
tion  [70]. 

Definition  4.7.1  For  any  x — (x\,  x2,  ■ ■ ■ , xn)  G lZn,  let 

X[l)  > X[2]  > ■ > X[n] 

denote  the  components  in  x in  descending  order,  called  the  order  statistics  of  x. 

Definition  4.7.2  For  any  x = (27,  x2,  ■ • • , xn),  y = (t/i,  2/2,  • • ■ , yn)  € T^n.  we  say  tl iat  x 
is  majorized  by  y ( or  y majorizes  x)  if 

Si=l  x[%]  — Si=l  V[i\i  k = 1,  2,  • • • , n — 1 

En  sr^n 

i= 1 X[i\  ~ 1 V{i\ 


(4.53) 
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Majorization  makes  precise  the  vague  notion  that  the  components  of  a vector  x are  “less 
spread  out”  or  “more  nearly  equal”  than  the  components  of  a vector  y by  the  statement  that 
x is  majorized  by  y.  A simple  example  of  majorization  is  that  (A,  A,  • • • , A)  is  majorized 
by  (1,  0,  • ■ • , 0). 

Lemma  4.7.3  ([70],p  218)  If  H is  an  n x n Hermitian  matrix  with  diagonal  elements 
(h\,  h2,  • ■ • , hn)  and  characteristic  roots  (Ai,  A2,  • • ■ , An),  then  (Ai,  A2,  • • • , An)  majorizes 
(h\ , h2,  j hn) 


Lemma  4.7.4  ([70], p 220)  If  (hx,h2,  ■■■ ,hn ) and  (Ai,  A2,  • • • , Xn)  are  2 n numbers  such 
that  fAi,  A2,  • • ■ , An)  majorizes  (hi,  h2,  ■ • • , hn),  then  there  exists  a real  symmetric  matrix 
H with  diagonal  elements  (hi,  h2,  • • • , hn)  and  characteristic  roots  (Ai,  A2,  ■ • • , An). 


Tet  \ = 52£=k*+iPk 
a — N_k, 


7Et=fc»+i  m 


When  the  constraint  in  (4.9)  is  satisfied,  it  is 


straightforward  that 


Pk'  + l 


e(y) 


7 2-^i=k*+i  hi 


N(  1 + 7)  — Ky  — k " 
= e(7)(A  + r)fc*+i) 

7 Vk-+i  ~ A 


+ Vk*+ 1 


= A + 
< A. 


1+7 


(4.54) 


With  the  assumption  of  a decreasing  order  of  the  effective  noise  densities,  that  is,  fji  > 
7)2  > ■ • • > f}K,  and  the  classification  of  users  in  the  Proposition  4.2.1,  the  power  allocation 
given  in  (4.10)  satisfies  px  > p2  > ■ ■ ■ > Pk ■ Thus, 

(A,  — , A,  0,  — , 0)  majorizes  (p**+i,Pfc*+2,  • ■ • ,Pk)- 

can  be  justified  by  showing  that 

k 

J2  Pi  < (k-k*) A,  k = k*  + l,k*  + 2,  - ■ • ,N, 

i=k*+l 
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k 


Y Pr  < ( N-k*)X , k = N + 1,N + 2, ■■■, K — 1, 


i=k*  + 1 
K 


J2  K = (N-k‘)\. 


(4.55) 


i=fc*+l 


From  lemmas  4.7.3  and  4.7.4,  there  exists  a (K  — k*)  x ( K — k *)  hermitian  matrix  B 
with  diagonal  elements  (pfc*+i,Pfc*+2,  • • • ,Pk)  and  eigenvalues  (A,  • • • , A,  0,  • • • , 0).  Take 
eigenvalue  decomposition  of  B as  B = UAU77.  Because  B is  hermitian,  matrix  U is 
unitary.  Let  ui,  u2,  • • • , UAr_fc.  be  the  orthonormal  eigenvectors  of  B corresponding  to 
(. N-k *)  nonzero  eigenvalues  A.  Then  the  ( N-k *)  x(K  — k*)  spreading  sequence  matrix 
A can  be  constructed  as 


4.7.3  Proof  of  Proposition  4.2.3 

We  obtain  the  lower  bound  of  the  ETSC  for  two  cases  when  K < N and  K > N . 

4.7.3. 1 Case  1:  K < N 

Given  the  transmission  powers  pk,  for  k — 1, 2,  ■ • • , K,  as  defined  in  (4.8),  the 
ETSC  is  lower  bounded  by 


(4.56) 


It  is  straightforward  that  AA1  — 


N^l+^-Ky—k 


7 Ej=fc»+i  Vl 


|^I.andA«Ar  = B. 


K K 


ETSC 


k= 1 1=1 
K 


K 


k= 1 


k=  1 


K 


k=l 

K 


(4.57) 


Employment  of  orthogonal  sequences  achieves  this  lower  bound  when  K < N . 
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4.7.3.2  Case  2:  K > N 


Given  the  fixed  sequence  power  allocation  pk,  for  1 < k < K,  defined  in  (4.10), 
we  can  construct  an  N x K matrix  A such  that  A ^ A is  a A x K square  matrix  with 
diagonal  elements  {pi,P2>  • • • , Pk}-  The  ETSC  can  be  expressed  in  terms  of  the  nonzero 
eigenvalues  of  AH  A as 


ETSC  = tr(AAH  AAh) 

= tr(UAUHUAUH) 

N 

= EA-- 


(4.58) 


71  = 1 


Let  A = (Ai,  A2,  • • • , Aat,  0,  0,  • • • , 0)  and  p = (pi,p2,  ■ • • ,Pk)-  From  theory  of  majoriza- 

K-N  zeros 

tion  [70],  we  know  that  the  precise  relationship  between  the  eigenvalue  vector  A and  the 
diagonal  elements  vector  p is  that  A majorizes  p.  Thus  ETSC  minimization  problem  can 
be  written  as 

N 

min^A*  (4.59) 

71=1 

subject  to 

A € TV, 


where  M = {A  = (A1}  A2,  • • • , A#,  0,  • • • , 0)  € U*  : A majorizes  p}.  With  the  spreading 
sequences  Ar  constructed  in  (4.15),  Af  Ar  has  diagonal  elements  vector  p,  and  eigen- 
value vector 


A*  = 


= (A;,A*,-.-,A^,o,o,---,q) 


K-N  zeros 


7E(=fc*+i  Vi 


= (7771,  •••,7%*, 


7E/=fc*+i  Vi 


N(l  + 7)  — K'y  — k*  ’ ’ N(1  + 7)  — Kj  — k* 


o1^1o  ) 

K-N  zeros 


= {Pw-Pk-, 


E;=fc»+ 1 Pi 
N-k*  ’ 


E;=fc»+ 1 Pi 
N-k*  ’ 


(4.60) 


K-N  zeros 

It  is  straightforward  that  X*  e Af  and  A]  > AT]  > • • • > A *N  from  the  classification  of 
overfaded  users  and  nonoverfaded  users  and  the  proof  of  the  Proposition  4.2.2.  Notice  that 
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\2n  is  a convex  function  of  An,  and  hence  the  symmetric  convex  map  /(A)  = J2n=i  A n 
is  Schur-convex  [70],  A useful  property  of  a Schur-convex  function  is  that  if  Ai  6 J\ f 
majorizes  A2  G A f,  then  f{A\)  > /( A2).  If  we  can  prove  that 

A majorizes  A*,  VA  £ AT,  (4.61) 


the  eigenvalue  vector  A*  and  the  sequence  allocation  scheme  in  (4.15)  gives  the  minimum 
ETSC. 

Now  we  proceed  to  justify  (4.61)  using  an  approach  similar  to  the  one  given  by 
Viswanath  et  al.  [49],  Let  Ptot  = Y,k=\  Ah-  Without  loss  of  generality,  we  assume  that 
Ai  > A2  > • • • > Ajv  in  A.  First  we  notice  that 

k k k 


E^E^E**'  * = 1.2 


i=l  i— 1 i=l 

We  complete  the  remaining  of  the  proof  by  hypothesis.  Suppose 

k k 


for  some  k*  < k < N.  Since 


and 


Ea‘^Ea.' 

i= 1 i= 1 


N k 

y ] \ — Ptot  y ' A}., 

i=k+ 1 


Afc+i  > Afc+2  > • • • > Aat, 


(4.62) 


(4.63) 


we  have 


Hence, 


Afe+i  > 


Ptot  — Yli= i Aj 

N-k 


fc+i 


EA<  ^ 


i— 1 


Ptot  — X)j=i  Aj 

N-k 


i= 1 


> -A i.  + E" 


N — k N — k 


i— 1 


(4.64) 
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D AT  L 1 


k 


i— 1 


(4.65) 


This  is  true  for  all  k = k*,  k*  + 2,  • • • , N — 1.  Hence,  the  proof  is  complete. 

4.7.4  Proof  of  Proposition  4.3.1 

The  user  capacity  constraint  can  be  proved  using  a derivation  similar  to  the  one 
provided  in  Ulukus  [53], 

Lemma  4.7.5  ([71])  A necessary  and  sufficient  condition  for  a solution  x (x  > 0,  f 0)  to 
the  equations 


to  exist  for  any  c > 0,  f 0 is  that  s > r,  where  r is  the  largest  (also  called  Perron- 
Frobenius)  eigenvalue  of  matrix  T.  In  this  case,  there  is  only  one  solution  x,  which  is 
strictly  positive  and  given  by 


From  Lemma  4.7.5,  we  know  that  a feasible  solution  for  (4.26)  exists  if  and  only  if 


where  pc  is  the  Perron-Frobenius  eigenvalue  of  matrix  C.  This  means  that  given  a fixed 
set  of  signature  sequences  s*,  for  k = 1,  2,  • • • , K,  the  range  of  the  uniform  target  SIR 
achieved  by  all  K users  is  determined  by  the  Perron-Frobenius  eigenvalue  of  the  squared 
cross  correlation  matrix  C,  which  depends  only  on  the  signature  sequences  of  the  users. 
It  is  hard  to  relate  the  Perron-Frobenius  eigenvalue  pc  to  the  signature  sequences  clearly, 
however,  pc  is  shown  [53]  to  be  bounded  by 


(si  - T)x  = c 


X = 


(si  — T)-1c 


1 


(4.66) 


(4.67) 


77 


We  can  see  that  pc  is  lower  bounded  and  upper  bounded  by  a factor  of  the  TSC.  Since  we 
want  to  minimize  pc,  it  is  reasonable  to  minimize  the  lower  bound  over  the  space  of 
all  the  possible  signature  sequences.  To  do  so,  we  need  two  lemmas  presented  below. 

Lemma  4.7.6  ([48])  (Welch ’s  lower  bound  on  total  squared  correlation) 


Lemma  4.7.7  ([72])  (The  Unifonnly-Good  Property)  If  the  sequences  sk,fork  = 1,  2,  • • • , K 
are  such  that  equality  holds  in  (4.68)  then 


We  are  trying  to  minimize  pc,  then  the  best  we  can  do  is  to  choose  signature  se- 
quences so  that  (4.67)  achieves  the  minimum  of  the  lower  bound  with  equality.  Lemmas 
4.7.6  and  4.7.7,  combined  with  (4.67),  indicate  that  the  set  of  signature  sequences  that  min- 
imizes the  TSC  are  also  the  signature  sequences  that  minimize  pc  and  the  smallest  possible 
Perron-Frobenius  eigenvalue  is  pc  — jj-  From  (4.66),  we  know  that  the  maximum  number 
of  users  that  can  be  supported  in  the  system  is 


4.7.5  Proof  of  Proposition  4.3.2 

Summing  up  all  constraints  in  (4.24)  and  making  use  of  the  uniformly-good  prop- 
erty of  the  WBE  sequences  (Lemma  4.7.7)  that  the  sum  of  each  column  of  C is  we  can 
obtain  a lower  bound  on  Pwbe'- 


(4.68) 


A;=l  i=l 


(4.69) 


(4.70) 


e(y)PWBE^  + e(7)X^fc 


(4.71) 
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This  implies 


We (7)  Ef=i  Vk 


WBE 


> 


(4.72) 


N-Ke{i)  ‘ 

Lemma  4.7.8  (Lemma  3.3  in  Zhang  and  Chong  [56])  Suppose  A and  B are  N x N sym- 
metric matrices,  and  A > B.  Then  for  any  x,zg  (x7  Az)2  > (x7  Bz)2. 


From  the  proof  of  user  capacity,  we  note  that  the  largest  eigenvalue  of  C is 
Employing  Lemma  4.7.8,  we  can  obtain  an  upper  bound  on  Pwbe 

Pw BE ~ l^P  = 

< 


where  1 is  a K x 1 column  vector  with  all  elements  as  1.  Combining  (4.72)  and  (4.73),  we 
obtain  (4.29). 

4.7.6  Proof  of  Convergence  of  Algorithm  4.5.1 


1 H 


1 H 


1 


n -1 


e(7)' 
1 . 


C 

K 


_e(y)  N 

We (7)  Ef=  1 m 
N — Ke(  7) 


-1 


'H 


(4.73) 


Algorithm  4.5.1  is  a generalization  of  the  TSC  reduction  algorithm  by  Ulukus  [51] 
to  the  system  with  unequal  transmission  powers  and  effective  noise  densities.  Here,  con- 
vergence of  Algorithm  4.5.1  is  proved  using  the  idea  in  Ulukus  [51]. 

First  we  prove  that  the  ETSC  decreases  monotonically  at  each  iteration  using  Al- 
gorithm 4.5.1.  K sequences  are  updated  sequentially  starting  with  the  first  sequence.  For 
the  kth  sequence,  we  note  that 


AtJj]A  fji]  = a*[j]a  «\j]  + A„[j]A  f [j],  (4.74) 

where  At(j]  = (a,  [7  + 1],  • ■ • , a t-i\j  + 1],  at+1[j],  ■ ■ • , aK[j]|.  The  resulting  change  in  the 
ETSC  when  updating  the  kth  spreading  sequence  using  Algorithm  4.5.1  is 


A*  [7]  — tr[(afc[j  + l]a^ \j  + 1]  + Ak[j]Ak  [j])2] 
tr[(a*[7]afb1  + Ak[j]A%  [j])2] 
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= 2a f [j  + 1]R*  [j] ak  [7  + 1]  - 2af  [j]Rfc [j]afc [j] 


= 2pfc- 


- 2afc  \j}TLk[j]ak[j]. 


(4.75) 


af[i]Rfc2[7]afe[i] 

Since  Rfc[)]  is  of  full  rank  and  ak[j]ak[j]  = pk,  we  can  employ  Cauchy-Schwartz  inequal 


ity  to  write 

pi  = ^m^ml[j]ak[j})2 
< l|af[i]R7[j]iri|af[y]Rl[j]||2 

= (af[7]Rfe-1[7]a,[i])(af[i]Rfe[7']afe[i]))  (4.76) 


and 


(affijR^bla^i])2  < ||R,-1[7]afe[7]||2||afe[y]||2 

= (af[j]Rfc2[j]afe[y])pfc.  (4.77) 


Combining  the  two  inequalities  above  and  (4.75)  yields  A k[j]  < 0,  for  all  j,  k = 1, 2,  • • • , K. 
Moreover,  A k[j]  = 0 for  f < k < K if  and  only  if 

a k[j]  = Pk^ljWlj],  (4-78) 


for  1 < k < K.  An  equivalent  expression  of  (4.78)  can  be  obtained  using  the  matrix 
inversion  formula: 

a*  {]}  = 5/tR^1  b']a*  [j]  - (4.79) 

Here,  f3k  and  gk  are  nonzero  scalars  the  values  of  which  are  decided  by  the  system  and 
channel  parameters.  This  indicates  that  the  ETSC  keeps  decreasing  until  each  spreading 
sequence  is  the  same  as  the  MMSE  receiver  weight  vector  up  to  a scaling  factor.  Since  the 
ETSC  decreases  strictly  at  each  iteration  except  at  the  fixed  point  described  above  and  the 
ETSC  is  lower-bounded  given  fixed  sequence  powers  allocation,  convergence  of  the  ETSC 
is  guaranteed. 

At  a fixed  point,  the  spreading  sequences  do  not  change.  Now,  we  analyze  the 
fixed-point  property  of  Algorithm  4.5.1.  Let  us  denote  a fixed  point  of  Algorithm  4.5.1  by 
the  N x K sequence  matrix  AT  — [ai,  • • ■ , a*,  • • • , a^r]. 
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For  the  case  of  K < N,  inserting  = ATA % + fjkI  into  (4.79),  we  can  rewrite 

it  as 

ATAHvak  = (gk  - rjk)ak  = k = l,2,---,K.  (4.80) 

Pre-multiplying  (4.80)  with  af , for  1 < l < K,  and  writing  the  resulting  equations  in 
matrix  form,  we  get 

A%AtA”At  = A“A  r$,  (4.81) 

where  $ = diagtyi,  02,  ■ • • , Referring  to  Corollary  1 in  Ulukus  [51],  if  Algo- 

rithm 4.5.1  starts  with  a full  rank  matrix  Ax[0],  A^A r is  invertible.  Multiplying  both 
sides  of  (4.81)  with  (A^A r)_1  yields 

A"Ar  = $.  (4.82) 

Thus  A^At  is  diagonal,  which  indicates  orthogonality  of  the  spreading  sequences  at  the 
fixed  point.  The  sequence  powers  are  fixed  during  the  adaptation,  thus  we  have  diag[ Af/,AT\ 
bi.P2,  • • • ,Pk]T  = 

For  the  case  of  K > N,  take  eigenvalue  decomposition  on  ATA%  = UAUH 
and  let  ak  = UHak  = [aM,  ak>2,  ■ ■ • , aktN]T.  We  note  that  the  elements  ak>n,  for  n = 
1,  2,  • • ■ , N in  the  sequence  vector  are  the  projections  of  the  spreading  sequence  ak  onto 
the  channels  defined  by  the  eigenvalue  decomposition.  Based  on  (4.80),  we  get 

Aak  = (f)kak,  k = 1,2,-**,  K.  (4.83) 

Eq.  (4.83)  indicates  that  (j)k,  for  1 < k < K , are  eigenvalues  of  ATA%.  Pre-multiplying 
(4.83)  for  User  k with  a/7,  and  pre-multiplying  (4.83)  for  User  l with  a j? , we  obtain  these 
two  equations 

af1  Aak  - (j)ka”ak, 
af  Aa,  = <f>ia.kai. 

3 We  use  the  notation  diag[4>  \ , 02,  • ■ • , 4>k]  to  indicate  ait  x K diagonal  matrix  with  <j>i , 02, ' ' ' i <Pk  on 
the  main  diagonal. 
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Since  A is  diagonal,  we  must  have 

bktf&k  = Mfa*,  (4.84) 

for  all  /,  k — 1,  2,  • • • , K,  which  indicates  that  (j)k  = </>/  when  af  a*  ^ 0.  Eq.  (4.84)  has 
several  implications: 

(i)  Define  the  spreading  sequence  set  A = {ai , a2,  • ■ • , a k}  at  the  fixed  point.  If  A does 

not  split  into  two  or  more  orthogonal  subsets,  then  <j)k  = cj)  for  all  k = 1, 2,  • • • , K, 
which  means  that  AtAj.  is  an  identity  matrix  with  equal  eigenvalues  Ai  = ■ • • = 
\N  = = 4 ).  If  A splits  into  several  orthogonal  subsets,  <£*  within  each  subset 

are  equal. 

(ii)  There  may  be  different  fixed  points  with  different  ways  of  space  partitioning,  that  is, 

the  ETSC  may  converge  to  local  minimum.  The  space  partitioning  scheme  described 
in  Section  3.2  gives  the  global  minimum  ETSC. 

Having  established  the  monotonicity  and  convergence  of  Algorithm  4.5.1,  we  can 
prove  that  Algorithm  4.5.1  can  always  converge  to  a fixed  point  although  it  may  not  be 
the  fixed  point  giving  the  the  global  minimum  ETSC.  Simulation  results  show  that  Algo- 
rithm 4.5.1  always  converges  to  the  fixed  point  giving  the  global  minimum  ETSC  starting 
with  random  sequences. 

4.8  Summary 

Two  power  and  sequence  allocation  schemes  are  proposed  to  support  users  with  a 
uniform  target  SIR  for  a single-cell  forward-link  CDMA  system.  The  first  scheme,  referred 
to  as  the  minimum-ETSC  scheme,  is  proposed  for  use  of  MMSE  signal  reception.  With 
this  scheme,  orthogonal  sequences  are  allocated  when  K < N.  When  K > N,  users  in 
the  system  are  classified  into  two  classes,  namely  overfaded  users  and  nonoverfaded  users, 
similar  to  the  idea  of  user  classification  in  Viswanath  et  al.  [50].  Overfaded  users  are  al- 
located orthogonal  “channels”  and  nonoverfaded  users  share  the  remaining  “channels”.  In 
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both  cases,  the  sequences  allocated  are  the  class  of  sequences  that  minimizes  the  ETSC. 
The  second  scheme,  called  minimum-TSC  scheme,  is  proposed  for  use  of  matched-filter 
signal  reception.  With  this  scheme,  the  unit-energy  signature  sequences  of  users  are  the 
class  of  sequences  that  minimize  the  TSC.  Power  efficiencies  of  the  proposed  schemes 
are  examined  through  both  theoretical  analysis  and  numerical  method,  which  indicate  that 
the  minimum-ETSC  scheme  has  better  power  efficiency  than  the  minimum-TSC  scheme. 
In  addition  to  its  decent  power  efficiency,  another  important  advantage,  which  makes  the 
minimum-ETSC  scheme  desirable  from  an  application  viewpoint,  is  that  the  set  of  se- 
quences can  be  easily  constructed  iteratively  using  the  proposed  adaptive  algorithms.  A 
joint  transmitter-receiver  adaptive  structure  combined  with  channel  estimation  is  also  pre- 
sented to  facilitate  the  use  of  the  minimum-ETSC  scheme  in  real-life  CDMA  forward  links. 


CHAPTER  5 

SEQUENCE  OPTIMIZATION  OVER  MULTIPATH  FADING  CHANNEL 

Spreading  sequence  optimization  in  CDMA  systems  with  multipath  is  attempted  by 
Rajappan  and  Honig  [57]  and  Ulukus  [58]  using  numerical  adaptation.  However,  no  theo- 
retical optimal  solution  is  established  in  the  current  literature.  In  this  chapter,  we  consider 
the  problem  of  sequence  optimization  in  terms  of  supporting  all  users  in  the  system  with  a 
minimum  amount  of  total  transmission  power  in  a DS-CDMA  point-to-point  transmission 
over  multipath.  In  Section  5.1,  we  define  the  system  model  and  give  a brief  discussion  of 
optimal  MMSE  signal  reception.  In  Section  5.2,  we  discuss  user  capacity  and  sequence 
optimization  problem  of  a CDMA  point-to-point  transmission  with  multipath  in  detail.  In 
Section  5.3,  a suboptimal  approach,  which  uses  eigenvectors  of  path  matrix  as  orthogonal 
spreading  sequences,  is  presented.  Comparison  of  this  suboptimal  approach  to  the  opti- 
mal solution  is  made  via  computer  simulation.  Proofs  of  the  propositions  are  provided  in 
Section  5.4.  Finally,  conclusions  are  drawn  in  Section  5.5. 

5.1  System  Model 

We  describe  the  system  model  of  a point-to-point  system  over  multipath  in  this 
section.  In  point-to-point  transmission,  receivers  of  all  users  are  co-located.  Hence,  re- 
ceived signals  at  all  receivers  undergo  the  same  multipath  channel  response,  and  the  power 
spectral  densities  (PSD)  of  the  thermal  noise  at  all  receivers  are  the  same,  that  is,  that  is, 
Hfc  = H and  r]k  = r)  for  A:  = 1,  2,  • • • , K.  With  these  assumptions,  the  general  signal 
model  described  in  Chapter  2 can  be  simplified.  The  (N  + L — 1)  samples  at  the  output  of 
the  chip  matched  filter  can  be  arranged  into  a column  vector 

K 

r = H ^ &om)am  + n,  (5.1) 

771  — 1 
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Accordingly,  the  optimal  MMSE  receiver  weights  of  the  A;th  receiver  are  given  by 


wfc  = R/Hafc, 


(5.2) 


where  Rr  is  the  total  correlation  matrix  observed  at  the  receivers,  which  is  given  by 

Rr  = E[rrH]  = HATAf  Hh  + rjl.  (5.3) 

In  above,  AT  = [ai  a2  • • -a^j  is  an  N x K matrix  formed  by  grouping  all  K spreading 
sequence  vectors.  We  note  that  the  total  correlation  matrices  observed  by  all  the  users  are 
the  same  according  to  the  system  model  we  adopted.  With  the  optimal  MMSE  weights, 
the  receiver  output  SIR  is  given  by 

SIRfc  = af  H^R-'Ha,,  (5.4) 


for  A;  = 1,  2,  • • • , K,  where  Rfc  is  the  noise-plus-interference  correlation  matrix,  which  is 
defined  as 

Rfc  = Rr  - Ha*af  Hff.  (5.5) 


An  equivalent  expression  of  the  SIR  using  the  matrix  inversion  formula  can  be  derived  as 
well: 


SIRfc 
1 + SIRfc 


afH^R^Hafc. 


5.2  Sequence  Optimization 


(5.6) 


5.2.1  User  Capacity 


The  sequence  optimization  problem  is  characterized  by  choosing  a set  of  spreading 
sequences  so  that  the  SIR  targets  are  met  for  all  data  streams  in  the  system  with  a minimum 
amount  of  total  transmission  power.  For  simplicity,  we  assume  a uniform  SIR  target  7 for 
all  users.  Mathematically,  the  optimization  problem  can  be  expressed  as 


K 

min  llafc 

fc=i 


subject  to 


a''H"Rt-'Ha*  > 7 


(5.7) 
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for  1 < k < K.  Equivalently,  the  constraints  can  be  written  as 


af  H^R^Ha*  > e(7). 


(5.8) 


Here,  we  follow  the  notation  in  Tse  and  Hanly  [35]  to  denote  e( 7)  = The  quantity 
e(7)  represents  the  effective  bandwidth  [35]  that  a user  takes  when  achieving  the  MMSE 
receiver  output  SIR  target  7.  The  constraints  in  (5.8)  can  be  written  in  matrix  form  as 


Due  to  the  limited  system  capacity,  there  may  not  be  a feasible  allocation  if  there 
are  too  many  users  in  the  system  or  if  the  SIR  target  of  the  users  is  too  high.  We  define  the 
system  capacity  as  the  maximum  number  of  users  the  system  can  support  with  a certain 
target  SIR.  Proposition  5.2.1  gives  a necessary  and  sufficient  condition  for  the  number  of 
users  that  can  be  supported. 

Proposition  5.2,1  The  necessary  and  sufficient  condition  for  the  existence  of  a feasible 
solution  of  this  optimization  problem  is  that  the  processing  gain  N,  the  number  of  users 
K,  and  the  SIR  target  7 satisfy  the  constraint: 


There  is  a similar  interpretation  for  this  constraint  as  in  the  case  of  the  AWGN 
channel  and  the  forward-link  transmission  over  slow  flat  fading.  All  K users  are  admissible 
if  and  only  if  the  sum  of  their  effective  bandwidth  is  less  than  the  spreading  gain  of  the 
system. 

5.2.2  Optimal  Spreading  Sequences 

A general  solution  for  sequence  optimization  with  colored  additive  noise  is  pre- 
sented in  Wong  and  Lok  [46].  The  optimization  problem  (5.9)  can  be  solved  using  a 
similar  procedure. 


diag  [A^H^R^HAt]  > e(7)I. 


(5.9) 


I<e{ 7)  < N. 


(5.10) 
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Let  H = U£Vff  be  the  singular  decomposition  of  the  path  matrix  H with  the 
singular  values  sn  arranged  in  a descending  order  in  the  diagonal  matrix 


Sl 


S = 


s2 


(5.11) 


Sn 


L J (N+L-l)xN 

In  the  following  part,  we  assume  that  H is  of  full  rank,  that  is,  sn  > 0 for  n — 1,2,**-,  iV. 
The  results  can  be  easily  extended  to  the  case  when  H is  not  a full-rank  matrix. 


Proposition  5.2.2  Suppose  that  Ke( 7)  < N,  let  Pn,  for  n = 1, . . . , min  (A,  N),  be  the 
solution  of  the  optimization  problem 


mm 


rain  {K,N) 

EPn 

ET 

71  = 1 U 


subject  to 


min  (K,N) 

E 


n= 1 


Pn 

1 + Pn 


Ke{  7) 


Pi  > P2  > ‘ ' ' > Pmin(K,N)  > 0. 


(5.12) 


For  the  case  of  K < N,  construct  a K x K unitary  matrix  W such  that  the  diagonal 
elements  of 


W diag 


Pi  P2  Pk 

1 + Pi  ’ 1 + P2 ' 1 + Pk. 


are  all  equal  to  e(y).  Then, 


A T = y/rjV  diag 


y/Pi  VP2 

Sl  S2 


sk 


(5.13) 


is  a solution  to  the  optimization  problem  in  (5.9)  when  K < N,  where  V is  an  N x K 
matrix  taken  from  the  first  K columns  (keeping  the  original  order ) of\.  For  the  case  of 
K > N,  construct  a K x K unitary  matrix  W such  that  diagonal  elements  of 


W diag 


Pi  P2  Pn 

1 + Pi  1+^2  1 + Pn 


WH 
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are  all  equal  to  e(7).  Then, 


A T = \/W  diag 


Vfil 

si 


VP 2 

S 2 


0 


VPn 

SN 


w 


(5.14) 


is  a solution  to  the  optimization  problem  in  (5.9)  when  K > N.  The  resulting  minimum 
total  transmission  power  for  both  cases  is 

K min  (N,K) 


P-min  — ^ ' 1 1 3-fe  1 1 V '/’  ' 


% (5.15) 

k= 1 n=l 

We  note  that  the  matrix  W described  above  can  be  obtained  by  an  iterative  pro- 
cedure described  in  Wong  and  Lok  [46].  Solution  for  the  scalar  optimization  problem  in 
(5.12)  is  given  in  the  Proposition  5.2.3. 


Proposition  5.2.3  Let  K*  denote  an  integer  such  that 

0 < K*  - Ke{j)  < 1, 

then  there  exists  a unique  integer  K*  in  {K* , K*  + 1,  . . . , min(Ar,  K)}  such  that 

Yj->K,-  Ke{i) 


(5.16) 


for  n = 1, . . . , A*,  and 


^ n - Ke^) 


(5.17) 


(5.18) 


for  n = K*  + 1,  • • • , min(Ar,  K).  Then, 

{V'^*  in 
2-k=i  sk 

K,-Ke{  7) 

0 if  K\  + 1 < n < min  (A,  K) 

is  a solution  to  the  optimization  problem  in  (5.12). 


Km  Sjl 

- 1 if  1 < n < 


(5.19) 


If  we  denote  the  eigenvectors  of  the  matrix  G = HffH  as  the  channel  space 
spanned  by  the  DS-CDMA  system  with  path  matrix  H,  construction  of  the  optimal  se- 
quences A t stated  above  has  the  flexibility  and  tendency  of  choosing  the  good  channels 
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with  large  eigenvalues,  while  discarding  or  allocating  less  users  into  the  bad  channels  with 
small  eigenvalues.  Solution  of  the  scalar  optimization  problem  (5.12),  which  is  provided 
in  Propostion  5.2.3,  gives  the  optimal  strategy  when  choosing  channels.  In  the  absence  of 
multipath,  such  as  the  channel  model  of  DS-CDMA  over  AWGN,  the  channel  matrix  is 
an  identify  matrix,  that  is,  H = I and  G = I,  which  means  that  the  channel  space  are 
“equally  good”.  In  this  case,  it  is  intuitively  true  that  orthogonal  sequences  are  allocated  as 
spreading  sequences  when  K < N,  and  the  whole  channel  space  is  evenly  occupied  when 
K > N,  which  is  consistent  with  the  results  by  Wong  and  Lok  [46,  52]  and  Viswanath 
[50], 


Neither  spreading  sequences  AT  nor  the  received  sequences  HA7’  are  orthogonal, 
in  general,  with  the  optimal  scheme  presented  in  Section  5.2.  Multipath  convolution  de- 
stroys the  orthogonality  among  the  received  sequences  eventhough  transmitted  sequences 
are  orthogonal,  however,  with  one  exception  that  we  will  address  in  this  section. 

5.3.1  Eigenvectors  as  Signature  Sequences 

Take  eigenvalue  decomposition  of  the  matrix  G = HffH  as  G = VAV^,  where 
A = £H£  = diag  [sf,  s\, . . . , s2N]  with  eigenvalues  s £ arranged  in  a descending  order  at 
the  diagonal  of  A.  Let  vn  be  the  eigenvector  corresponding  to  the  eigenvalue  s£,  and  let 
r„  = Hv„,  for  n = 1,  2,  • • • , N,  we  have 


Hence,  orthogonalities  among  the  received  sequences  are  maintained  when  unit-energy 
eigenvectors  of  G are  used  as  signature  sequences.  Spreading  sequence  of  the  kth  user 
is  assigned  as 


where  is  the  transmission  power  of  the  fcth  user,  for  k = 1,  2,  • • • , K.  Only  K < N users 


5.3  A Suboptimal  Scheme  with  Orthogonal  Sequences 


for  n = m, 
for  n ^ m. 


(5.20) 


cifc  \/Pfevfc) 


(5.21) 


are  admissible  in  the  system  using  this  scheme  due  to  the  limited  number  of  eigenvectors. 
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With  orthogonal  sequences,  MMSE  receiver  weights  are  the  same  as  those  of  the  matched- 
filter  receiver.  Receiver  output  SIR  of  the  kth  user  with  the  matched- filter  receiver  can  be 
derived  as 


SIR*  = 


Pk  |vfH"Hvfcl2 

Pks'j 

7 

V 


(5.22) 


for  k = 1, 2, . . . , K.  Hence,  transmission  power  of  the  kth  user  should  be  allocated  as 


„ _ rn 

Pk  2 

4 


(5.23) 


to  satisfy  the  target  SIR  7.  From  (5.23),  we  note  that  the  eigenvectors  with  larger  eigenval- 
ues should  be  assigned  to  users  to  lessen  the  transmission  power.  With  a descending  order 
of  the  eigenvalues  s) \,  we  should  always  choose  the  left  K columns  of  V as  the  signature 
sequences  for  all  K users.  The  total  transmission  power  required  by  this  scheme  is 

K 


n K - 1 


k= 1 


(5.24) 


5.3.2  Power  Efficiency  Evaluation 


Power  efficiency  of  this  scheme  is  evaluated  using  computer  simulations.  We  as- 
sume a sample  system  with  spreading  gain  N = 64  and  the  number  of  resolvable  paths 
L = 4.  There  are  K = 64  users  in  the  system.  Channel  gain  for  each  path  within  h is  a 
Raleigh  distributed  random  variable.  We  conduct  100  independent  realizations  and  com- 
pare the  total  transmission  power  PT  in  (5.24)  required  by  this  suboptimal  scheme  to  the 
minimum  power  Pmin  in  (5.15)  required  by  the  optimal  scheme.  For  ease  of  comparison, 
we  normalize  PT  with  Pmin . Simulation  results  are  shown  in  Fig.  5.1.  We  can  see  that 
the  eigenvector  allocation  scheme  consumes  more  power  for  all  the  realizations,  which  is 
consistent  with  the  theoretical  analysis. 
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Figure  5.1:  Power  consumption  of  the  schemes:  N=64,  K=64,  7=4  dB 

5.4  Proof  of  Propositions 
5.4.1  Proof  of  Proposition  5.2.1 

Note  that 


A^H^R^HAr 


= A*H"  [HArA"Hff  + r/l]  1 HA 

-i 


AfHH 


HtjjH 


HAr A"H 

L Vv  Vv 


hat 


= I 


' A^HffHAT 
1 + — 


A"(HffH) 

y/V 


(H//H)^ArAf(HffH)l 


-l 


(HffH)5AT 


(5.25) 


Let  B = we  can  rewrite  the  transmitter  optimization  problem  in  Section  5.2  in 

terms  of  the  matrix  B: 


mintr  [Bff(HffH)_1B] 
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subject  to 

diag  [I  - (I  + B^B)-1]  > e(7)I. 
Equivalently,  the  constraints  can  be  written  as 

diag  [Bff (I  + BBff)_1B]  > e(7)I. 


(5.26) 


(5.27) 


The  necessary  and  sufficient  condition  states  that  if  the  target  SIR  is  set  as  7,  then  at  most 
K = L^yJ  users  can  be  admitted  into  the  system. 

We  first  prove  the  necessity.  If  a feasible  solution  exists,  the  constraints  in  (5.27) 
have  to  be  satisfied.  In  particular,  the  sum  of  the  diagonal  elements  on  both  sides  have  to 
satisfy  the  corresponding  inequality: 


Keirf)  < tr  [BH(I  + BBff)_1B] 

- tr  [BBff(I  + BBh)_1] 

= tr  [UAUff(I  + UAU77)-1] 

N \ 

= E " 


71=1 

< N, 


l + Ki 


(5.28) 


where  we  have  performed  spectral  decomposition  for  BB77  into  UAUH.  Because  BBH 
is  Hermitian  and  semi-positive  definite,  all  the  eigenvalues  Xn  > 0,  for  1 < n < N,  and  U 
is  a unitary  matrix. 

Now,  we  prove  the  sufficiency.  We  can  construct  an  N x K Welch  bound  equality 
(WBE)  sequence  matrix  S = [si  s2  • • • sK]  composed  of  K unitary  norm  column  vectors. 
Mathematically,  the  WBE  sequence  matrix  is  characterized  by 


SSH  = % I 
diag  [S*S]  = I. 


(5.29) 


The  WBE  sequence  matrix  is  scaled  by  a constant  g to  give  the  desired  matrix  B = gS 
with  the  property  that 

diag  [B*(I  + BB^-’B]  = g2  I.  (5.30) 
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When  Ke{ 7)  < N,  the  equation  9\K  = e^)  has  a solution.  Hence  the  desired  matrix 
B can  be  constructed. 

5.4.2  Proof  of  Proposition  5.2.2 


With  the  singular  value  decomposition  H = U£Vff  in  (5.1 1),  we  have  the  spectral 
decomposition 


(HhH)  1 = V diag 


1 1 1 

3 ’ 3 ’ ' ‘ ' ’ 3" 
*1  *2 


J 


V 


H 


(5.31) 


of  the  matrix  (H^H)  \ with  eigenvalues  . . . , 4-  arranged  in  an  ascending  order. 

To  proceed  we  need  to  make  use  of  the  following  result,  the  proof  of  which  can  be  found, 
for  example,  in  Marshall  [70,  pp.  249], 


Lemma  5.4.1  Suppose  X and  Y are  two  Hermitian  N x N matrices.  Arrange  the  eigen- 
values Xi  of  X in  a descending  order  and  the  eigenvalues  yl  ofY  in  an  ascending  order. 


Then 


N 

tr[XY]>J2xiyi-  (532) 

i— 1 

Let  Pi,P 2,  ■ • ■ , Pmin (n,k)  be  the  eigenvalues  of  BffB  arranged  in  descending  order, 


and  apply  this  lemma  to  tr[BK(HwH)  1B],  we  have 


tr[B  (H  H)  B 


tr[(HffH)_1BBff] 

min(Af,K’) 


> 


X 

n— 1 


Pn 


(5.33) 


Now,  consider  the  optimization  problem  in  (5.12)  in  the  statement  of  Proposition  5.2.2. 
Based  on  (5.33),  we  see  that  the  optimization  problem  in  (5.12)  concerns  the  minimization 
of  the  lower  bound  of  tr[Bff(HffH)-1B]  under  a relaxed  constraint: 


tr[I  - (I  + B"B)-1]  = ]T  ffT  = Ke(7).  (5.34) 

-L  t Pn 

71=1 

Assuming  that  /3X,  - • - > Pmin (n,k)  is  the  solution  of  the  optimization  problem  in  (5.12), 

we  proceed  to  justify  the  validity  of  sequence  matrix  construction  scheme  in  (5.13)  and 


(5.14). 
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For  the  case  of  K < N,  choose  V and  W as  stated  in  Proposition  5.2.2,  the  choice 


of 


B = V diag  [ y/pi,  y/fo,  ■ ■ ■ , 


W 


H 


(5.35) 


gives  a solution  to  the  optimization  problem  in  (5.27)  when  K < N.  First,  notice  that  from 
direct  substitution,  tr[B^  (H^Ft)  1 B]  = ^n=i  Hence  this  choice  of  B attains  the 
lower  bound  in  (5.33).  Second,  for  this  choice  of  B,  the  lower  bound  of  tr[B  H (H^H)  1 B] 
is  minimized  since  the  eigenvalues  of  B^B  are  f}x,  /32, . . . , Pk,  which  is  the  solution  to  the 
minimization  of  the  lower  bound.  Combining  these  two  observations,  we  see  that  the 
choice  in  (5.35)  solves  the  optimization  problem  in  (5.27)  with  the  relaxed  constraint.  If 
we  can  show  that  this  choice  of  B also  satisfies  the  original  constraint,  it  will  also  be  a 
solution  to  the  optimization  problem  in  (5.27).  This  is  evident  from  the  construction  of  W 
as 


diag  [I  - (I  + B^B)"1] 

Pi  (3 2 


diag 

e(7)I. 


W diag 


(3k 


1 + Pi  ’ 1 + /?2  ’ ’ 1 + (3k  . 


W 


H 


(5.36) 


H 


For  the  case  of  K > N,  choose  V and  W as  stated  in  Proposition  5.2.2,  it  is  easy 
to  show  that  the  choice  of 

VA 

\Z/02 

B = V diag  . 0 W 

v/Av 

attains  the  minimum  of  the  lower  bound  of  tr[B/7  (H^H)  1 B],  and  the  SIR  constraints 
are  satisfied  from  the  construction  of  W : 


(5.37) 


diag  [I  - (I  + B^B)-1] 

fh  fo 


diag 

e(7)I 


W diag 


(3n 


1 + Pi  1 + /?2  1 + 0N 


> o, . . - , o 


w 


H 


(5.38) 
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Hence,  the  choice  of  B in  (5.37)  gives  a solution  to  the  optimization  problem  in  (5.27) 
when  K > N.  We  note  that  the  W described  above  can  be  obtained  by  an  iterative 
procedure  described  in  Wong  and  Lok  [46], 

Obviously,  with  Ar  = ^(H^HQ^B,  the  original  optimization  problem  in  (5.9) 
is  also  solved  as  indicated  in  (5.13)  and  (5.14). 

5.4.3  Proof  of  Proposition  5.2.3 


First  notice  that 


for  n = 1,  2, . . . , K*. 


Y,->K'~  Ke(i  )• 

' Sk 

i i ^ 


Now,  let  an  = for  n = 1, . . . , min  (TV,  K),  and  rewrite  (5.12)  as 


mm 


min(N,K ) 

E 


1 Qtr ] 


n=l  1 


subject  to 

min(N,K) 

^2  an  = Ke{~f) 

n=l 

1 ^ ^ Oi 2 ^ ' 2 Chnir^iVj/C)  0. 


(5.39) 


(5.40) 


Writing  /(a)  = J2n=iN'K^  ‘l ’s  easY  t0  see  ^at  /(a)  *s  convex  an^  continuously 

differentiable  in  the  region  in  TZmm(N’K)  defined  by  the  constraints  in  (5.40).  Based  on  this 
observation,  the  minimization  problem  in  (5.40)  can  be  solved  by  applying  Theorem  4.4.1 
in  Gallager  [73,  pp.  87],  which  essentially  states  that  a'  is  a solution  to  (5.9)  if  and  only  if 
there  is  a /r  such  that 

for  all  n that  a!n  > 0 (5.41) 

for  all  n that  a!n  = 0.  (5.42) 

We  claim  the  following  choices  of  a'  and  /j,: 

1 - if  1 < n < IU 

0 if  /F*  + 1 < n < min(Ar,  K ) 


df{a') 

dan 

d/(a') 

dan 


> M 


(5.43) 
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A 


1 


K, 

Er 


2 


(5.44) 


K*  ~ Keffl)  “ sfc 

It  is  easy  to  see  that  the  choice  in  (5.43)  satisfies  the  constraints  in  (5.40).  Moreover,  it  can 
be  easily  proved  that  the  two  conditions  (5.41)  and  (5.42)  are  indeed  satisfied.  To  see  this, 
notice  that 


df(a')  = 1 1 

dan  si  (1  - an)2 

( H if  1 < n < A'* 

| -y  >f  K*  + 1 < n < min(iV,  K). 

v sn 


(5.45) 


Hence  the  proof  will  be  completed  if  we  can  show  that  4-  > ^ for  K*  < n < min (N,  K ). 
It  suffices  to  show  that  \ > n or  equivalently  > K*  — Ke( 7),  for  K*  + 1 < 

n < K.  This  can  be  easily  proved  by  the  definition  of  A'*.  Hence  (5.43)  is  a solution  to  the 
minimization  problem  in  (5.40). 

Transforming  back  to  f3n,  the  second  statement  in  the  proposition  that  (5.19)  is  a 
solution  to  the  original  minimization  problem  in  (5.12)  is  proved. 


5.5  Summary 

In  this  chapter,  we  investigate  the  sequence  optimization  problem  in  CDMA  point- 
to-point  transmission  with  multipath.  Solution  of  sequence  optimization  with  additive  col- 
ored noise  for  the  case  when  the  number  of  desired  users  is  smaller  than  or  equal  to  the 
spreading  gain  is  suggested  in  Wong  and  Lok  [46],  Following  an  approach  similar  to  the 
one  adopted  in  Wong  and  Lok  [46]  as  well  as  extending  the  results  in  Wong  and  Lok  [46] 
to  the  case  when  the  number  of  users  is  larger  than  the  spreading  gain,  we  solve  the  opti- 
mization problem  of  supporting  users  with  a uniform  target  SIR  using  the  minimum  total 
transmission  power.  Comparisons  of  optimal  solution  to  a suboptimal  approach,  which 
uses  eigenvectors  of  the  path  matrix  directly  as  orthogonal  spreading  sequences,  are  made 
via  computer  simulations. 


CHAPTER  6 

A HYBRID  ARQ  PROTOCOL  OVER  OPTIMAL  SEQUENCE  CDMA  LINKS 


The  capacity  of  a CDMA  link  is  determined  by  the  maximum  transmission  power, 
the  spreading  gain,  and  the  noise  power  of  the  system.  Through  the  employment  of  the 
optimal  spreading  sequences  proposed  in  Wong  and  Lok  [46,  52]  and  Viswanath  et  al.  [50], 
transmission  power  can  be  utilized  more  efficiently.  In  this  chapter,  we  propose  a type-I 
hybrid  ARQ  protocol  that  can  integrate  the  usage  of  optimal  sequences  in  a CDMA-based 
wireless  link.  In  Sections  6.1  we  describe  the  type-I  hybrid  ARQ  protocol  at  the  data-link 
layer.  Under  the  proposed  scheme,  packets  in  the  transmission  buffer  are  divided  into  two 
classes,  namely  retransmission  packets  and  newly  arrived  packets.  Packets  in  each  class  are 
assigned  optimally  selected  spreading  sequences  and  a uniform  target  SIR.  In  Sections  6.2, 
we  provide  a brief  introduction  of  the  structure  of  the  optimal  sequence.  In  Section  6.3,  we 
construct  a resource  allocation  algorithm  to  determine  the  target  SIRs  and  the  numbers  of 
the  two  classes  of  packets  to  be  transmitted  in  each  time  slot.  In  Section  6.4,  we  model  the 
dynamic  behavior  of  the  system  by  a time-homogeneous  discrete-time  Markov  process.  In 
Section  6.5,  performance  of  the  proposed  resource  allocation  algorithm  is  evaluated  using 
both  Monte  Carlo  simulation  and  Markov  analysis.  Also,  we  investigate  the  optimization  of 
the  resource  allocation  algorithm  to  give  the  maximum  possible  throughput.  Comparison 
is  made  between  the  proposed  type-I  hybrid  ARQ  protocol  together  with  the  optimized 
resource  allocation  and  an  equivalent  FDMA  system  that  has  the  same  bandwidth,  the 
same  transmission  power,  and  a similar  ARQ  protocol.  The  proof  of  the  upper  bound  of 
the  retransmission  delay  is  presented  in  Section  6.6.  Finally,  conclusions  are  drawn  in 
Section  6.7. 
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6.1  Data-Link  Layer:  Type-I  Hybrid  ARQ  Protocol 

We  model  the  system  as  a slotted  packet-based  system  over  a point-to-point  link, 
and  employ  the  synchronous  DS-CDMA  system  as  the  underlying  physical  layer  signaling 
technique.  Hence,  multiple  packets  can  be  transmitted  in  a time  slot  due  to  the  multiple 
access  capability  provided  by  the  CDMA  technique. 

At  the  data-link  layer,  the  packet  arrival  process  is  modeled  as  a Poisson  process, 
that  is,  the  number  of  packets  that  arrive  in  a time  slot  is  Poisson  distributed  and  the  num- 
ber of  arrivals  is  independent  from  time  slot  to  time  slot.  The  arrival  rate  of  the  packets  is 
constant  during  the  whole  transmission  period.  We  assume  that  the  transmitter  employs  a 
type-I  hybrid  ARQ  scheme  [74]  with  a buffer  that  can  hold  a maximum  of  B packets.  A 
packet  of  size  L bits  arriving  at  the  buffer  is  admitted  to  the  system  (stored  in  the  buffer)  if 
the  buffer  is  not  full  at  the  time  of  arrival.  Otherwise,  the  packet  is  rejected.  Each  packet  is 
equipped  with  an  error-correcting  code  that  is  capable  of  correcting  e errors  in  the  received 
packet.  For  simplicity,  we  assume  that  the  error-detecting  capability  required  in  the  ARQ 
scheme  is  provided  by  the  error-correcting  code  and  the  receivers  can  always  detect  any 
error  in  a received  packet.  If  there  are  more  than  e errors  in  the  received  packet,  the  receiver 
will  send  a negative  acknowledgment  back  to  the  transmitter  requesting  for  retransmission 
of  the  packet  in  error.  Otherwise,  the  receiver  will  send  a positive  acknowledgment  indi- 
cating that  the  packet  is  successfully  received.  Again  for  simplicity,  we  assume  that  the 
feedback  channel  is  error-free,  the  propagation  delay  is  small,  and  the  acknowledgments 
are  very  short.  Hence  the  transmitter  knows  the  status  of  the  transmitted  packets  almost 
immediately  after  the  transmission  is  completed. 

Packets  in  the  transmitter  buffer  are  transmitted  on  a first-come-first-serve  (FCFS) 
basis.  If  a packet  is  successfully  transmitted  in  the  current  time  slot,  it  will  be  removed 
from  the  buffer,  and  the  packets  behind  will  shift  one  position  forward  in  the  buffer.  If  the 
transmission  of  a packet  fails,  it  will  remain  in  the  buffer  in  the  same  position  relative  to 
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the  other  packets  before  the  transmission1.  In  this  way,  we  can  divide  the  packets  in  the 
buffer  into  two  classes,  namely  the  newly  arrived  packets  and  the  retransmission  packets. 
The  FCFS  criterion  described  above  guarantees  that  the  retransmission  packets  are  always 
in  front  of  the  newly  arrived  packets  in  the  buffer.  Based  on  the  optimal-sequence  CDMA 
signaling  technique  and  the  resource  allocation  algorithm  to  be  described  in  Sections  6.2 
and  6.3,  the  transmitter  transmits  Kr  retransmission  packets  and  Kn  newly  arrived  packets 
in  a time  slot.  Moreover,  different  power  levels  are  assigned  to  the  transmitted  packets  so 
that  the  received  SIR  of  the  Kr  retransmission  packets  and  the  Kn  newly  arrived  packets  are 
7r  and  jn,  respectively.  We  note  that  because  of  the  FCFS  criterion,  all  the  retransmission 
packets  in  the  buffer  are  transmitted  if  Kn  is  larger  than  zero  in  a time  slot. 

6.2  Physical  Layer:  Optimal-Sequence  CDMA 

At  the  physical  layer,  each  packet  to  be  transmitted  is  assigned  a spreading  se- 
quence of  period  N,  where  N is  the  spreading  gain.  The  data  bits  in  a packet  are  BPSK- 
modulated,  and  then  spread  using  assigned  the  spreading  sequence.  The  sum  of  the  spread 
spectrum  signals  carrying  information  of  the  packets  is  transmitted.  Linear  MMSE  re- 
ceivers [28,  30]  are  employed  to  demodulate  the  CDMA  signal. 

6.2.1  Characterization  of  the  Optimal  Sequence 

The  set  of  sequences  (and  implicitly  the  transmission  powers  of  the  packets)  are 
chosen  according  to  the  optimal  scheme  described  in  Lok  and  Wong  [52]  and  Viswanath 
et  al.  [50].  We  assume  that  there  are  K packets  with  target  SIRs  71,72,  • • - ,7 k to  be 
transmitted.  The  underlying  wireless  channel  is  modeled  as  an  AWGN  channel,  and  the 
power  spectral  density  of  the  noise  samples  is  a2  at  the  output  of  the  chip-matched  filter  in 

the  MMSE  receiver  [28,  30],  All  the  K packets  are  admissible  into  the  system  if  and  only 

'This  packet  may  move  forward  in  the  buffer  if  some  of  the  packets  in  front  of  it  are  successfully  trans- 
mitted. 
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if  the  following  capacity  constraint  is  satisfied: 

K 

^e(7 k)<N,  (6.1) 

k=i 

where  e(7*)  = represents  the  effective  bandwidth  [50]  the  &th  packet  takes  with  tar- 
get SIR  7^.  When  constraint  (6.1)  is  satisfied,  there  is  an  optimal  choice  of  sequences  that 
requires  a minimum  amount  of  total  transmission  power  to  support  all  packets  with  target 
SIRs  [50,  52],  For  completeness  of  illustration,  we  give  a brief  introduction  of  the  struc- 
ture of  the  optimal  sequences  for  the  general  case  of  arbitrary  target  SIRs,  which  is  a set  of 
generalized  Welch  bound  equality  (GWBE)  sequence  [50],  This  optimal  sequence  alloca- 
tion scheme  adopts  a user  classification  concept  according  to  their  target  SIR  requirements. 
There  exist  a unique  k*  £ {0, 1,  • • • , N — 1}  such  that 


K 


(■ N - k*)e{ 7fc.)  > e(7j)  ^ (N  ~ k*)e(7fc*+i)> 

j=k*+ 1 


(6.2) 


with  an  assumption  of  a descending  order  of  the  target  SIRs  71  > 72  > • • • > 7 k-  With  this 
constraint,  packets  1, 2,  ■ • • , k*  are  denoted  as  oversized  packets  and  packets  k*  + l,k*  + 
2,  • • • , K are  denoted  as  nonoversized  packets.  Under  the  GWBE  sequence  allocation 
scheme,  oversized  packets  are  allocated  independent  channels  with  transmission  power 
pk  — (727*.,  for  k = 1, 2,  • • ■ , k*,  and  nonoversized  packets  share  the  remaining  channels 


with  transmission  power  pk 


(N-k*)a2e{  7fc) 
^-fc*-Ei=fc*  + ie(7() 


, for  k = k*  + 1,  k*  + 2,  • • • , K.  We 


denote  the  spreading  sequence  of  the  kth  packet  as  an  iV  x 1 column  vector  ak,  and  form 
an  N x K matrix  AT  = [aia2  • • - a*-]  by  grouping  all  spreading  sequences  vectors.  The 
GWBE  sequences  are  completely  characterized  by  a specific  eigenvalue  distribution  of 


ArA": 


(P!,-"Pfc*,  N_k>  N_k*  h 


(6.3) 


given  the  power  allocation  described  above.  The  optimal  sequences  reduce  to  a set  of 
Welch  bound  equality  (WBE)  sequence  of  equal  power  [52]  when  all  packets  require  an 
uniform  target  SIR,  and  reduces  to  a set  of  orthogonal  sequences  when  K < N.  Several 
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attractive  features  of  the  GWBE  sequences  make  it  feasible  for  implementation  in  practical 
system.  The  first  advantage  is  that  it  can  be  easily  constructed  using  iterative  algorithm 
[50].  Another  important  property  is  that  the  MMSE  receivers  reduce  to  the  matched  filter 
receivers  for  all  packets  when  the  GWBE  sequences  are  employed.  This  property  makes 
implementation  complexity  of  using  GWBE  sequences  the  same  as  the  current  practical 
system  provided  that  the  set  of  sequences  currently  employed  are  known  at  the  receiver. 
6.2.2  Optimal  Sequences  with  Two  Classes  of  Packets 


Packets  in  the  buffer  are  divided  into  two  classes  as  illustrated  in  the  data-link  layer. 
A set  of  GWBE  sequences  is  chosen  such  that  the  target  SIRs,  and  7„,  of  the  retrans- 
mission and  the  newly  arrived  packets  are  achieved.  We  assume  that  both  the  transmitter 
and  the  receiver  are  equipped  with  the  same  algorithm  to  determine  the  optimal  sequences. 
At  the  beginning  of  each  time  slot,  the  transmitter  sends  a header  to  the  receiver  to  pass 
the  information  needed,  e.g.,  Kn , Kr , 7n,  jr,  to  generate  the  set  of  optimal  sequences  for 
this  time  slot.  We  assume  that  this  header  is  very  short  and  can  always  be  successfully 
received.  Hence,  its  effect  will  be  ignored  in  the  following  development.  Suppose  that  the 
maximum  total  transmission  power  is  fixed  to  P.  Let  the  ratio  of  the  total  signal  power  to 
the  noise  power  be 


5 = 


P 

No*' 


(6.4) 


The  following  constraints  are  necessary  and  sufficient  for  the  existence  of  the  set  of  optimal 
sequences  as  described  under  different  conditions  based  on  the  results  in  Lok  and  Wong 
[52]  and  Viswanath  et  al.  [50]: 


(i)  for  Kr  + Kn  < N, 


KrJr  + I<n 7„  < NS] 


(6.5) 


(ii)  for  Kr  < N,  Kr  + Kn  > N,  and  ( N - Kr)e{ jr)  > Kne( 7n), 

(N  - Kr)Kne(ln) 


Kr^r  + 


N — Kr  — Kne(  7„) 


< NS] 


(6.6) 
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(iii)  for  Kn  < N,  Kr  + Kn  > N,  and  (N  - Kn)e( 7„)  > Kre( 7r), 

(N  - Kn)Kre( 7r) 


N-a^fn 


N — Kn  — Kre{lr) 


< NS- 


(6.7) 


(iv)  for  Kr  + Kn  > N and  none  of  the  above, 


Kre{ 7r)  + 7f„e(7n)  < Ne(S). 


(6.8) 


The  transmitter  determines  suitable  values  of  Kr,  Kn,  jr,  and  7„  in  each  time  slot 
so  that  the  constraints  (6.5)— (6.8)  are  satisfied.  The  choices  of  the  parameters  are  obtained 
by  a resource  allocation  algorithm  that  maximizes  the  throughput2  of  the  link  while  trying 
to  keep  the  delay  experienced  by  each  packet  within  a reasonable  value. 


6.3  Resource  Allocation  Algorithm 


The  resource  allocation  algorithm  is  responsible  for  selecting  Kr,  Kn,  jr,  and 
in  each  time  slot,  so  that  the  link  throughput  is  maximized  while  the  average  retransmis- 
sion delay3  is  guaranteed  to  be  smaller  than  a predetermined  value  Dr.  The  design  of  the 
resource  allocation  scheme  is  significantly  simplified  by  considering  the  retransmission 
delay  instead  of  the  total  delay. 

6.3.1  Selection  of  Kr  and  -yr 


We  restrict  the  choice  of  the  SIR  of  the  newly  arrived  packets  7„  to  be  lower 
bounded  by  %.  From  constraints  (6.5)— (6.8),  the  maximum  number  of  newly  arrived  pack- 
ets, K (7„),  that  can  be  transmitted  in  a time  slot  is 


^ f [Ne(S)/e(%)\  if  0 < % < S 

K [ln)  - \ [NS/%\  if  7»  > 


(6.9) 


given  7„  and  S.  To  guarantee  that  the  average  retransmission  delay  is  smaller  than  the 

predetermined  value  Dr,  the  resource  allocation  algorithm  assigns  a higher  priority  to  the 

2Here  throughput  is  defined  as  the  average  number  of  successfully  transmitted  packets  per  time  slot. 

Retransmission  delay  is  defined  as  the  number  of  time  slots  elapsed  from  the  moment  that  a newly 
arrived  packet  fails  its  first  transmission  to  the  moment  that  it  is  successfully  transmitted. 
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retransmission  packets  in  the  sense  that  it  allocates  system  resource  to  the  retransmission 
packets  before  considering  the  newly  arrived  packets.  Moreover,  we  impose  the  restriction 
that  the  SIR  of  the  retransmission  packets  satisfies  7r  > %.  We  recall  that  the  constraints 
(6.5)— (6.8)  have  to  be  satisfied  when  the  total  transmission  power  is  limited  to  P.  Hence 
it  is  possible  that  there  are  no  newly  arrived  packets  being  transmitted  ( Kn  = 0)  and  all 
the  transmitted  packets  are  retransmission  packets  in  some  time  slots.  From  the  FCFS 
principle,  we  know  that  a retransmission  packet  will  remain  in  the  buffer  and  keep  being 
retransmitted  every  of  the  following  time  slots  until  the  transmission  succeeds.  In  order  to 
guarantee  that  the  average  retransmission  delay  is  smaller  than  Dr,  we  restrict  the  choice 
of  SIR  of  the  retransmission  packets  7r  as  follows.  Given  jr,  the  maximum  number  of 
retransmission  packets  that  can  be  transmitted  in  a time  slot  is 


lNeiS)/e(lr)\  if  0 < 7r  < S 
l7rj-t  [NS/ 7rJ  if  7 r>S. 


(6.10) 


With  the  lower  bound  jn  on  the  SIR  of  the  newly  arrived  packets,  we  can  obtain  an  upper 
bound  on  the  average  retransmission  delay  experienced  by  a packet  by  considering  the 
following  worst  case  scenario: 


(i)  The  SIR  of  the  newly  arrived  packets  is  set  to  %. 

(ii)  The  transmitter  transmits  K( %)  newly  arrived  packets  and  all  the  K( %)  packets 
fail.  Then  there  will  be  K( %)  retransmission  packets  in  the  buffer  after  the  trans- 
mission. It  is  obvious  that  K (7„)  is  the  maximum  possible  number  of  retransmission 
packets  in  the  buffer. 

(iii)  The  packet  considered  is  the  last  retransmission  packet  in  the  buffer. 

Since  at  most  K{ jr)  retransmission  packets  can  be  transmitted  in  a time  slot,  the  last  re- 
transmission packet  in  the  buffer  will  not  be  retransmitted  until  the  first  (K(%)  — K( jr)) 
packets  are  successfully  transmitted.  Once  the  retransmission  process  of  this  packet  starts, 
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it  will  be  retransmitted  in  every  one  of  the  following  time  slots  until  the  transmission  suc- 
ceeds. Hence  the  retransmission  delay  is  the  sum  of  the  waiting  time  and  the  time  needed 
to  successfully  transmit  the  packet.  The  average  retransmission  delay  of  this  worst  case 
packet  (hence  an  upper  bound  on  the  average  retransmission  delay)  is  given,  as  a function 


of  7 r and  7„,  by 


Dr  (7 rt  Tn)  — < 


P»( 7r) 

B 

KMPsilr) 

Mjn) 

KMPsM 


if  B < Krfrr)  < K{%) 
if  A'(7r)  < B < K{%) 
if  K{yr)  < K{%)  < B, 


(6.11) 


where  Ps(7r)  is  the  probability  of  successfully  transmitting  a retransmission  packet  given 


that  the  SIR  is  set  to  7 r.  With  the  assumption  that  the  MAI  component  at  the  output  of  the 


MMSE  receiver  is  Gaussian  [75]  and  the  bit  errors  are  independent  within  a packet,  the 


packet  success  probability  with  the  target  SIR  7 is  given  by 


p.(  7)  = E (L)  w>M]ii  - nwi1'*,  (6.12) 

where  Pb( 7)  = Q{^/ 7)  is  the  bit  eiror  probability.  To  obtain  (6.11),  we  have  also  made 
the  assumption  that  the  error  events  are  independent  from  packet  to  packet.  To  illustrate 
the  result  above,  we  consider  a sample  system  with  L = 1023,  e = 5,  N = 15,  S = 9dB, 
and  % = 3dB.  Plots  of  the  upper  bounds  on  the  average  retransmission  delay  Dr( jr,%) 
versus  the  choice  of  7r  are  given  in  Fig.  6.1  for  B = 15,  30,  and  60,  respectively.  We  note 
that  the  curves  for  the  latter  two  cases  overlap.  Using  the  bound  in  (6.11),  the  resource 
allocation  algorithm  chooses  7r  to  satisfy  the  requirement  that  the  average  retransmission 
delay  is  bounded  above  by  Dr.  For  example,  if  Dr  is  set  to  4 slots,  we  should  choose 
7.86  dB  < 7,.  < 13.78  dB  to  meet  retransmission  delay  requirement  when  B — 30  or  60 
from  the  observation  to  Fig.  6.1.  We  notice  that  it  may  not  always  be  possible  to  satisfy 
an  arbitrary  retransmission  delay  requirement.  For  example,  it  is  indicated  in  Fig.  6.1  that 
Dr  will  be  larger  than  1.32  slots  for  the  system  considered  no  matter  how  we  choose  jr. 
The  purpose  of  resource  allocation  algorithm  is  to  choose  the  value  of  7r  to  maximize  the 
link  throughput.  We  are  going  to  analyze  the  effect  of  different  choices  of  % and  7,.  on  the 
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Figure  6.1:  Upper  bound  on  average  retransmission  delay 

system  performance  in  Section  6.5.  For  the  moment,  we  assume  that  the  exact  choices  of 
7„  and  have  been  determined  such  that  the  retransmission  delay  requirement  is  satisfied. 
Once  7r  is  fixed,  the  choice  of  Kr  becomes  explicit  (because  of  the  FCFS  transmission 
criterion): 

Kr  = min{iF(7r),  kr},  (6.13) 

where  kr  is  the  number  of  retransmission  packets  currently  stored  in  the  buffer,  that  is,  we 
should  transmit  as  many  retransmission  packets  as  allowed  by  the  power  constraint. 

6.3.2  Selection  of  Kn  and  7n 

After  Kr  and  7r  are  chosen  in  a time  slot,  any  “left-over”  capacity,  as  specified  by 
the  constraints  (6.5)-(6.8),  can  be  used  to  transmit  the  newly  arrived  packets.  The  values 
of  Kn  and  7„  are  chosen  to  maximize  the  throughput  of  the  link. 

For  a fixed  Kn , the  average  number  of  successfully  transmitted  newly  arrived  pack- 
ets in  the  current  time  slot  is  KnPs( jn).  Since  the  packet  success  probability  Ps( jn)  is  a 
monotone  increasing  function  of  7„,  the  throughput  of  the  link  will  be  maximized  if  we 
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choose  7„  to  be  the  largest  possible  value  satisfying  the  constraints  (6.5)-(6.8).  It  can  be 
shown  that  such  a choice  of  7„  (as  a function  of  Kn ) is  given  as  below: 


(i)  for  Kr  < N and  1 < Kn  < N — Kr, 


7 n{Kn)  = - KrlrY, 


(ii)  for  Kr  < N,Kn>  N - Kr  + 1,  and  e(7r)  < ^f*e(S), 


7 n{Kn)  — ts 

n. 


NS 


(N  - Kn)Kre(lr) 
N-Kn-Kre{lr)\  ’ 


(iii)  for  Kr  < N,  Kn  > N - Kr  + 1,  and  e{S ) < e(7r)  < e{S), 


7n(-^n) 


Ne(S)  — Kre(  7r) 


Kn  - Ne(S ) + Kre{^r)  ’ 

(iv)  for  Kr  < N,  Kn  > N — Kr  + 1,  and  e(7r)  > e(S), 

1 

7n(-^n) 


K, 


( i 1 1 ) - 1 ’ 

n \ N-Kr  ^ NS-KrJrJ 


(v)  for  Kr  > N and  e(7r)  < K ne(S), 


7 n{Nn)  — __ 


(VS 


(iV  - X„)AVe(7r) 

N — Kn  — Kre(ir)\  ’ 


(vi)  for  Kr  > Af  and  e(7r)  > ^*e(S), 


7™  ( Kn ) 


Ne(S)  — Kre{^r) 

Kn  - Ne(S)  + Kre{ir)' 


(6.14) 


(6.15) 


(6.16) 


(6.17) 


(6.18) 


(6.19) 


Given  and  jr,  the  optimal  number  of  newly  anived  packets  to  be  transmitted 
Kn,  in  the  sense  of  maximizing  the  throughput,  is  given  by 


Kn  = argrnax  KPS(  jn(K)),  (6.20) 

0 <K<kn;  7n(^0>7n. 

where  7 n{K)  is  defined  in  (6. 14)— (6. 19)  and  kn  is  the  current  number  of  newly  anived 
packets  in  the  buffer.  Once  Kn  is  chosen,  the  target  SIR  of  the  newly  arrived  packets  is  set 
to  be 


7 n — 7n  ( ) • 


(6.21) 
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6.4  Performance  Analysis 
6.4.1  Two-dimensional  Markov  Model 

The  numbers  of  newly  arrived  packets  and  retransmission  packets  and  the  SIR  of 
each  packet  vary  from  time  slot  to  time  slot  due  to  the  statistical  nature  of  transmission  er- 
rors and  retransmission  events.  Since  transmission  decision  in  the  next  time  slot  made  by 
the  ARQ  protocol  is  based  only  on  the  packets  currently  stored  in  the  buffer,  the  dynamic 
behavior  of  the  system  can  be  modeled  by  a time-homogeneous  discrete-time  Markov  pro- 
cess. Each  buffer  slot  can  be  in  one  of  three  different  modes.  It  can  be  either  occupied  by  a 
newly  arrived  packet,  occupied  by  a retransmission  packet,  or  empty.  Let  kn,  kr  and  ke  be 
the  number  of  newly  arrived  packets,  the  number  of  retransmission  packets,  and  the  num- 
ber of  empty  slots  in  the  buffer,  respectively.  We  note  that  0 < kn  < B,  0 < kr  < B — kn, 
and  ke  = B — kn  — kr.  To  model  the  dynamics  of  the  buffer,  it  suffices  to  define  a two- 
dimensional  state  vector  u — (kn,kr)  because  only  two  of  the  three  values  above  are 
independent.  It  is  easy  to  see  that  the  system  has  (g+1>jg+2)  different  states  in  total.  We 
denote  the  state  transition  probability  that  the  system  moves  from  one  state  ut  = {kn\,  kr\ ) 
at  fth  time  slot  to  the  new  state  u>t+ 1 = (kn 2,  K2)  at  (t  + l)th  time  slot  as 

Pt(knu  kri]  kn2i  kr 2)  — Pr  (u7_|-i  — (kn 21  kr2)\ujt  ( kn\,  Ayi))  • (6.22) 

The  two-dimensional  time-homogeneous  Markov  chain  u gives  a sufficient  description  of 
the  dynamics  of  the  system  for  the  purposes  in  this  chapter.  For  instance,  we  can  evalu- 
ate the  stationary  distribution  of  the  Markov  chain,  and  then  obtain  the  steady  state  link 
throughput,  the  transmission  delay  and  the  packet  rejection  probability  of  the  system.  Sup- 
pose that  the  state  of  the  system  at  the  beginning  of  the  ft h time  slot  is  uit  — (kni,  kr\). 
The  numbers  of  the  two  types  of  packets  to  be  transmitted  in  the  fth  time  slot  are  functions 
of  B , N,  S,  7„,  7r,  kni,  and  kr\.  For  notational  simplicity,  we  can  denote  the  number 
of  newly  arrived  packets  and  the  number  of  retransmission  packets  to  be  transmitted,  re- 
spectively, as  Kn(kn  1,  kr  1)  and  Kr(kn l5  kr\)  since  all  other  parameters  are  fixed  during  the 
transmission  process.  The  resource  allocation  algorithm  decides  the  values  of  Kn(kn\ , kr\) 
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and  Kr(kn\,  kr i),  and  the  system  transmits  the  corresponding  packets  at  the  beginning  of 
the  fth  time  slot.  The  new  state  of  the  system  ut+i  for  the  (t  + l)th  time  slot  is  evaluated 
at  the  end  of  the  ith  time  slot.  State  transition  event  of  the  buffer  from  one  time  slot  to  the 
next  time  slot  consists  of  two  steps,  namely,  packet  transmission  and  packet  arrival.  As  a 
result,  construction  of  the  state  transition  matrix  can  also  be  decomposed  into  two  steps. 

Step  1.  At  the  beginning  of  the  fth  time  slot,  Kn(kn i,  kr i)  newly  arrived  packets 
with  SIR  7 n and  Kr(kni,  kr i)  retransmission  packets  with  SIR  7r  are  transmitted.  Packets 
that  are  successfully  transmitted  are  cleared  from  the  buffer,  while  the  failed  packets  stay 
in  the  buffer.  A newly  arrived  packet  will  become  a retransmission  packet  after  its  first 
failure  and  a retransmission  packet  will  continue  to  be  a retransmission  packet  if  it  fails 
again.  We  can  get  the  packet  transmission  success  probabilities  Ps^n)  of  the  newly  arrived 
packets  and  Ps(7r)  of  the  retransmission  packets  from  (6.12)  with  the  knowledge  of  the 
target  SIRs  7n  and  7r,  the  packet  size  of  L bits,  and  the  error  correction  ability  of  e bits. 
Suppose  that  in  out  of  the  Kn(kni,  kr i)  newly  arrived  packets  succeed,  and  ir  out  of  the 
Kr(kn\,  kr\)  retransmission  packets  succeed  during  the  transmission  in  the  ith  time  slot. 
The  transient  state  of  the  buffer  after  the  packet  transmission  is  cu<  = (kn i,  kr\)  with  kn\  = 
kn i - Kn(kn i,  kr i)  and  krl  = kr\  — ir  + Kn(kni,  kr i)  - in.  With  the  assumption  that  the 
transmission  of  each  packet  is  independent,  the  probability  that  the  system  moves  from  the 
original  state  (kn i,  A:ri)  to  the  transient  state  (kn i,  kr i)  is 

P\  {kn 1 , krl  i kn i , ~kri)  = B(Kn(kni,  kri),in,  Ps('1n))B(Kr(kni,kri),ir,  Ps(ir)),  (6.23) 

where  B(K,  i,p)  denotes  the  binomial  distribution,  that  is,  for  0 < i < K,  and  0 < p < 1, 

B{K,i,p)=  (^p^l-p)^.  (6.24) 

We  amange  the  (B+1Hg+2)  states  into  a row  vector 

[(0,  0),  • • • , (0,  B),  (1,  0),  ■ • • , (1,  B — 1),  • • • , {B,  0)].  (6.25) 

Then  the  element  P\(kni,  A;rl ; kn i,  kr\ ) of  the  state  transmission  matrix  Pi,  with  respect  to 
the  ordering  in  (6.25),  can  be  obtained  according  to  (6.23). 
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Step  2.  We  assume  that  the  number  of  packets  that  arrive  the  system  during  the 
tth  time  slot  is  Poisson  distributed.  When  a newly  arrived  packet  arrives,  it  is  admitted  into 
the  system  if  the  buffer  is  not  full.  Otherwise,  it  is  rejected.  If  we  assume  that  jn  newly 
arrived  packets  are  admitted  into  the  buffer,  the  final  state  of  buffer  at  the  end  of  the  tth 
time  slot  is  cot+ 1 = [kn 2,  kr2),  where  kn2  = kn\  + jn  and  kr2  = kr\-  The  state  transition 
probability  of  this  event  can  be  obtained  through  the  Poisson  arrival  distribution: 


f e~x\^ 

in'- 

if 

0 ^ in  < P kn  1 

kr  1 

P2ikv.i1  krh  kn 2,  kr 2)  S 

' 1 ^B-kn\-kr\-\  e~x \in 

I 1 L>jn=  0 jn\ 

if 

jn  — B kn  1 kr\ 

1 

[0 

if 

jn  > B kn\  kr\. 

(6.26) 

Using  (6.26),  we  can  obtain  the  elements  P2(kniikri',  kn2,  kr2)  of  the  state  transmission 
matrix  P2,  with  respect  to  the  ordering  in  (6.25). 

The  state  transition  probability  Pt{kn i,  kr i;  kn2,  kr2 ) from  the  state  c ot  = (kn i,  kr\) 
to  the  state  u)t+\  = {kn2,  kr2)  can  be  obtained  via  a combination  of  Steps  1 and  2. 

B B—kni 

Pt{kn\i  kr i,  kn2,  kr 2)  E E P 1 {kn  1,  kr  1,  kn  1,  kr\)P2(kn\  , kr\ , kn2t  kr 2 )■  (6.27) 

kn  1—0  kr  1 — 0 

Here,  Pt{kn  1,  kr  1;  kn2 , kr2)  is  the  element  of  state  transition  matrix  with  respect  to  the 
ordering  in  (6.25).  It  is  obvious  that  Pt  satisfies 

P(  = PiP^  (6.28) 

Since  the  Markov  chain  ui  is  irreducible,  a stationary  distribution  n exists  [76],  According 
to  the  ordering  in  (6.25), 

7T  = [7r(0,  0),  • • • , 7r(0,  H),7t(1,  0),  • • • , 7r(l,  — 1),  • • • , 7r(5, 0)],  (6.29) 

where  7r(£:n,  kr)  is  the  stationary  probability  of  the  state  ( kni  kr ).  Moreover,  7 r can  be  found 
by  solving  the  following  set  of  equations: 

f 7T  = 7rPj 

l Zlfen=0  ^2kr=0  kr)  — 1. 


(6.30) 
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6.4.2  Link  Throughput,  Packet  Rejection  Probability  and  Transmission  Delay 

System  throughput,  average  transmission  delay  and  packet  rejection  probability  of 
the  system  with  a fixed  load  and  transmission  power  can  be  obtained  readily  from  the 
stationary  distribution  tv  of  the  Markov  chain  evaluated  above. 

6.4.2.1  Link  Throughput 

Throughput  is  defined  as  the  steady  state  average  number  of  successfully  transmit- 
ted packets  per  time  slot,  which  can  be  expressed  mathematically  as 

B B—kni 

T = y y;  {Kn(kni,  kr\)P3{^n)  + Kr[kn  1,  fcri)fi(7r)Mfcwi,  fcri).  (6.31) 

kjl  1 — 0 kr  1 =0 

The  transient  state  of  the  buffer  is  Cbt  = (fcnl,  kr i)  after  the  packet  transmission.  The  buffer 
admits  at  most  ( B — knX  — kr\)  newly  arrived  packets,  and  rejects  any  extra  packets  if  the 
number  of  newly  arrived  packets  during  the  fth  time  slot  exceeds  (B  — kn\  — kr i).  The 
average  number  of  packets  rejected  at  the  transient  state  (kn i,  kr i)  is  then 

B-knl—kr\  \ 

Nr{knUkri)  = \-(B-kni-krl)+  y (B  - kni  - krl  - w)  — ■ (6.32) 

n= 0 

According  to  the  ordering  in  (6.25),  we  define  a row  vector 

Nr  = [iVr(0, 0),  • • ■ , Nr{0,  B),  Nr{l,  0),  • • • , Nr{l,  B - 1),  • • • , Nr{B,  0)].  (6.33) 

6.4.2.2  Packet  Rejection  Probability 

Packet  rejection  probability  is  defined  as  the  steady  state  percentage  of  the  newly 
arrived  packets  that  are  rejected  from  entering  the  transmitter  buffer.  Numerically,  it  can 
be  calculated  as  the  ratio  of  the  average  number  of  packets  rejected  by  the  system  in  a time 
slot  to  the  packet  arrival  rate  A: 

^ B B—kni  B B—kn\ 

Pr  = A ^ ^ ^{knli  krl)  E E P\  {knit  k rli  kn i,  krl)Nr{knUkn) 

kn  1=0  fcrl=0  knl=0  kr  1=0 

7rPiNrr 


A 


(6.34) 
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6.4.2.3  Average  Total  Transmission  Delay 

Average  total  transmission  delay  is  defined  as  the  average  number  of  time  slots 
elapsed  from  the  moment  that  a newly  arrived  packet  is  admitted  into  the  system  to  the 
moment  that  it  is  successfully  transmitted.  We  can  apply  Little’s  formula  [77]  to  obtain  the 
steady  state  total  transmission  delay  Dt.  According  to  Little’s  formula,  the  average  number 
of  packets  in  the  buffer  is  equal  to  the  product  of  the  average  arrival  rate  and  the  average 
total  transmission  delay.  Taking  the  rejected  packets  into  consideration,  the  steady  state 
total  transmission  delay  is  given  by 

p _ Efcnl=0  Ylkri=0  (^"1  + krl)^(knU  &rl) 

A " (633) 

6.4.2.4  Average  Retransmission  Delay 


Retransmission  delay  is  defined  as  the  number  of  time  slots  elapsed  from  the  mo- 
ment that  a newly  arrived  packet  fails  its  first  transmission  to  the  moment  that  it  is  success- 
fully transmitted.  The  average  retransmission  delay  Dr  can  also  be  obtained  by  Little’s 
formula.  We  note  that  the  arrival  rate  of  the  retransmission  packets  is  the  transmission 
failure  rate  of  the  newly  arrived  packets.  Hence,  it  is  straightforward  that 

EL=o  E?~= V kriK{knUkrl) 


Dr  x ft  ^B-kn 


E*nl= 0 E*rl= o1  n(kni,kri)Kn{knl,krl)(l  - Ps( 7„)) 

6.5  Numerical  Results 


(6.36) 


In  this  section,  we  compare  the  analytical  results  developed  above  with  results  ob- 
tained from  Monte  Carlo  simulations  for  the  proposed  type-I  hybrid  ARQ  protocol  and 
resource  allocation  algorithm  over  the  optimal-sequence  CDMA  link.  We  also  examine 
the  system  performance  for  different  choices  of  the  retransmission  packet  SIR  and  the 
lower  bound  jn  of  the  newly  arrived  packet  SIR  to  further  optimize  the  resource  allocation 
algorithm. 


Ill 


Figure  6.2:  Throughput  of  CDMA  system  at  different  load 


6.5.1  Analytical  Results  vs  Simulation  Results 

We  consider  a sample  system  with  parameters  L — 1023,  e — 5,  N = 15,  B = 30, 
and  S = 9dB.  The  upper  bound  on  the  average  retransmission  delay  is  set  to  Dr  = 4 
slots.  We  set  % to  3dB  and  employ  the  smallest  7r  that  satisfies  the  retransmission  delay 
requirement  (see  Fig.  6.1).  Plots  of  the  link  throughput  T,  packet  rejection  probability 
Pr,  average  total  transmission  delay  Dt,  and  average  retransmission  delay  Dr  as  functions 
of  the  packet  arrival  rate  (load)  A obtained  by  the  Markov  chain  analysis  and  Monte  Carlo 
simulations  are  given  in  Figs.  6.2,  6.3,  6.4,  and  6.5,  respectively.  We  can  see  that  the  results 
obtained  by  the  Markov  chain  analysis  and  the  Monte  Carlo  simulation  match  well. 

We  observe  that  the  maximum  throughput  (link  capacity)  achieved  is  slightly  less 
than  14  packets/slot  from  Fig.  6.2.  we  observe  that  the  throughput  increases  linearly  when 
the  input  load  is  smaller  than  the  link  capacity,  and  finally  levels  off  at  the  maximum  value 
after  the  input  load  exceeds  the  link  capacity.  The  resource  allocation  algorithm  limits 
the  maximum  number  of  packets  transmitted  in  each  time  slot  so  that  the  increased  load 


average  total  delay  (slots)  C package  rejection  probability 
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6.3:  Packet  rejection  probability  of  CDMA  system  at  different  load 


0 5 10  15  20  25  30 

packet  arrival  rate  (packets/slot) 


Figure  6.4:  Average  total  transmission  delay  of  CDMA  system  at  different  load 
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Figure  6.5:  Average  retransmission  delay  of  CDMA  system  at  different  load 

does  not  cause  excessive  congestion  on  the  CDMA  link.  Hence,  maximum  throughput 
is  maintained  even  when  the  load  exceeds  the  link  capacity.  The  difference  between  the 
input  load  and  the  link  capacity  translates  to  the  fact  that  some  newly  arrived  packets 
are  rejected  from  entering  the  transmitter  buffer  as  indicated  in  Fig.  6.3.  It  is  obvious 
that  the  packet  rejection  probability  increases  as  the  input  load  increases.  Moreover,  the 
design  of  the  proposed  resource  allocation  algorithm  is  verified  by  Fig.  6.5  that  the  average 
retransmission  delay  is  smaller  than  the  predetermined  bound  of  4 slots. 

6.5.2  Effect  of  yr  and  % on  System  Performance 

We  know  that  Dr( 7r,  7„)  is  a function  of  and  % from  the  derivation  of  the  upper 
bound  on  the  average  retransmission  delay.  As  discussed  in  Section  6.3,  there  is  a range  of 
values  of  that  satisfy  the  retransmission  delay  requirement.  The  value  of  7„  decides  the 
maximum  number  of  newly  arrived  packets  that  can  be  put  into  transmission  within  a time 
slot,  which  affects  the  system  performance.  A natural  question  is  what  the  best  choice  of 
7r  and  % would  be  in  the  sense  of  maximizing  the  link  throughput  given  a retransmission 
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delay  requirement.  In  this  section,  we  employ  Markov  chain  analysis  to  evaluate  the  system 
performance  under  different  choice  of  and  %. 

6.5.2. 1 Effect  of  7r 

There  is  a range  of  values  of  that  satisfy  the  retransmission  delay  requirement. 
We  search  the  optimal  7*  within  this  range  that  maximizes  the  link  throughput  of  the  sam- 
ple system  discussed  in  the  previous  section,  assuming  a fixed  input  load  of  A = 30  pack- 
ets/slot. We  can  see  from  Fig.  6.1  that  the  search  range  of  is  from  7.86dB  to  13.78dB. 
The  system  performance  results  vs  the  choice  of  -yr  obtained  from  the  Markov  chain  anal- 
ysis are  shown  in  Figs.  6.6-6.9.  We  observe  that  there  exists  an  optimal  value  of  7*  = 9dB 
that  achieves  the  maximum  throughput.  We  note  that  this  choice  of  7,.  also  minimizes  the 
packets  rejection  probability  and  the  total  transmission  delay.  Incidentally,  this  is  also  the 
value  of  7 r that  gives  the  smallest  value  of  the  upper  bound  of  the  average  retransmission 
delay  in  Fig.  6.1.  Fig.  6.9  shows  that  the  average  retransmission  delay  levels  off  at  the 
value  of  1 slot  when  the  jr  is  large.  This  is  reasonable  since  the  transmission  success  prob- 
ability Ps( jr)  approaches  1 when  7r  is  large,  and  hence  the  retransmission  packet  is  always 
successfully  transmitted  in  one  time  slot.  The  existence  of  such  an  optimal  choice  of  jr 
has  been  consistently  observed  for  other  choices  of  the  system  parameters.  To  maximize 
the  throughput,  the  resource  allocation  algorithm  should  set  jr  to  the  value  that  gives  the 
smallest  upper  bound  of  the  average  retransmission  delay. 

6.5.2.2  Effect  of  % 

Selection  of  the  lower  bound  jn  on  the  newly  arrived  packet  SIR  affects  the  re- 
transmission delay  bound  DT( 7r,  7n),  and  hence  affects  the  selection  of  the  retransmission 
packet  SIR  7,..  Here,  we  examine  the  effect  of  % on  the  link  throughput.  We  change  % 
from  ldB  to  lOdB.  For  each  value  of  %,  the  retransmission  packet  SIR  is  set  to  7*,  which 
gives  the  smallest  value  for  the  average  retransmission  delay  bound,  as  well  as  the  maxi- 
mum throughput  as  discussed  above.  It  is  shown  in  Fig.  6. 10  that  the  selection  of  the  lower 
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Figure  6.6:  Throughput  of  CDMA  system  at  different  7 r 


Figure  6.7:  Packet  rejection  probability  of  CDMA  system  at  different  7r 
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6.8:  Average  total  transmission  delay  of  CDMA  system  at  different  7r 


Figure  6.9:  Average  retransmission  delay  of  CDMA  system  at  different  7r 
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Figure  6.10:  Throughput  of  CDMA  system  at  different  7n 


bound  7„  of  the  newly  arrived  packet  SIR  has  minimal  effect  on  the  link  throughput,  as 
long  as  7„  is  not  too  large  (less  than  9dB  for  the  sample  system).  When  jn  is  too  large, 
the  throughput  decreases.  The  reason  is  that  the  resource  allocation  algorithm  requires 
7n  > 7n  and  7r  > Hence,  a large  value  of  limits  the  number  of  packets  that  the 
system  transmits  in  each  time  slot. 

6.5.3  Comparison  with  an  Equivalent  FDMA  System 

It  is  illustrative  to  compare  the  performance  of  the  proposed  type-I  hybrid  ARQ 
protocol  over  the  optimal  CDMA-link  with  other  “equivalent”  systems.  For  this  purpose, 
we  compare  the  performance  of  the  proposed  CDMA  system  to  that  of  a frequency  divi- 
sion multiple  access  (FDMA)  system  with  the  same  system  bandwidth  and  the  same  total 
transmission  power.  Since  the  spreading  gain  of  the  CDMA  system  is  N,  the  “equivalent” 
FDMA  system  would  have  N parallel  narrowband  channels,  each  supporting  the  same  bit 
rate  as  that  supported  by  a single  code  in  the  CDMA  system.  We  assume  that  a standard 
type-I  hybrid  ARQ  protocol  [74]  (at  most  one  packet  is  transmitted  per  slot  per  channel) 
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with  the  same  error  correcting  capability  and  the  same  packet  size  as  the  one  used  in  the 
CDMA  system  is  independently  applied  on  each  of  the  N parallel  narrowband  channels. 
Moreover,  we  assume  the  same  power  spectral  density  of  the  noise  in  the  CDMA  and 
FDMA  systems. 

We  compare  the  maximum  throughputs  of  the  proposed  CDMA  system  and  the 
equivalent  FDMA  system.  For  the  FDMA  system,  we  can  choose  to  evenly  distribute 
the  total  transmission  power  over  a subset  of  the  N channels.  In  this  case,  the  maximum 
throughput  achievable  by  the  FDMA  system  is  upper  bounded  by 

TpoMA  ~ max  KPs(NS/K).  (6.37) 

1 <K<N 

We  choose  a sample  system  with  parameters  L = 1023,  e = 10,  N = 15,  B = 30, 
Dr  = 4 slots  and  % = 3dB.  The  SIR  of  the  retransmission  packets  is  set  as  7*,  which 
gives  the  smallest  value  for  the  average  retransmission  delay  bound  A-(7r>7n)-  The  link 
capacity  of  the  proposed  type-I  ARQ  protocol  over  the  optimal  CDMA  link  versus  the  total- 
signal-power-to-noise  ratio  S is  plotted  in  Fig.  6.1 1.  The  upper  bound  on  the  maximum 
throughput  of  the  equivalent  FDMA  system  is  also  plotted  on  the  same  figure.  We  observe 
that  the  proposed  CDMA  system  achieves  the  same  maximum  throughput  as  the  upper 
bound  of  the  “equivalent”  FDMA  system  when  the  transmission  power  is  limited.  When 
the  transmission  power  is  abundant,  it  is  clear  Fig.  6.1 1 that  the  proposed  CDMA  system 
provides  a significant  capacity  improvement  over  the  FDMA  system.  For  example,  when 
S = 20dB,  the  maximum  throughput  achieved  by  the  proposed  CDMA  system  is  about 
17  packets  per  slot,  which  is  13.3%  higher  than  that  of  the  FDMA  system.  This  is  easy 
to  understand  from  the  different  mechanisms  of  the  FDMA  and  CDMA  systems.  The 
FDMA  system  can  transmit  at  most  N packets  within  one  time  slot  over  its  N narrowband 
channels.  Hence  the  maximum  link  capacity  in  the  FDMA  system  is  N no  matter  how 
large  the  transmission  power  is.  The  CDMA  system  enables  us  to  transmit  more  than 
N packets  within  one  time  slot  simultaneously.  With  the  proposed  resource  allocation 
scheme  and  the  error  correcting  capability  provided  by  the  type-I  hybrid  ARQ  protocol, 


119 


Figure  6.11:  Link  capacity  of  CDMA  system  vs  FDMA  system 

it  is  possible  to  attain  a link  capacity  larger  than  N packets  per  time  slot.  It  is  worth 
to  mention  that  the  transmission  power  in  the  system  is  used  in  the  most  efficient  way 
through  the  employment  of  the  optimal  CDMA  signaling  and  proposed  ARQ  scheme.  The 
system  achieves  its  maximum  possible  capacity  given  the  total  transmission  power. 

6.6  Proof  of  Upper  Bound  of  the  Average  Retransmission  Delay 

In  this  section,  we  show  that  the  function  Dr( ryr,‘ln)  defined  in  (6.11)  gives  an 
upper  bound  on  the  average  retransmission  delay. 

To  guarantee  that  the  average  retransmission  delay  is  smaller  than  the  predeter- 
mined value  Dr,  the  resource  allocation  algorithm  assigns  a higher  priority  to  the  retrans- 
mission packets  in  the  sense  that  it  allocates  system  resource  to  the  retransmission  packets 
before  considering  the  newly  arrived  packets.  The  transmitter  is  allowed  to  transmit  at 
most  K ('7n)  packets  within  one  time  slot,  given  the  lower  bound  7„  on  the  target  SIR  of 
the  newly  arrived  packets.  Similary,  K^r)  is  the  maximum  number  of  retransmission 
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packets  that  can  be  transmitted  within  one  time  slot,  given  that  the  target  SIR  of  retrans- 
mission packets  is  set  to  jr.  We  always  have  k(jr)  < E(jn),  with  the  constrain  jr  > %. 
The  scenario  that  the  transmitter  transmits  the  maximum  number  of  newly  arrived  packets 
and  all  the  K( %)  transmitted  packets  fail  generates  the  maximum  number  of  retransmis- 
sion packets  in  the  buffer.  There  are  three  different  cases  to  consider  in  obtaining  the  upper 
bound  on  the  average  retransmission  delay,  depending  on  the  buffer  length  B,  SIR  lower 
bound  7n  of  newly  arrived  packets  and  the  target  SIR  jr  of  the  retransmission  packets. 

Case  1:  K{ jr)  < K(%)  < B.  Under  the  worst  case  scenario,  there  are  K( %) 
retransmission  packets  in  the  buffer.  To  obtain  the  upper  bound  on  the  average  retransmis- 
sion delay,  we  consider  the  last  retransmission  packet,  that  is,  the  iT(7n)th  retransmission 
packet.  This  packet  has  to  wait  until  the  first  K( %)  - K(jr)  retransmission  packets  are  all 
successfully  transmitted  before  its  retransmission.  Once  this  retransmission  packet  starts  to 
be  retransmitted,  it  will  keep  being  retransmitted  until  the  transmission  succeeds.  Hence, 
the  retransmission  delay  includes  two  parts,  namely,  a waiting  time  Tw  and  a transmission 
time  Ts.  Assume  that  the  transmission  events  of  the  packets  are  independent.  Then  the 
number  of  successfully  retransmitted  packets  in  the  fth  time  slot,  denoted  by  Xt,  is  a bino- 
mial random  variable  ranging  from  0 to  K( jr)  with  mean  E[ Xt\  = K('yr)Ps(,yr).  Also,  Xt 
is  independent  from  time  slot  to  time  slot.  Suppose  that  the  K (7„)th  retransmission  packet 
waits  Tw  time  slots  before  starting  to  be  retransmitted.  Then  it  is  straightforward  that 

Tw 

STw  = = k{ %)  - K( 7r).  (6.38) 

t= l 

Making  use  of  Wald’s  identity  [78],  the  average  waiting  time  is  given  by 

E[T„]  = (6-39) 

Since  the  packet  success  probability  for  the  retransmission  packets  is  Ps( jr),  the  average 
transmission  time  of  a retransmission  packet  can  be  obtained  easily  as  E[TS\  = ■ 

Hence  the  upper  bound  on  the  average  retransmission  delay  in  this  case  can  be  expressed 
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as 


DriPri^ln) 


Case  2:  K( jr)  <B<  K{%). 

B newly  arrived  packets  within  one  time  slot.  Hence,  the  worst  case  scenario  is  that  there 
are  B retransmission  packets  in  the  buffer.  To  obtain  the  upper  bound  on  the  average 
retransmission  delay,  we  consider  the  B th  retransmission  packet.  This  packet  has  to  wait 
until  the  first  B - K( jr)  packets  are  all  successfully  transmitted  before  its  retransmission. 
Similar  to  the  previous  case,  we  have 

St.  = y>,  = B - A’( 7r),  (6.41) 

t= 1 

and  E[Xt\  = k[^r)Ps{^r).  Making  use  of  Wald’s  identity  as  above,  the  upper  bound  on 
the  average  retransmission  delay  in  this  case  can  be  expressed  as 


E[TW\  + E[TS\ 

| i 

K(%) 


(6.40) 


K{jr)Ps{lr) 

In  this  case,  the  transmitter  can  transmit  at  most 


A-(7r,7n)  = E[Tw]  + E[TS] 

B - K{ 7r)  | 1 

K{lr)Ps{lr)  Ps(  7r) 

B 

K(lr)PsM 


(6.42) 


Case  3:  B < K( 7r)  < K( %).  In  this  case,  the  transmitter  can  transmit  at  most 
B newly  arrived  packets  within  one  time  slot.  Hence,  the  worst  case  scenario  is  that  there 
are  B retransmission  packets  in  the  buffer.  To  obtain  the  upper  bound,  we  consider  the  B th 
retransmission  packet.  Because  B < K( jr),  the  Bth  packet  in  the  buffer  is  retransmitted 
in  the  next  time  slot  right  after  its  first  failure  without  waiting.  The  upper  bound  of  the  total 
retransmission  delay  in  this  case  is  just  the  transmission  time  Ts 


A-(7r,7n)  = E[T„ \ 


1 


P.hrY 


(6.43) 
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6.7  Summary 

We  presented  a type-I  hybrid  ARQ  protocol  and  an  accompanying  resource  allo- 
cation algorithm  to  make  use  of  the  flexibility  provided  by  the  optimal-sequence  CDMA 
link.  A two  dimensional  Markovian  model  was  formulated  to  model  the  dynamic  behavior 
of  the  proposed  system.  Based  on  this  model,  Markov  analysis  was  employed  to  evalu- 
ate the  system  performance,  leading  to  closed-form  expressions  for  the  link  throughput, 
packet  rejection  probability,  and  packet  transmission  delay.  Numerical  results  indicate  that 
when  the  transmission  power  is  limited,  the  proposed  CDMA  system  achieves  the  same 
link  capacity  as  an  equivalent  FDMA  system  having  the  same  bandwidth,  the  same  trans- 
mission power,  and  a similar  ARQ  protocol.  When  the  transmission  power  is  abundant, 
the  proposed  system  outperforms  the  FDMA  system  in  terms  of  achieving  a larger  link 
throughput.  In  this  sense,  the  proposed  protocol  makes  efficient  use  of  the  available  trans- 
mission power.  The  power  management  unit  of  a wireless  link  can  utilize  the  proposed 
protocol  to  adaptively  change  the  transmission  power  in  order  to  cater  for  the  time-varying 
demand  on  the  throughput  of  the  link. 


CHAPTER  7 

CONCLUSIONS  AND  FUTURE  WORK 


The  main  focus  of  this  work  is  to  explore  the  potential  of  the  physical  layer  tech- 
nique of  spreading  sequence  optimization  in  wireless  CDMA  communication  systems.  The 
objective  of  the  sequence  optimization  is  to  support  a desired  set  of  users  in  the  system  with 
the  target  SIRs  using  a minimum  amount  of  transmission  power.  This  work  is  prompted 
by  the  increasing  demand  of  supporting  high  data  rates  at  a variety  of  quality-of-service 
(QoS)  levels  for  future  multimedia  communications  when  the  usable  frequency  spectrum 
and  transmission  power  are  limited. 

7.1  Conclusions 

In  order  to  design  power-efficient  and  QoS-assuring  wireless  networks,  advanced 
processing  techniques  need  to  be  developed  at  both  the  physical  layer  and  upper  layers. 
Chapters  3-6  present  the  contributions  of  this  work  at  both  the  physical  layer  and  the  data 
link  layer. 

Chapter  3 starts  with  a brief  summary  of  the  sequence  optimization  problem  in 
CDMA  systems  over  an  AWGN  channel  developed  in  Wong  and  Lok  [46,  52]  and  Viswanath 
[50].  Different  scenarios  of  user  admission/departure  in  practical  systems  are  cast  as  three 
different  optimization  problems,  namely  single-user  optimization,  multi-user  optimization, 
and  all-user  optimization.  Thereafter,  two  joint  transmitter-receiver  adaptation  schemes  are 
proposed  to  utilize  the  optimal  spreading  sequences.  In  both  schemes,  training  symbols 
are  sent  in  blocks  so  that  spreading  sequences  and  MMSE  receiver  weights  are  adapted 
simultaneously.  In  the  first  scheme,  spreading  sequences  are  updated  at  the  receivers  in  a 
distributed  manner,  and  are  fed  back  to  the  base-station  to  spread  the  next  block  of  data. 
In  the  second  scheme,  spreading  sequences  are  updated  at  the  base-station  in  a centralized 
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manner,  using  a minimum  amount  of  feedback  information  from  the  receivers.  Simulation 
results  show  that  the  second  scheme  gives  better  performance  in  terms  of  faster  conver- 
gence and  less  oscillation. 

Chapter  4 focuses  on  power  and  sequence  allocation  in  a CDMA  forward  link.  We 
assume  that  receivers  of  different  users  observe  independent  slow  flat  fading.  Two  allo- 
cation schemes,  a minimum-ETSC  scheme  for  MMSE  and  matched-filter  signal  reception 
and  a minimum-TSC  scheme  for  matched-filter  signal  reception,  are  proposed  to  support 
users  in  the  system  with  a uniform  target  SIR.  Both  theoretical  derivations  and  Lagrangian 
multiplier  search  results  indicate  that  the  minimum-ETSC  scheme  exhibits  higher  power 
efficiency  than  the  minimum-TSC  scheme.  In  addition,  there  is  a strong  indication  that 
the  minimum-ETSC  scheme  may  be  the  optimal  scheme  that  consumes  the  minimum  to- 
tal pow'er  when  the  matched-filter  receiver  is  used.  The  minimum-ETSC  scheme  is  not 
optimal  when  the  MMSE  receiver  is  used.  However,  Lagrangian  multiplier  searching  re- 
sults indicate  that  the  minimum-ETSC  scheme  is  effective  and  consumes  only  slightly 
more  power  than  the  minimum  power  searched  under  most  channel  conditions.  Adap- 
tive algorithms  and  structures  are  also  presented  to  facilitate  the  implementation  of  the 
minimum-ETSC  allocation  scheme  in  real-life  communication  systems. 

Chapter  5 considers  the  problem  of  sequence  optimization  in  CDMA  systems  with 
multipath.  We  simplify  the  general  model  in  Chapter  2 to  a point-to-point  link,  in  which 
all  receivers  are  co-located  and  all  signals  undergo  the  same  multipath  channel.  We  solve 
the  optimization  problem  of  supporting  users  with  a uniform  target  SIR  using  the  minimal 
total  power.  Eigenvectors  of  the  path  matrix  are  treated  as  the  channel  space  stated  above 
spanned  by  the  system.  The  optimal  solution  distributes  users’  spreading  sequences  into 
the  channel  space  in  the  optimal  manner.  Comparisons  of  optimal  solution  to  a suboptimal 
approach,  which  uses  the  eigenvectors  of  the  path  matrix  directly  as  orthogonal  spreading 
sequence,  are  made  using  computer  simulation. 
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Chapter  6 aims  to  design  a data  link  protocol  to  allow  simultaneous  transmission 
of  multiple  packets  over  a CDMA  link  using  the  optimal  spreading  sequences.  Under  the 
proposed  type-I  hybrid  ARQ  protocol,  transmitted  packets  are  classified  into  two  classes, 
namely  newly  arrived  packets  and  retransmission  packets.  A resource  allocation  algorithm 
is  constructed  to  determine  the  numbers  and  the  target  SIRs  for  each  class  of  packets  to 
be  transmitted  in  each  time  slot  to  maximize  the  system  throughput  under  a delay  con- 
straint. We  model  the  dynamic  behavior  of  the  system  by  a time-homogeneous  discrete- 
time Markov  process  and  employ  this  Markovian  model  to  analyze  the  performance  of  the 
system.  Numerical  results  obtained  from  both  the  Markov  analysis  and  the  Monte  Carlo 
simulation  show  that  the  proposed  type-I  hybrid  ARQ  protocol  with  CDMA  signaling  out- 
performs an  equivalent  FDMA  system  that  uses  the  same  bandwidth,  the  same  transmission 
power,  and  a similar  ARQ  protocol,  particularly  when  the  power  is  abundant. 

7.2  Future  Work 

Spreading  sequence  optimization  in  CDMA  systems  over  an  AWGN  channel  is 
thoroughly  analyzed  and  understood.  However,  most  of  wireless  communication  systems 
work  over  more  complex  channels.  In  this  work,  we  have  attempted  to  address  spreading 
sequence  design  in  a CDMA  forward  link  over  slow  flat  fading,  and  proposed  a minimum- 
ETSC  scheme  that  has  good  power  efficiency  in  most  channel  conditions.  However,  more 
research  needs  to  be  done  to  obtain  the  structure  of  the  optimal  spreading  sequences  that 
require  the  minimum  transmission  power. 

We  have  also  investigated  the  sequence  optimization  problem  for  CDMA  systems 
with  multipath,  and  solved  the  optimization  problem  for  the  case  of  point-to-point  trans- 
mission. However,  inter-symbol  interference  (ISI),  which  is  omitted  in  the  signal  model 
in  Chapter  2,  will  compromise  the  receiver  output  SIR  when  the  optimal  sequences  are 
employed.  One  of  the  standard  methods  of  dealing  with  ISI  is  to  insert  a cyclic  prefix  in 
signals.  The  cyclic  prefix  creates  a guard  space  between  neighboring  transmit  blocks  in 
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the  time  domain  so  that  ISI  can  be  completely  eliminated.  In  this  application  scenario, 
the  cyclic  prefix  elements  are  taken  from  the  last  (L  — 1)  elements  of  the  original  spread- 
ing sequence  a^,  and  are  copied  to  the  beginning  of  it  to  form  a new  spreading  sequences 
ak  = [a^_L+1, . . . , a^_v  a(0k\  a[k\  . . . , a^1]T.  For  convenience,  we  use  a D transform 
to  denote  the  Ith  spread  spectrum  symbol  of  the  kth  user 

sk(D)  = 6ifc)(a^L+1JD°  + ...  + a^lJD-i+2  + a[)fc)D-z'+1  + ...  + aS;L122-L-;v+2),  (7.1) 

where  D represents  a delay  of  one  chip.  The  multipath  channel  response  can  also  be 
formulated  in  the  D transform  as  hk{D ) = hk,\D°  + hkpD~l  + • • • + hk,L,D~L+l . The 
received  signal  is  passed  through  a chip  matched  filter  and  the  filter  output  is  sampled  every 
chip  interval.  Due  to  multipath  convolution,  there  are  (N  + 2L  — 2)  samples  for  every  data 
symbol,  which  overlap  with  the  samples  of  the  neighboring  symbols  by  (L  — 1)  samples. 
Overlapped  samples  are  discarded  at  the  receiver,  while  the  N samples  in  between,  which 
we  arrange  into  an  N x 1 column  vector  rfc,  are  fed  into  an  N-tap  MMSE  receiver  to  form 
decision  statistics.  Similar  to  the  signal  model  in  Chapter  2,  rk  can  also  be  expressed  as  a 
multiplication  of  the  channel  matrix  with  the  spreading  sequences: 

K 

ffc  = Hfc^6i(m)am-FnA:,  (7.2) 

m—  1 

where  b and  am  are  the  transmitted  bit  and  the  spreading  sequence  of  the  mth  user, 
respectively, 
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and  hk  is  an  AWGN  noise  vector  with  elementwise  variance  r]k.  We  note  that  channel 
matrice  Hfe,  for  A:  = 1, 2,  • • • , K,  are  circulant  matrices.  It  is  well  known  that  any  circu- 
lant  matrix  can  be  diagonalized  by  the  discrete  Fourier  transform  (DFT)  matrix,  which  is 
defined  as 

F = [uij  = — (7.4) 

where  u = . Diagonlizing  the  channel  matrices  with  the  DFT  transform  = 

FAkFH,  and  letting  a.k  = FHak,  we  can  rewrite  the  sequence  optimization  problem  as 

K 

min  ||afc||2 

k= 1 


subject  to 


■ h 


' I< 

E 


-i  -i 


Vk(FHF) 


v-l 


h > e(7), 


(7.5) 


_m= 1 

for  1 < k < K.  We  note  that  the  effect  of  multipath  channel  is  transformed  to  additive 
colored  noise  matrices  in  (7.5).  Hence,  the  original  sequence  optimization  is  simplified 
through  the  utilization  of  cyclic  prefix.  For  the  case  of  point-to-point  transmission,  optimal 
solution  can  be  readily  obtained  using  the  same  technique  as  was  adopted  in  Chapter  5.  For 
the  case  of  CDMA  forward-link  transmission,  this  simplified  formulation  makes  it  easier 
for  us  to  attack  the  problem  than  the  original  one. 

In  DS-CDMA  communication  systems,  spreading  sequences  of  all  users  are  de- 
sired to  have  small  cross-correlations.  The  total  squared  correlation  (TSC),  which  indi- 
cates the  mutual  interference  among  the  spreading  sequences  of  all  users,  forms  a gauge 
of  performance  for  the  spreading  sequences.  Optimality  of  the  WBE  sequences  and  the 
GWBE  sequences  in  CDMA  systems  over  an  AWGN  channel  is  proved  by  Massy  [47]  and 
Viswanath  et  al.  [50].  An  important  feature  of  the  classes  of  TSC  minimization  sequences 
is  that  they  can  be  constructed  iteratively  using  the  TSC  reduction  algorithms.  In  CDMA 
systems  over  fading  channels,  construction  of  the  optimal  sequences  may  be  difficult  due 
to  the  complexity  of  the  model.  Instead  of  obtaining  the  exact  optimal  solution,  we  can 
design  “good”  spreading  sequence  sets  that  consume  only  slightly  more  power  than  the 
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minimum  power  required.  The  minimum-ETSC  scheme  and  minimum-TSC  proposed  in 
Chapter  4 are  examples  of  this  design  approach.  Similar  approaches  may  also  be  used  in 
CDMA  systems  with  multipath  by  choosing  a set  of  spreading  sequences  that  minimizes 
a variant  expression  of  TSC  among  transformed  sequences.  Numerical  method,  such  as 
Lagrangian  multiplier  search,  can  be  developed  to  solve  the  optimization  approximately, 
and  be  used  to  evaluate  the  resulting  total  transmission  power  required  by  the  designed 
sequences. 

In  the  type-I  ARQ  protocol  proposed  in  Chapter  6,  we  only  consider  one  type  of 
traffic.  We  model  the  system  by  a time-homogenous  discrete-time  Markov  process  and 
obtain  a stationary  probability  distribution  by  calculating  the  state  transition  probabilities. 
The  complexity  of  the  system  is  reasonable  for  us  to  analyze  the  system  performance  ac- 
curately using  the  method  of  Markov  analysis.  Third-generation  wireless  communication 
systems  support  heterogeneous  traffics,  consisting  of  voice,  video,  data,  and  etc.  There  are 
different  QoS  requirements  for  different  services.  A future  research  direction  is  to  develop 
data-link  layer  protocols  that  support  multimedia  communications.  The  system  can  still  be 
modeled  exactly  by  a multi-dimensional  Markov  process  in  the  existence  of  multiple  ser- 
vices. However,  the  model  will  be  too  complex  to  allow  us  to  obtain  a stationary  probability 
distribution  and  to  study  the  system  performance  accurately.  This  analytical  difficulty  can 
be  circumvented  by  utilizing  the  equilibrium  point  analysis  (EPA)  technique  that  makes 
use  of  the  equilibrium  point  to  compute  the  system  performance  approximately.  An  equi- 
librium point  is  defined  as  a point  at  which  the  expected  change  on  modes  of  the  system 
is  zero.  EPA  assumes  that  the  system  is  always  working  at  an  equilibrium  point  when  it 
reaches  the  steady  state.  We  can  easily  analyze  the  multi-dimensional  Markov  chains  by 
the  EPA  technique  since  it  does  require  calculation  of  the  state  transition  probabilities. 
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